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FOREWORD 

This report was prepared by personnel of the Flight Mechanics 
Division, Air Force Flight Dynamics Laboratory, and the Digital 
Computation Division, Aeronautical Systems Division. The report was 
prepared under Project ,431, "Flight Path Analysis," Task 143109, 
"Trajectory and Motion Analysis of Flight Vehicles." The formulation 
and interim documentation were completed by Major Urban H. D. Lynch. 
Programming was accomplished by Mr, Fay 0. Young of the Digital 
Computation Division (ASVCP) of the Aeronautical Systems Division 
Computer Science Center (ASV), This report, prepared by Mr. John 
J. Dueweke of the High Speed Aero Performance Branch ( FXG ), combines 
the applicable portions of FQL-TDR-64-1, Part I, Volume 1, with the 
interim documentation prepared by Major Lynch. 

This report is divided into four parts: 

Part I: Capabilities of the Takeoff and Landing Analysis Computer 
Program 

Part II: Problem Formulation 

Pa-'t III: User 1 s Manual 

Part IV: Programmers Manual 

This document was submitted by the authors September 1971, 

This report has been reviewed and is approved. 

PHILIP P. ANT0NAT0S 

Chief, Flight Mechanics Division 

Air Force Flight Dynamics Laboratory 
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ABSTRACT 

A well-defined integration of the various aspects of the aircraft 
takeoff and landing problem is presented in the form of a generalized 
computer program. Total aircraft system performance ib evaluated 
during the glide slope, flare, landing roll, and takeoff. 

The flight dynamics of a generalized, rigid body, aerospace 
vehicle are formulated in six degrees of freedom; a flat, nonrotating 
Earth is assumed. The independent equations of motion of up to five 
oleo-type landing gears are also formulated. 

A control management formulation is developed to automatically 
adjust control variables to correct errors in the vehicle's dynamic 
state. Stability in the small is used to maintain stability ' i n the 
large. 

The equations of motion are integrated using a generalized variable- 
step Runge-Kutta technique. 

The formulation is programmed for the CDC 6000 and Cyber 70 Computer 
Systems. The program is programmed in Fortran Extended using the Scope 
3.4 operating system. 
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V 
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Total applied vector moment 

Vector moment of the ground reaction about the k*'*' strut axle 
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Constant braking moment array, the use of which is determined 
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Braking moment applied along k th strut axle - ft-lbs 
Braking moment lower limit array - ft-lbs 
Braking moment upper limit array - ft.-lba 
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Total vector moment of all gear ground reactions about the 
aircraft mass center 

Vector moment of the strut ground reaction about 

the aircraft mass center 
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Body axes components of the vector - ft-lbs 


Total vector moment about the nominal mass center of all 
the forces applied to the vehicle 

Total mass - slugs 

Total mass of k tJl body - slugs 

Total vector moment of applied force* to body about 
oth body refarence point 

t.h 

Secondary piston mass in k L strut - slug 
Total mass of K + 1 bodies - slug 
Throttle setting 

Actual throttle setting rate array - per second 
Actual throttle setting array 

y.etual £.i.ii>>.cle settings lor engines A and B in Common Two 
Engine Logic 

Engine reverse throttle constraint array 
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dC 
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Desired throttle setting array 


Desired throttle settings for engines A end B In Coosaon Two 
Engine Logic 

Desired throttle setting in Cotnon Engine Logic 
Fixed throttle setting array 

Constant engine throttle setting rate magnitude used in 
control response - per second 

Component of ground reaction normal to axis of '■*’ strut-pounds 
Landing reveres throttle setter array 


Landing reverse throttle settings for engines A and B in 
Coonou Two Engine Logic 

Landing reverse throttle setting In Common Engine Logic 
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Nj(TN) Engine yaw moment array ~ ft-lbe 
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Two Engine Logic 
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P ok 
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PS 


m 


PS. 


PS. 

p 
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No landing roll Indicator - used t< stop 
program at impact 

Body axes system aerodynamic normal force, in the (-1 ) 
direction-lbs 

th 

Number of tires on k strut axle 
2 

Pressure - lbs/ft 
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tH 

Air compression force acting on the k tn strut - lbs 
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Upper air chamber pressure - lbs/ft 
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Lower air chamber pressure - Ibs/ft 

Landing tail-down constraint angle - 
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2 

Nominal tire foot print pressure - lba/ft 
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2 

Preload pressure of upper air chamber - lbs/ft 

2 

Preload pressure of lower air chamber - lbs/ft 

Phugold control sensitivity - deg/ft 

Roll (aileron) overcontrol constant - deg/deg 

Angle of attack overcontrol constant in glide slope - deg/deg 

Angle of attack overcontrol conatant for LAP > 2 -deg/dag 

Sideslip (rudder) overcontrol constant - deg/deg 

Euler jaw angle ovcrcor.tval constant - deg/deg 

Body axes components of the vector m - rad/sec 

2 

Required dynamic pressure in glide elope - lbe/ft 
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Dynamic pressure - Ibs/ft^ 

Pitch rate in flare based on a, and flare acceleration - rad/sec 
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Vector displacement of aircraft mass center from origin 
of the inertial axes systo" 

2 

Inertial acceleration of the platform origin - ft/sec 


Inertial axes components of the vector sum (R^) + r^ - ft 


CH 


Vector displacement of parachute attachment point measured 
from aircraft mass center 


Body axes components of the vector sum (R^) + r^ - ft/sec 


_ 


gR 


Inertial axes components ot the vector sun 
R - R r + (R^ + r^, the velocity of the 


th 


k strut axle as seen by the runway axes origin - ft/sec 

Angle of attack rate feedback constant in glide slope - sec 

Angle of attack rate feedback constant for IAP >_ 2 - sec 

Sideslip angle rate feedback constant - sec 

Euler roll angle rate feedback constant - sec 

Euler yaw angle rate feedback constant - sec 

Range from the ecarting point - nautical miles 

Displacement vector of runway origin from origin 
of inertial axes system 

Magnitude o£ the vector R R ; the component of the vector 

K , in the 1 direction 8 (there are no other components) - ft 
gR xg 
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RF 


f 


Total vector displacement of the k* strut 
mess center from the origin of the inertial 
axes system 

Inertial vector acceleration of the secondary piston 

Vector position ^ the body reference 
point from the o body reference point 


Body axes components of the vector (R^ “ ft 


Runway length - ft 

Vector displacement from the aircraft mass 

center to the point of application (l.e., the 

tire footprint) of the k 1 * 1 strut ground reaction force 

Vector displacement of the aircraft mass 
center from the origin of the runway axes a/stem 

Vector displacement of the k th strut axle 
from the origin of the runway axes system 


k c ^ strut axes components of the vector sum 
) + r kc> the vector displacement of the 
strut mass center from the aircraft mass cantar 


Total vector velocity of the bottom surface 
of the k c ^ r. "rut tires as seen by the 
origin of the runway axes system 

Vertical plane glide slope position rate feedback 
constant - sec 

Horizontal plane glide slope position rate feedback 
constant - aec 

Vector displacement of the point P(x,y,x) from the 
origin of the body - fixed axes system 
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Vector position of mass center of the K + i 
bodies from the eth body reference point 


Fully extended position of strut axle 
from origin of k fc h strut axle system - ft 

Instantaneous vector displacement of k** 1 strut axle 
from origin cL k th strut axes system 


Vector displaecraen 
from origin of k*" 


t of k C ^ strut mass center 
strut axes system 


Veci r displacement of k C strut secondary 
platan upper surface from origin of 
k^ otrut axes system 


Undeflected tire outer radius - ft 

2 

Aircraft reference area - ft 

Sum of the forcea resisting the k 1 * 1 strut 

movement(that portion of the kth strut 

ground reaction transmitted to the main airframe) - lbs 

2 

Parachute reference area - ft 


A scalar variable 


Constant acceleration of j 
Final value of a 
Final value of a 
Initial value of e 
Initial value of a 

Displacement of k** strut from fully 
extended axle position - ft 

Maximum allowable displacement of k C ^ strut - ft 

Distance betveen k th strut axle and k C ^ 
strut mass center - ft 

Displacement of k C ^ strut secondary piston 
from Its extension stop - ft 

Distance betveen k 1 ^ strut axle and k 6 * 1 
strut secondary piston extension stop - ft 
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®k2STOP k** 1 secondary piston acceleration as seen by the k** 1 

axes system which would exist if the secondary piston were 
against either of its stops (i.e., the value of Equation 
(58 of Appendix II with f , _ removed, used to obtain the 
stop contact force, 

s k2T Maximum displacement of k C ^ strut secondary piston (distance 

between secondary piston compression and extension stops less 
piston height) - ft 

®Rk X zk component of inertial acceleration of k** 1 strut axes 

system - ft/sec 2 

f Thrust vector in flare 


T 

T(IN) 


T d (!N) 


dA 


Mb 


*0 


Magnitude of the vector T - lbs 

Actual thrust array - lbs 

Desired thrust in Common Engine Logic - lbs 

Desired thrust vector in glide slope 

Desired thrust - lbs 

Desired thrust array - lbs 

Desired thrust for engines A and 3 in Common Two Engine 
Logic - lbs 

Total thrust required in Common Two Engine Logic - lbs 
Engine thrust obtained from TTFS subprogram - lbs 
Impact time - sec 

Lower thrust limit allowed in flare - lbs 
Time after impact to stop landing roll - sec 
Upper thrust limit allowed in flare - lbs 



Body axes components of net thrust vector - lbs 
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TF(I),TF(IN) Fixed throttle Indicator array 


t 


Time - aec 


Sk 

Ski 

Sk2 


C ch 

S 

t 

r 

t 

rv 


rl 


C r2 


Time after Impact to start braking - aec 

First time after impact to change array - sec 

Second time after impact to change Ig array - aec 

Time after impact to daploy parachute - aec 

Time required to perform the hg state change in flare - sec 

Time after impact - aec 

Time after Impact to reverse engines - aec 

First time after Impact for sequencing engine 
conditions in landing roll - aec 

Second time after impact for sequencing engine 
conditions in landing roll - aec 


t Elapsed time to perform constant acceleration 

change fro*^initial state (s,> s* ) to final 
state (Sf, ) of a scalar variable, a - sec 

t Time after Impact to set spoiler aerodynamics 

* p staging Indicator - sec 


t #fc Time after Impact at which nose over begins - sec 

t^ Time required ta perform the state change in. fiar« _ Me 


u 



Body x>axis component of inertial 
velocity of aircraft mass center - ft/sec 

t 2 

Body x-axla component of the vector V - ft/sec 


Body axes components of the vector V - ft/sec 


Body axes components of the vector V - ft/aed* 
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V 



V. __ 

hTD 

V 

■ 

V 

stall 


Body axes components of the vector V w - ft/sec 


Inertial velocity vector of mesa center (in body fixed 
axes system) 

Velocity vector rate as seen from the body fixed axes 
system 

3 

Volume - ft 

Velocity of piston pushing fluid through an orifice 
Airspeed vector 

Magnitude of the vector V - ft/sec 

ft 

Desired airspeed vector 

Component parallel to runway of k th strut axle Inertial 
velocity - ft/sec 

Airspeed for takeoff - ft/sec 

Desired inertial vector velocity down the glide slope 
Magnitude of the vector - ft/sec 

Inertial velocity magnitude error in glide slope - ft/sec 

Inertial velocity vector 

Magnitude of the vector - ft/sec 

Component parallel to runway of velocity of k** 1 strut tire 
footprint as seen by the runway coordinate system 

Initial valus of In flare - ft/sec 

Speed parallel to runway to stop landing roll - ft/sec 

Aircraft stall speed - ft/sec 
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P 

Z 

P 


LIST OF SYMBOLS (Contd) 

Tire footprint velocity of i C ^ wheel - ft/eec 
Runway exes component of the vector - ft/sec 

Runway axes components of the vector R^, - ft/sec 

Wind vector 

Initial value of In flare - ft/sec 

Preload volume of upper air chamber - ft"* 

3 

Preload volume of lower air chamber - ft 
Variable-step Runge-Kutta(lntegratlon technique) 
Weight vector 

Width of a tlra ou the k*** strut - ft 

Inertial axes components of the vector R - ft 


Inertial axes components of the vector - ft/sec 


Inertial axes components of the starting point of the 
vehicle center of mass - ft 


Components of the vector^ in a platform coordinate system (in this 
formulation Identical to the Inertial frame - ft) 
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Distance down runway to stop landing roll - ft 


Body axes components of the vector r - ft 


Center-of-gravity position from the body axes 
origin in the l x direction - ft 

Reference center-of-gravity position from the 
body axes origin in the 1^ direction - ft 


Body axes components of the vector R R - ft 


Aircraft acalar distance down runway - ft 



Runway axes components of the vector (R-R gK ) ~ ft 



*111 

*RF 

*RI1 


*RF2 

\k) 

T tk 

*ok 


Runway axes components of vector (R-R ) - ft/sec 

gx 


Distance down runway for touchdown - ft 

Landin’ speed parallel to runway for touchdown - ft/sec 

First distance down runway for sequencing 
engiue conditions in takeoff roll - ft 

Second distance down runwey for sequencing 
engine conditions in takeoff roll - ft 


Runway axes components of the vector - £t 


£ 
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*TD 

T 


y N CiK) 

V 1 *) 


Initial value of in flate; distance down runway before 
which the aircraft must not touch down - ft 

Body axea system aerodynamic side force, in the 1 
direction - lbs * 

Body y-axla component engine position array - ft 

Body z-axis component engine position array - ft 
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Angie of attack - deg 

Angle of attack rate - rad/sec 
Desired angle of attack - deg 

Desired angle of attic*, rate - deg/sec 

Preceding value of desired angle of attack - deg 

Maximum allowed value of - deg/sec 

Upper limit on angle of attack - deg 

Lover limit on angle of attack - deg 

Angle of attack position error - deg 

Total angle of attack error - deg 

Angle of attack rate error ~ deg/sec 
Angle of attack for takeoff - deg 

Sideslip angle - deg 
Sideslip angle rate - rad/sec 

Desired sideslip angle - deg 

Sideslip angle position error - deg 
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Elevation angle of relative to earth - deg 
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Elevation angle of relative to runway - deg 

Desired flight path angle of V ^ -deg 
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AC 

AV 


Ax 


Ax 


C.G. 


Ao 

hi 

a 

AA 


Generalized force inpuC (body axes system) - lbs 


Gyroscopic momenta due to rotation rates p,q,r and angular 
momentum of rotating machinery (body axes system) - ft-lbs 


Generalized moment input (body axes system) - ft-lbs 


Jet damping moment vector 

Current integration Interval determined by VSRX - sac 
Allowed glide slope inertial velocity magnitude error - ft/sec 


Distance from aerodynamic reference station to aircraft mass 
center - ft 

Canter-of-gravlty position from the reference center of gravity 
in the 1^ direction; (* c G - G ) 

*ref 

Allowed magnitude of total angle of attack error, 0 ^-deg 
Allowed magnitude of total sideslip angle error, 8 #T - deg 


Allowed magnitude of total Euler roll angle error, - deg 
Allowed magnitude of total Euler yaw angle error, ♦ ^ - deg 
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Am 


5 hs 


Allowed percentage of magnitude for 

wheel angular speed errort - %/(100Z) (i.e., nondimensional) 

Roll control surface (aileron) deflection rate 
magnitude used in control response - deg/sec 

Pitch control surface (elevator) deflection rate 

from 6 _ to d _ in landing roll - deg/aec 
qi. qr 

Pitch control surface (elevator) deflection rate 
magnitude used in control response - deg/aec 


Altitude normal to runway above 
to begin ’’hold" mode - ft 






+d. 


6 

P 


Sum of the quantities r^, and e 
Tire deflection of a tire on the strut - ft 
Actual roll control surface (aileron) deflection - deg 



5 

P 


Actual roll control surface (aileron) deflection rate 
deg/aec 


6p d 


Desired roll control surface (aileron) deflection - deg 


Lower llmlc on roll control surface (aileron) deflection - deg 


6p * 


Nominal roll control surface (aileron) deflection - deg 


5 P 0 Upper limit on loll control surface (aileron) deflection - deg 


Aq Actual pitch control surface (elevator) deflection - dig 

6q Actual pitch control surface (elevator) deflection rate - deg/sec 


«q c 

6 «c2 


Pitch control surface (elevator) deflection for "bang-bang" 
control in glide slope - degrees 

Pitch control turface (elevator) deflection for "be.ig-bang" control 
for IAP 2 2 - deg 


Aq y 


Pinal value of pitch control surface (elevator) 
deflection in lauding roll - deg 
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T.ultlal (Impact) value of pitch control surface 
(tJevs r) deflection in landing roll - deg 

1 

■ 



Lover limit on pitch control surface (elevator) 
deflection deg 

a 

i • 


5q N 

Nominal pitch control surface (elevator) deflection - 

deg 

■ 


dq T0 

Takeoff pitch control surface (elevator) deflection - 

deg t 



6q u 

Upper lindt on pitch control surface (elevator) deflection - 

deg 1 



5 rd 

"aw control surface (rudder) deflection rate 
magnitude used in control response - deg/sec 

i 

| 

* 


fir 

Actual yaw control surface (rudder) deflection - nag 

r 

ji 

► 


4r 

Actual yaw control surface (rudder) deflection rate - 

deg/sec I 



6r d 

Desired yaw control surface (rudder) deflection - deg 




5r L 

Lover limit on yav control surface (rudder) 
deflection - deg 


k 


5C N 

Nominal yav control surface (rudder) deflection - dag 


i 

i ^ 


ar o 

Upper limit on yav control surface (rudder) 
deflection - deg 




Ac 

Fixed angular perturbation of glide slope In 
vertical plane - rad 

r 
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4o 

Fixed angular perturbation of glide slope In 
horixontal plane - rad 





Arbitrary runway profile - ft 




*GS 

Nominal glide slope angle - deg 




*1 

Error multiplier ror 





Incremental error In 
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Error multiplier for 





Incremental error In C. 
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Error multiplier for 
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Incremental error in 
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Error multiplier for 


c 8 

Incremental error in 


e 9 

Error multiplier for C # 


C 10 

Incremental error in C 

a 


C 11 

Error multiplier for C 

A 


e 12 

Incremental error in C 

A 


e 18 

Incremental error in C.G, 

. location - ft 

f 19 

Incremental error in Z 

at 

- alug-ft 2 

c 20 

Incremental error in I 

rr 

- slugrft 2 

C 21 

Incremental error in I 

» 

- alugrft 2 

C 22 

Incremental err :>r in . 

*7 

- alug-ft 2 

C 23 

Incremental error in I 

XI 

- alug-ft 2 

*24 

Incremental error in I 

y* 

- alug-ft 2 


C Intermediate i-uli direction 
n Intermediate y-axla direction 
0 ' ) 




♦ ’ 

8 " 

♦ " 
T" 


) 


Pitch-yev-roll sequence of cotetlon engine 


Pitch-roll-yaw sequence of rotation engine 
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Pitch angle of k^ strut axes relative to body axes - deg 
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Pitch-roll-yaw sequence of rotation angles for body axes relative 
to Inertial platform (not programmed) - deg 

Pitch angle of aircraft from horizon - deg 

Pitch angle of rotating machinery axis relative to body x-axls-deg 

Pitch rate of rotating machinery axis relative to body x-axis-rad/sec 

Coefficient of friction between a tire on the k th strut and 
the runway 

Coefficient of friction at wing gear support for k th strut 
Intermediate x«axla direction 

Radius vector from origin of Inertial axes system to the point 

*(*• 7 .*) 

Atmospheric density - slug/ft^ 

Magnitude of desired roll angle in glide slope - deg 
Desired Euler roll angle in glide slope A flare - deg 
Euler roll angle position error - deg 
Total Euler roll angle error - deg 

Actual Euler roll angle (based on pltch-yaw-roll direction cosine 
sequence) - deg 

Actual Euler roll angle rate (baaed on pltch-yaw-roll direction 
coaine sequence) - rad/sec 
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Horizontal flight-path angle; heading angle - deg 
Euler yaw angle position error - deg 
Total Euler yaw angle error - deg 

Actual Euler yaw angle (based on pitch-yaw-roll direction cosine 
sequence) - deg 


Yaw-pitch-roll sequence of rotation angles 


Yaw-pitch-roll sequence of rotation angles for body axes relative 
to inertial platform (not programed as such - used as names 

for T' e* +') - deg 
P P P 

Actual Euler yaw angle rate (based cn pitch-yaw-roll direction 
cosine sequence) rad/aec 

Inertial rotattm rate vector of body axes 

Rate of ch of Z vector; inertial angular acceleration 
vector of \ axes 

Absolute magnitude of constant control angular acceleration for 
any wheel - rad/aec* 

Inertial rotation rate vectcr of k th body axes 

Rotation rate of machinery within the body; about the machinery 
spin axis - revolutions par minute (R.P.M.) 

Wheel angular speed error array - rad/sec 

Rotational vector velocity of the tires on the k*^ strut 

Magnitude of the vector T ^ component of the vector 

(there are no other components) - rad/sec 

ch 2 

Wheel angular acceleration on k strut axle - rad/sec 

Desired (required) wheel speed erray to obtain desired "percent 
skid" - rad/aec 
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Desired (required) wheel acceleration array to maintain 
desired "percent skid" - rad/sec^ 



Body axes components of the vector w - rad/sec 


Inertial rotation rate vector of o^ body; inertial rotation 
rate vector of body axes system 



Unit vector in the direction of the vector 
for the k c ^ strut 


xlll 





AFFDL-TR-71-155 
Part II 


Notational Conventions 


( * ) 
(•* ) 
ill 

dt 

( ") 

<*) 

[ 1 

£ 

k 

(“ ) 
(“) 


"as seen by the k** 1 coordination system" 

First derivative of ( ) with respect to time 

Second derivative of ( ) with respect to time 

Total derivative of ( ) with respect to time 

The vector ( ) 

The tensor ( ); the matrix ( ) 

A matrix 

K 

Summation over all the gears of the vehicle; I 

k-1 

Total derivative of the vector ( ) 

( ) vector rate as seen from a rotating ares syst 
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Axes Systea 



Body-fixed exes system, origin at vehicle 
nominal mass center (Also, coordinate 
system of o c ^ body) 


Earth-fixed axes system, origin at sea 
level (Inertial frame), 1 normal to 
flat-Earth, positive down*® 





1 


xR 




Body-fixed strut axes system, directed 
downward along k fc k strut, origin 
along this line of action 


Runway axes system, origin at sea 
level, axes system fixed relative 
to inertial frame in present formulation 


Direction Coalnea 


12 3 

*1*2*3 

n l*2*3 


Matrix of direction cosines. Used to transfer 
quantities from Earth-fixed (Inertial Axes 
Systea to Body-fixed Axes System 


*kia • *k*J 

• 1 • 

*k31 9 Vjj 

*ai # l ci3 

• 1C 

_*C31 # *C33 


Matrix of direction cosines. Used to transfer 
quantities from Body-fixed Axes Systea to k** 1 
Strut Axes System 


Matrix of direction cosines. Used to transfer 
quantities from Earth-fixed Axes System to 
Runway Axes System 
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*U1 *L12 *L13 
*L2l *L22 *L33 

_*L31 *L32 *133 


Matrix of direction cosines. Used to transfer 
quantities from Runway Axes Systea to Body-fixed 
Axes Systea 


*Illk *U2k *I13k 


*X21k *I22k 
*I31k *!32k 


*I23k 

*I33k 


Matrix of direction cosines. Used to 
quantities from Runway Axes System to 
Axes Systea 


transfer 
k 1 ** Strut 
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SECTION I 
INTRODUCTION 

In the design cf an aircraft, the engineer is confronted with the 
problem of takeoff and lancing and the design of aircraft systems and 
techniques to perform this function. The final evaluation of these 
systems lies in the answer to the question: How does the aircraft and 
its system perform as a unit? The Takeoff and Landing Analysis (TOLA) 
Computer Program is the result of an attempt to generalize the aircraft, 
the main aircraft control systems, and the landing-takeoff situation 
into a single comprehensive calculation to answer this question. 

Various analyses and simulations have been developed which are 
rigorous and thorough on a particular aspect of the landing-takeoff 
situation (e.g., References 20 through 24). Other are cited in 
Reference 25, which covers virtually every aspect of the problem. 

The TOLA simulation answers the above question in the form of a 
well-defined integration of the various aspects of takeoff and landing. 
In the equations of motion the assumption is made that the main aircraft 
framr is rigid; however, the dynamic efforts of up to five independent 
landing gears are included in the equations. The position and velocity 
of each strut and secondary piston are obtained by numerical integration 
subject to position constraints (for example, the main strut must 
move within the limits of the fully extended position and strut 
bottoming position). The same form of solution applies to the 
aircraft itself. 

The purpose of this report is to unify and summarize the complete 
formulation of the TOLA Computer Program. TOLA is a FORTRAN IV 
modification to Option 2 (SDF-2) of the Six-Degree-of-Freedom Flight 

Path Study Generalized Computer Program of References 1 and 2. 
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The SDF-2 modification resulting in the TOLA Computer Program 
was undertaken in two parts which are documented in FDMG TM 68-5, 
"Derivation of the Equations of Motion for the Landing Gear and 
Ground Reaction Modification to SDF-2," and in FDMG TM 68-11, 
"Autopilot Equations and Logic for the Takeoff and Landing Analysis 
Modification of SDF-2." 
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SECTION II 

TAKEOFF AND LANDING ANALYSIS COMPUTER PROGRAM 


This section presents a general description of the Takeoff and 
Landing Analysis (TOLA) Computer Program formulation. The following 
will be discussed in order: the original six degree of freedom 
equations, the landing gear and ground reaction equations, and the 
autopilot equations. 

1. ORIGINAL SDF-2 FORMULATION 

The Takeoff and Landing Analysis (TOLA) Computer Program is based 
on two modifications to SDF-2, which is the second calculation option 
of References l and 2. SDF-2 has six degrees of freedom, but assumes 
a flat, nonrotating earth. Those portions of the original SDF-2 
formulation retained by TOLA are presented in Appendix I. Figure 1 
contains a summary of the steps performed by SDF-2 in TOLA, and most of 
the equations used to perform these steps are found in Appendix I. 

The exceptions are the equations for the landing gear forces and 
moments which are discussed in Appendix II. 


2. LANDING GEAR AND GROUND REACTION FORMULATION 


When each landing gear of a moving aircraft comes in contact with 
the runway, it is subjected to a force, which is the ground 
reaction between the tire and the runway. The point of application 
of this force, namely the tire footprint, is located by a vector 


*Pk from the aircraft center of gravity (eg). Therefore, F 
generates a moment about the eg , 
the k struts from 1 to K, 


TRk 


ftj-Rk = ^Pk x ^TRk ‘ Sun,nin 9 over 


■? 


TR ' L r TRV = F TRA T *o +F TRB , yo r ''TRC l zo 


l.„, + Fi 


rT, 


(I) 


M TR = IMtrx - m T*^* 0 ♦ M Tyi’yo + M Tz 1 zo 


(2) 


in the body-fixed axes system. 
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Figure 1. SDF-2 Summary 
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If the strut were a rigid portion of the total airframe, then 


TR 


and 


M_ D would be the total ground reaction force and moment 
I K 


transmitted to the airframe. However, the strut is able to telescope, 
with the result that part of each Fy^ is used to accelerate the strut 
relative to the airframe; the remainder of F 
(or "felt" by) the airframe. 


TRk 


is transmitted to 


Likewise, since ^TRk = R Pk x F TRk * then part of each My^ 
will be evident as a moment of the inertia force portion of Fy^; 
the remainder of Wy^ is transmitted to the airframe. 


Let F , F, .and F „ be the body-axes components of that portion 
xm ym’ zm J r 

of the total ground reaction force transmitted to the main airframe, 

and L , M . and N m be the body-axes components of ^hat portion of the 
m hi m 

total ground reaction moment transmitted to the main airframe. 


Then, 


and 


F xm 


f tra 


£m k S k a k3) 

k 

F ym 

r 

f TRB 

+ 

0 

F zm 


f trc 


Z™ k S k d k33 
k 




Mr* 


Z m k S k °kll R ky 


= 

My y 

+ 

~Z m k S k RRkCGx 

k 

Nm 




Z m k Sk a k l 3 R ky 
k 


( 3 ) 


( 4 ) 


where m^ is the strut mass, is the strut acceleration relative to 


the main airframe, a k ^ j are direction cosines relating the strut 
orientation to the body axes system, and and are 

moment arms of the strut inertia forces. 


£ 
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Equations 3 and 4 are included in the summation of forces and 
moments as follows: 


V 


T x -fl + Mg x + AF X 


F xm 

F y 

= 

Ty + y + + AF y 

■f 

F ym 

Fz 


T z + Np + Mg z + AF Z 


F zm 


L ^ 


l t +al t +£ 



M 


My 4- AMy +m 

+ 


N 


Ny + ANy +n 




The derivation of the strut acceleration, S^, requires the 
definition and analysis of the forces internal to the strut. A detailed 
derivation of the various terms in Equations 3 and 4 will be found in 
Appendix II of this report. The flow charts in Appendix II are 
summarized in Figure 2, where the circled letters correspond to those 
found in the flow charts. The basic equations of motion used in 
Appendix II are derived iri Reference 6. 

3. CONTROL MANAGEMENT FORMULATION 

The purpose of the TOLA control management formulation is to 
determine appropriate values for Sq, Sr, Sp, N(IN), and M„(I), which 
are the control variables for pitch, yaw, roll, throttle, and 
braking, respectively. This determination is made in three major 
steps called the maneuver logic, the autopilots, and the control 
variable response. The respective tasks of these steps is to determine 
what needs to be done, how to do it, and how much can be accomplished 
within a given time step. 

a. Maneuver Logic 

The maneuver logic is that portion of the control management 
formulation that is specifically oriented toward the takeoff and 
landing problem. It determines what needs to be done in two steps. 
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Figure 2. Landing Gear Sumrwry 
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The first step is the determination of which phase of the analysis 
is to be examined at the present time step. This is done in the 
problem phase logic in which a sequence of tests is performed on the 
present position of the aircraft relative to the runway. In this 
manner, it is determined whether the aircraft is in the glide slope, 
flare, landing roll, or takeoff roll phase. Each phase has its own 
logic for determining the desired state of the aircraft. 

The second step performed in the maneuver logic is to determine if 
an error exi ts in the present position, velocity, or acceleration of 
the aircraft relative to the runway. If so, the required dynamics are 
solved for toe desired angle of attack, a d , the desired roll angle, 

, and the desired thrust, . 

The values of ^, and T^ obtained in the above two steps 
determine what needs to be done in terms of aircraft orientation and 
thrust to correct kinematic state errors which may exist, depending on 
the prob.em phase. This accomplishes the main purpose of the maneuver 
logic, however, prior to entry into the autopilots, checks are made 
for possible input system failures. If the aircraft is in the landing 
roll phase, the condition of the brakes is determined at the present 
time step for later use in the Brake Autopilot. The final step for all 
phases in the maneuver logic is to determine the condition of all engine 
in the engine failure logic. This is done at the present time step 
for later use ,- n the Throttle Autopilot. 

b. Autopilots. 

The purpose of the autopilots is to determine the means by 
which th? desired vehicle orientation, thrust, and braking are to be 
accompli shed. The results are obtained in the form of Sq^, fcr^, 6p^, 
(IN) and Mg(I), which are, respectively, the desired pitch, yaw, and 
roll control surface deflections, the desired throttle setting array, 
and the actual braking moment array. These variables are determined 
in the Pitch, Yaw, Roll, Throttle, and B^ake Autopilots, respectively, 
which are discussed in detail in Appendix III. 
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c. Control Variable Response 

The purpose of the control variable response is to get a first 
order approximatioi. of the effects of control surface and engine lags 
on aircraft performance. Constant time rates of change are used to 
determine the measure of 5q d , 5r d> Sp d> and N d (IN) attainable within a 
given time step. The results are the actual control variables Sq , 
6r, 8p, and N(IN). 


d. Control Management Formulation Summary 

A detailed development of the logic and equations for the 
control management formulation is presented in Appendix III. The 
basic steps in the control management development are summarized in 
Figure 3. 



Figure 3. Control Management Formulation Summary 
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SECTION III 

SUMMARY AND DISCUSSION 

The Takeoff and Landing Analysis (TOLA) Computer Program is 
summarized in Figures 4 and 5. Figure 4 shows the various capabilities 
of TOLA as depicted in a perspective view of an inclined runway. TOLA 
can simulate aircraft control and performance during qlide slope, flare, 
landing roll, and takeoff subject to changing winds, engine failures, 
aerodynamic ground effect, runway limitations, and control variable 
limitations and lags. The (eg) oscillation seen in the landing roll 
represents the multiple landing gear dynamics. Runway perturbations 
may be input as a function of distance down the runway, as exemplified 
by the "bump." The capabilities to cut power; to actuate spoilers and/or 
ur>g chute; to have brakes locked, off, or on with constant or controlled 
braking, and to reverse engines are shown as functions of time. This 
demonstrates how the TOLA simulation allows the order of these effects 
to be changed. The takeoff is accomplished as an acceleration run to 
a predetermined airspeed, at which a desired angle of attack is s<_nt to 
the pitch autopilot. The above versatility is maintained in as generalized 
a manner as possible. 

Figure 5 is a block diagram summarizing the flight mechanics 
performed in TOLA. The instantaneous dynamic state of the vehicle is 
determined by the six-degree-of-freedom routine, which determines the 
accelerations caused by the aerodynamic, thrust, and landing gear forces 
and moments. These accelerations are integrated twice by the executive 
routine, which then sends the kinematic state information in terms of 
position and velocity to the control management routine. In the 
Maneuver Logic, errors in kinematic state determine the desired angle 
of attack, thrust, and roll angle. These are processed as forcing 
functions by the Autopilots, which determine the desired control settings. 
The Control Variable Response then determines the actual control settings 
attainable during a given time step. 
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TOLA development grew out of a need for comprehensive, quantitative 
analysis of aircraft takeoff and landing performance. The simulation 
attempts to generalize conventional powered aircraft, the main control 
systems, and the landinq takeoff situation into a single comprehensive 
calculation. The simulation does not perform the design function; it 
simply takes input data and evaluates performance. 

The concept a single, comprehensive, quantitative simulation 
of total system performance has yielded very promising results (see 
paper, "Capabilities of the TOLA Simulation," presented July 1969 at 
the AIAA Aircraft Design and Operations Meeting, L. A., Calif.) for 
the takeoff and landing problem. To date, the capability of the 
program as a design tool to do tradeoff studies in major system 
component design has only just begun. Even with the cursory results 
received so far, many questions come to mind: 

(a) What effect does limited runway length, changing winds, and 
engine failure have on a go-around decision for a particular situation? 

(b) How does a change in.the control schedule for the landing roll 
affect maximum gear loads? 

(c) What limitations would have to be placed on the landing if on', 
strut failed to brake or failed to extend from the fuselage? 

(d) With multiple-engine aircraft and a reverse capability, is it 
safe to have some engines in reverse during landing in view of possible 
engine failure? 

These are just a few of the questions that are within the capability 
of study by the TOLA simulation. In its interest to develop better 
technology, the Air Force Flight Dynamics Laboratory will continue to 
Improve the TOLA simulation and use it as a tool to study the takeoff 
and landing problem. 
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APPENDIX I 

ORIGINAL SDF-2 EQUATIONS 


1. DERIVATION OF EQUATIONS OF MOTION 

This section presents the derivation of the equations of motion, 
of a body in "inertial" space, as required for use iri the Takeoff and 
Landing Analysis computer program. The equations of motion will form 
a portion of the computation loop which is unaffected by the libraries 
of interchangeable subprograms describing alternate control systems: 
airframe aerodynamics, atmospheres, and geophysical parameters, or the 
data-monitoring subprograms to be incorporated. The several coordinate 
transformations and velocity and angle resolutions, which complete 
this central portion of the problem, are described later in this 
appendix. 


Since the equations involving the moments of inertia, aerodynamic 
forces, and thrust forces are greatly simplified if expressed in body 
coordinates, this system of body reference will be used. The two 
basic equations which define the motion of a body are: 

F= ^-(frV) (7) 

M = jjj- (H) (8) 



Figure 6. Generalized Inertial and Body-Axes Coordinate Systems 



14 V-AfT QuhJL l2$0-r^k’ ' 
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Numerical analyses of these vector equations require their resolution 
into vector components and definition of the scalar coefficients. 

These manipulations are discussed in detail in many texts in mechanics 
(e.g., References 7 through 14. The essential steps of the derivation 
are reviewed here, however, for completeness. 

To determine the displacement accelerations, consider a point P 
displaced from the origin of coordinate system x-y-z such that the 
vector F designates the point. Figure 1 illustrates the system. 

Let the origin of the coordinate system x-y-z be displaced from the 
origin of a space-fixed coordinate system X^-Y^-Z^ b y an amount and 
direction given by R. Further, let the coordinate system x-y-z 
rotate in the X^-Y^-Z^ space such that the vector, tU, defines the 
rotation. 

Then _ 

r = xl x + yly + 2l r 

R = Xgr xg + Yg7y g + Zg7 zg (9) 

= CU* I X + OJyTy + UI Z I Z 

The coordinate system x-y-z will be recognized as the body axes and the 
coordinate system X g“ Y g'^ g are the nonmoving "inertial" or Newtonian 
axes. The total velocity of the point is given by 

"p = "R + ^ = *g<*g + Vyg+V*g + (*~ytUj +zu)y) l x 

OO) 

+ ( y + *uj z - zwx) ly + (i-xu>y ♦ yu % ) 1/ 

It is more convenient to express the velocity of the body-axes origin 
in body-velocity components than in velocity components coincident with 
the 'inertial" reference coordinates. The vector If can be written in 
any coordinate system, so 

"R 1 V - ^g , xg + Tg l yg + Zg I yg - * 0 7 x + Y 0 ly + ^o'z 


t 


17 



AFFDL-TR-71-155 
Part II 

and Equation 10 may be rewitten as 


« • _ _ _ 

? = V + r = x 0 T x + y 0 l y -+-z 0 l z + ( x p -y p OJ z + z p tu y ) l x 

+ (y p + XpW 2 -Zp(U x )T y +(z p -x p uJ y +y p w x )T z 


(I!) 


where the subscripts o and p have been added to distinguish between 
the velocity components of the origin and the relative movement of the 
point P with respect to the origin of the x-y-z coordinate system, 
respectively. Differentiating Equation 11 gives the relation for the 
total acceleration to be 


p 


= [Vo - y'o^Z - 1 - *0^y] T* + [Vo + *O^Z “*o"x]T y + [*o -X 0 Wy +y 0 ^x] 'z 

+ [x‘ p -2y p cv z +2z p a) y -x p {cu z + u y ) +y p ( cu x u» y -u> z ) + Zploiy + aj x cu 2 )|T x 

02 ) 

+ [y p +2x p uj z -2z p ou x -y p (uj x + co z )+z p (u> y a> z ~w x ) +x p (w z + w x ui y )]T y 
+ [z p -2x p o) y + 2y p ai x -z p (aj y + <u x ) + x p ( w z uj x -uj y ) + y p (oL' x + ai y uJ z )] T z 


This acceleration relation is completely general and applies to any 
point on the body. In developing the equations of motion, the point 
of interest, P (x,y,z), is the center of gravity. If the center of 
gravity is assumed to move, relative to the body, along the x-axis 
only, the following simplification can be made. 


Vp : z p - 0 



The components x , y , z and w , w , w. are more commonly known as 
ooo x y if 

u, v, w and p, q, r, respectively. The components u, v, and w are the 
velocities of the reference point on the body. Making the above 
substitutions gives 


F x - M | u - vr + wq +xp -xp(r 2 + q 2 ) j 

F y = M [ i + or - wp +2x p r + xp(r + pqjj (13) 

F z - M [w -uq + vd -2x p g + *p (r P -q)j 


IB 


# 






AFFDL-TR-71-155 
Part II 


Since most vehicles are designed to have small center-of-gravity 
travel, the acceleration and velocity of the center of gravity are 
both very small quantities and may be omitted from the problem 
formulation. If the reference point is further restricted to be the 
center of gravity, then x p and its derivatives may be omitted from the 
equations and the components u, v, and w are the velocities of the center 
of gravity. In matrix form the equations reduce to the following: 



Note that in the analysis of flight-test data, where the output of 
accelerometers mounted away from the center of gravity are used to 
record the motion of the body, the complete form of Equation 13 must be 
used. Note also that, although Equation 7 states Newton's Law as the 
time derivative of the momentum, a formal differentiation of 7//y, 
assuming/?/ to be a function of time, has not been performed in the 
derivation of Equation 14. Such a formal differentiation gives 

r .»jZ ♦ 

This differentiation leads to erroneous results, however, since the 
residual momentum of the expelled gases has not been accounted for by 
this procedure (see References 6 and 14). The equation should be 



when the residual momentum of the expelled mass is properly considered. 
Here, c Is the velocity of the expel’ed mass with respect to the 
continuing body. The contribution 7/7 2 Is the momentum-change portion 
of the thrust and is Included in the summation of external forces. 
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There are additional accelerations produced which are unique to 
configurations which have very large fuel flow rates and which have the 
thrust nozzle located a considerable distance from the center of gravity. 
These accelerations, linear and angular, are the so-called jet-damping 
contributions. The term is a correction to accelerations computed on 
the basis of only the externally applied forces (or moments) and 
accounts for the moment of momentum which is imparted to the fuel by 
the pitching velocity of the body. The derivation of this contribution 
is considered in greater detail in Reference 1. The principal 
contribution to the equations for linear acceleration are in the y- and 
z-directions and have been added to the expressions of Equation 14 to 
give the following result. 



-r q 
0 -p 
P 0 



V 

W 


0 

-Ztilrti 
+ 2fcq-t z 


(15) 


The relations expressing the rotational motion are obtained in a 
straightforward manner. The components considered in this anslysis 
come from three basic sources: the time rate of change of the moment 
of momentum, the gyroscopic moments which arise from the rotating 
machinery of the vehicle, and .he externally applied moments. The 
moment of momentum of a body (or angular momentum) about its center 
of gravity, in terms of its components, is given by 

"y 

"z 

or, sinco w , w ,and w are p, q, and r, respectively: 
x y z 





H, 


1 x x 1 x y ~'xz 

H y 

r. 

-I*y iyy ’ 1 y ? 

H z 


"’xz _, yz 'zz 

_ 


. 


H - 






+ [ • *yP + 1 yyP *yz'] 'y 

+ [ ~* xz P — l yiR +l zz , ]^z 


( 16 ) 
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★ 

The required differentiation of the moment of momentum gives 

F. = ['xxP + 'xxP + < ^2z -l yy5q ^ ~lyz(P ? ' _^^ ) _l xr( ,r+ P‘^ , 
-l xy (q-pr)-i X2 r -’ xy q]Tx 

+ [l yy <5 lyyO "*■ ( IxX ~ I ZZ ^ P r ’ ^ XZ I r “P ) ~ I Xy l P ‘FR 1 ’ ) 
-ly 2 l r_ pq) _ ixyP — lyz r ] ^ 

+ [*ZZ r + ^ZZ r + (lyy _l xx)p j ~ I xy(p 2-C l 2 >~ , yz*5 + P''/ 

-ixz(p-q r )-ixzP“*yz' 1 ]^z 


7 is general practice at this ooint in the derivation of the equations 
of >n ("0 assume that the reference . xes of the aircraft are principal 
axes and that the moments of ir ;rt a do not vary with time. This 
conveniently elimates the prod cts of nertia and the t me derivatives 
of the moments and products o ine tia, respectively. However, it is 
desired to have a more general applicability than this for the computer 
pro~ram being developed, and these terms will be retained. The 
inclusion of the time derivatives of the inertia implies that all moment 
of momentum has been remove' from the mass being lost by the boiy. This 
assumes that the gases have no swirl after they have left the body. 

Staging and dropping of di icrete masses from the body introduce 
discontinuities in the macs and inertia properties of the body. The 
solution must not proceed across these discontinuities. Therefore, the 
integration of the equati, is of motion will be interrupted when mass is 
dropped and automatically reestablished immediately thereafter (see 
Section 3, Stages and Staging). 


The time rate of change of inertia noted here refers tc that change 
occurring at constant mass only. 
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The .jet damping contribution to the expressions for angular 
acceleration (from Appendix I of Reference 1) is: 

AM 0 = T x -qm 1 ^ Ty - (18) 


The expression for the total angular acceleration due to the time rate 
of change of the moment of momentum, including jet damping, is 
conveniently given in matrix form: 



The torques due to precession and changes in rotational speed of 
rotating machinery abroad a vehicle which is free to gyrate in space 
can contribute significantly to the angular accelerations which the 
vehicle experiences. Appendix II of Reference 1 derives the torques 
generated by the precession of rotating machinery in general terms and 
simplifies these relations as required for the solution of the following 
problems: 

(a) The motion of an aircraft powered by an engine with a rotating 
mass which is fixed in its orientation with respect to the reference 
axis of the aircraft. 

(b) The motion of an aircraft powered by a rotatinq-mass engine 
which can be rotated in a plane parallel to the plane of symmetry 
(e.g., convertiplane which is in the transition from vertical flight 
to forward motion, or vice versa). 
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(c) The motion of a satellite in which motors are being operated 
(by the proper selection of reference ax^s). 

The gyroscopic moments due to the rotational rates p, q, and r and 
the angular momentum of the rotating machinery are approximated as 
follows: 

AL r ' -1 xr to r (q + $r) sin# r 

AM r = lxro>r (psin0 r +rcoi9r ) (20) 

AN r = -i xr (q +9 r )(jj r cos 9r) 

The complete rotational equations of motion are, therefore, from 
Equations 17, 18, and 19 

M = Ll X + Mly + Nl 2 

in which 

L = XxP'E'XXP^'i^ZZ - ) 

-l xz (f +pq)-l xy (q-pr)-i xz r - ( xy q 
-pmf-J - I xr co r (q + # r ) sin# r 

M - I y y Q ^lyyq^llxx - 1zzlp r_ *xz( r ~ P ) 

-l xy (p+qr) - l yz (r-oq)-i xy p-i yz r 
-qm I m + 'xr cu r ( p sin# r + r cos 9 r ) 

M " 1 zz r + ( zz r 1 ^yy _ *xx 1 PP — ^xy 1P ~P 1 

— lxz ( P •q r )- | yz'P + P r > -l xzP-iyzP 
- r ml* - (q + # r )u) r cos# r 
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These relations, 


L 

M 

N 


'xx “ 'xy 
>xy 1 yy 



written in matrix form, are: 
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-l xr w r (q +9 r ) sin# r 

1 

X 

*■< 

l 

*< 


q 

+ 

i xr u> r ( p sin@ r + r cos@ r ) 

'xz 'yz 'zz 
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-'xcCUrfq + 9 r ) cos 9 r 


Equations 15 and 21 constitute the general six-degree-of-freedom 
equations of motion which will be used in the computer program. The 
program instructions will provide for the removal of certain 
combinations of terms, as follows: 


(a) All product of inertia terms for the case where the body is 
inertial *y symmetrical about the x-axis. 

(b) The product of inertia terms I and I , which are zero when 

xy yz 

the x-z plane is a plane of symmetry. 

(c) The terms containing the time rates of change of inertia, 
products of inertia, and mass. 

(d) The gyroscopic contributions of rotating machinery. 

(e) The jet damping terms, both forces and moments. 

2. COORDINATE SYSTEMS AND COORDINATE TRANSFORMATIONS 

This section presents a description of the reference coordinate 
systems chosen for the Takeoff and Landing Analysis computer program. 

The coordinate transformations required to relate the various parameters 
of the computation to the several coordinate systems are also derived. 
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The coordinate transformations required in the program may be categorized 
as follows: 

(1) Transformations inherent in solving the basic equations of motion. 

(2) Transformations to provide input data to the guidance and 
autopilot. 

(3) Transformations to oresent readout data in the most desirable 
form and auxiliary transformations which may be required for the 
definition of certain special parameters. These transformations may be 
deleted from the program when they are not required, 

a. Coordinate Transformations for Basic Equations of Motion. 

The coordinate systems and transformations required to describe the rigid 
airframe motion in six degrees of freedom and the coordinate transfor¬ 
mations which relate the aerodynamic angles and velocities to ground- 
referenced velocities in the presence of winds are presented. 

(1) Body-Axes Coordinates 

The equations o f motion (Section 4) are solved in a body coordinate 
system (see Figure 6). The origin of this system is at the center of 
gravity of the aircraft, with the x-axis along the geometric longitudinal 
axis of the body. The positive direction of the x-axis is from the 
center of gravity to the front of the body. 

The y-axis is positive to the right, extending from the center of 
gravity in a water-line plane. The z-axis forms a right-hand orthogonal 
system. This coordinate system was chosen because inercial characteristics 
are thus made independent of attitude. 

Accelerations and velocities computed in the x-y-z body axis must be 
related to velocities arid accelerations referenced to a fixed point on 
the surface of the planet to (a) describe the motion which a fixed 
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observer would sense, and (b) compute the aerodynamic forces on the 
body immersed in an atmosphere with winds which are referenced to a 
point on the surface of the planet. 

(2) Inertial Coordinates. The inertial coordinates used in 
this analysis are the x ~“ Y g -7 g axes. X g and lie in the plane tangent 
to the earth's surface, and is the inward or downward normal vector 
tc this plane. The effects of the earth's curvature and rotation rate 
are truly negligible in the takeoff and landing phases of aircraft flight. 
Thus, the inertial coordinates defines a flat, nonrotating earth. 

(3) Direction Cosines. The direction cosines relating the 

body x-y-z axes to the inertial coordinate system X^-Y^-Z^ are obtained 

in the following manner. Let 1, 1 , 1 be unit vectors along the 

x y z — — — 

body axes, x, y, z, respectively, and let 1^, 1 Yg l Zg >. be unit 

vectors along the inertial axes, X g , Y g , Z g , respectively. The 
direction cosine matrix relating these two sets of unit vectors will be: 





$\ iz Is 


Ug 


*y 

n 

m | m 2 rn j 


>yg 


1 

•r' 

1_ 


n | n 2 rij 



Performing the matrix 

multiplication indicated 

gives 


'x T + Az'vg + ^3 l Zg 

ly = m ! (23) 

7z = n,Tx g +n 2 l Yg + n 3 T Zg 
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The derivatives of 1 ,1 , 1, with respect to time in terms of their 

x y z 

components in the inertial system are found by differentiating 

Equation 23. These derivatives are: 

• * • • 

'x ~i|IXg + / 2 ^Yg + <? 3 l Zg 

7 y = m|7xg+m 2 T Yg +m 3 Tzg 124) 

•z = n l~Xg +n 2 l Yg + ^a~Zg 


The derivatives of 1 , 1 , 1 

x y z 

on the change in direction of 


with respect to time are dependent only 
the unit vectors. Therefore, 


l x = tux i x = rly - qi z 

7 y = cuxly = pl z -rl x ( 25 ) 

T z = cux 7 z = qT x — pl y 


where 


w - pl x +ql y + rl 2 


Equating the relations for l x from Equations 24 and 25: 



+ i 


2 "y g + 4'i 


zg 


f lv * P'a 


Substituting the relationships for 1 and 1 , respectively, gives the 

y ^ 

relation: 

+ 4 ! Yg + ^3~Zg = r < m • *Xg + ^'Yg + ^s'Zg* 

<261 

-q(nj l Xg +n 2 1 Yg + n 3*Zg) 


By using the component properties of a vector, the relations 


f 1 = rm, - qn, 

(27) 

/ 2 rm 2 pn 2 

(28) 

i3 - rm 3 ~ q n 3 

(29) 


are obtained from Equation 26. 
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Performing the same 
the time derivatives of 


operation for the 1 and l z components defines 
the remaining direction cosines. These are: 


m, = pn, — r/. 

(30) 

nig — png r,4 

(31 1 

m 3 ' P n 3 -4 

(32) 

n, - q/, - pm. 

(33) 

n 2 = P'4 — P m 2 

(34) 

n 3 = p 4 — P m 3 

(35) 


Equations 27 through 35 are integrated to obtain the instantaneous 
values of the direction cosines. This method of calculating the 
direction cosines has been selected instead of the usual evaluation by 
means of the Euler angles because, regardless of the order of rotation 
selected, there are points at which certain Fuler angles become 
undefined. The direction cosines evaluated by this method are always 
defined (Reference 1). The method by which the orthogonality of the 
direction cosines is maintained is described in Appendix III of 
Reference 1. The Euler angles may be calculated from the direction 
cosines if desired; however, they are not required for component 
resolution. 

(4) Inertial Velocity 

The components of inertial velocity in the body coordinate 
system, u, v, and w, will be resolved into velocity components X^, Y , 
and 1 in the inertial coordinates. Since components of inertial 
velocity are known in body coordinates, a resoTution of components 
using the direction cosines given in Equations (27-35) will give 


<?7 
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components of inertial velocity in the inertial coordinate system, 
as follows: 


iz 

/, 


m * 


m *j 


m, 


(36) 


The planet-referenced velocity may be calculated from its components: 


Vg =/*l + *5 + Z\ 


(37) 


(a) Nine integrations are involved in the present method of 
computation, instead of the three that are normally required when the 
Euler-angle rates are integrated to give the Euler angles. However, 
a coordinate transformation is required to obtain the rates, and the 
sines and cosines of the angles must also be computed in the usual 
direction cosine computation. The machine time required for the two 
methods of computation is comparable. 


The flight-path angles are computed in the takeoff and landing problem: 



(38) 


r 



(39) 


The three inertial components of winds will be introduced in a 
tabular listing with altitude as the independent variable. Let the 
three components of wind be written as follows: 


V 


W 



(40) 
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The airspeed vector is given by: 

V Q = v g -7* (41) 


where V is the velocity relative to the atmosphere. The three 

a 

inertial components of airspeed are: 


V Q = 


iXg-XgJ 


1 Y g~ Xg* 1 


+ (z n -z. 


w ! ' z g 


(42) 


The elevation flight path angle of the airspeed vector is: 


Yo - s |n 


-(Z g -Zg r ) 


(43) 


Resolving inertial wind components to body-axis components requires the 
same direction-cosine matrix used in Equation 22, and the body components 
of winds are: 


u w 





1 

X- 

«o 

< 

_1 

v w 

• 

m, 

m # m. 


v 

"w 


n i 

n, n. 


Zgw 

_ 


_ 



— — 


(44) 


The body components of airspeed are determined by substracting the body 
components of wind from the body components of velocity. 

The body components of airspeed will be used to compute the angle of 
attack and sideslip. 


a ■ 

fi ■ 



(45) 

(46) 
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The definitions of angle of attack and sideslip are consistent with 
the aerodynamic data normally obtained from wind tunnel tests of sting- 
mounted models because of the manner in which the sting may be moved. 

If aerodynamic data as obtained from turntable-and-strut mounted models 
are used, an alternate definition may be required depending upon the 
procedure used in data reduction. 

b. Guidance and Autopilot Coordinate Transformations 

The vehicle attitude information taken from the gimbals of a 
stabilized platform and the outputs of platform-mounted accelerometers 
may be required in certain autopilot and guidance-system computations 
in the Takeoff and Landing Analysis computer program. This section 
presents the derivation of the equations relating accelerometer and 
attitude information to data computed in the central program. The 
method for deriving coordinate transformations for any gimbal arrange¬ 
ment is presented for reference. 

(1) Gimbal Arrangements and Rotation Sequences 

Three frequently used gimbal arrangements will be considered 
in this section. Each gimbal is equivalent to an intermediate 
coordinate system in a series of Fuler-angle rotations. Reading from 
the outer gimbal to the inner gimbal (and neglecting redundant gimbals) 
the arrangements considered are: 

(a) Yaw-Pitch-Roll 

(b) Pitch-Yaw-Roll (only this sequence is programmed in 

TOLA) 


(c) Pitch-Rol1-Yaw 

where the analogy between coordinate system rotations and gimbal 
movement is used. Other gimbal arrangements are possible but the 
three discussed in this section are the ones most frequently used. 
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Transformations for the alternate arrangements can be obtained with 
these same techniques. 

(2) Euler Angles. In the central program, the direction 
cosines relating the vehicle body-coordinate system to a fixed inertial 
system are calculated by integrating functions of the body angular 
l velocities, p, q, and r. The direction cosines relating the body and 

inertial systems are determined by the cosines of the angles between 
the various axes of the coordinate systems and are dependent only upon 
the position of the body coordinates referenced to inertial coordinates. 
That is, the order of rotation selected to arrive at a certain 
orientation does not alter the numerical valuer of the direction cosines 
for that orientation. 

Each individual direction cosine may, therefore, be defined in 
terms of the Euler angles from a given sequence of rotations. These 
definitions will provide the Euler angles of the body with respect to 
the platform coordinate system for the three rotational sequences 
selected. 

The direction cosines, in terms of three sets of Euler angles, 

% will be derived using the method of Reference 15. The technique used 

is to find the direction cosines for each individual rotation in a 
sequence and determine the complete transformation by multiplying the 
individual direction cosine matrices. The overall picture of the 
rotations is best observed on a unit sphere diagram. The points on 
the unit sphere represent the intersections of the coordinate axes 
with the surface of the sphere. 

The order of rotation and the axis about which rotation occurs can 
be described using the following diagram. 
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\ 


AXIS AND ROTATION ORDER 


X 

t 

x 

u 


I® 

n 


z 

$ ♦ 

z 

C 


1. 


2 . 


3- 


This diagram indicates that the first rotation is about the inertial 
Z-axis through the Euler angle 'V. The second rotation is about the 
intermediate axis n through the angle 0. The final rotation is about 
the body x-axis through the angle 


The derivation of each sequence of rotations will proceed in the 
following manner: (a) the order of rotation will be defined; 

(b) the unit sphere showing all three rotations will be presented; 

(c) the individual rotations will be shown in three separate diagrams 
that contain the plane perpendicular to the appropriate axis of 
rotation; (d) the direction cosines for each individual rotation will 
be written in this manner: 



C 5Y 

c*z 

°nx 

C hY 

c gz 

CxX 

C xY 

C xZ 


where C^. is the cosine of the angle between the i and j axes; and 
(e) the matrix of direction cosines relating the inertial and body 
coordinates will be determined by matrix multiplication. 
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(a) Yaw-Pitcli-Roll Rotation 



ORDER 


1 . 


2 . 


3- 


Figure 7. Unit Sphere for Yaw-i : ch-Roll Sequence of Rotation 
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THIRD ROTATION 

x 1 0 0 x 
y = 0 Cos Sin $ t) 
z 0 -sin Coo £ 

The transformation matrix is given by 

x X 
y = 0 o y * 
z z 



or, in terms of the planar rotation matrices, the intermediate axes are elim¬ 
inated by 


X 


1 

0 

0 


Cos 0 

0 

-Sin e 

Cos y 

Sin y 

0 

X 

y 

= 

0 

Cos 0 

Sin <f 


0 

1 

0 

-Sin y 

Cos y 

0 

Y 

z 


0 

-Sin $ 

Cos (f 


Sin 0 

0 

Cos 0 

0 

0 

1 

z 


The direction cosine elements of the transformation matrix are obtained by per¬ 
forming the indicated multiplication. For the yaw-pitch-roll rotational 
sequence 


| (Cos 0 Cos ♦) 

(-Cos $ Sin y 
+ Sin 0 Sin 0 Cos y) 


(Cos 0 Sin y) 

(Cos $ Coe y 
+ Sin ^ Sin 0 Sin V) 


(-Sin 8) X 

(Sin Cos 0) X 


(Sin / Sin y 
-* Coo $ Sin 0 Cos y) 


(-Sin 0 Cos y 
•f Coe Sin 0 Sin y) 


(Cos $ Cos ©) II Z 
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(b) Pitch-Yaw-Roll Rotation 



Figure 8. Unit Sphere for Pitch-Yaw-Roll Sequence of Rotation 


FIRST ROTATION 


l 


Cos 9’ 0 -Sin 9* 


X 

Y 

= 

0 10 


Y 

n 


Sin 0' 0 Cos ©» 


Z 



SECOND ROTATION 


X 


Cos v' Sin '| fl 0 


6 

c 


-Sin Cos 0 


Y 

n 


0 0 1 


n 
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THIRD ROTATION 


X 


1 0 

0 


X 

y 

= 

0 Cos 0' 

Sin 0* 



2. 


0 -Sin 0* 

Cos 0 1 


11 


The transformation matrix is given by 




or, in terms of the planar rotation matrices, the intermediate axes are elim¬ 
inated by 


X 


1 

0 

0 


Cos *' 

Sin V' 

0 


Cos ©' 

0 

-Sin ©• 

y 

- 

0 

Coo 0* 

Sin 0' 


-Sin *' 

Cos *• 

0 


0 

1 

1 

0 1 

z 


0 

-Sin 0* 

Cos 0' 


0 

0 

1 


Sin O' 

0 

Cos O' | 


The direction cosine elements of the transformation matrix nro obtuined by 
performing the Indicated nultipl 1 cat ion. For the pltch-y.iw-roll rotational 
sequence 


(Cos ♦' Cos O') 

(Sin 0' Sin O' 

-Cos 0 • Sin *' Cos O'? 

{Cos 0 ' Sin &' 

* Cos O' Sir. *' Sin y ') 


(Sin *•) 

(Cos 0' Cos *') 

( -Sin • Cos *•) 


(-Cos *' Sin 0) 

(Sin 0• Cos O' 

+ Cos 0* Sin O' Sin ♦') 

(Cos 0' C ■■ D' 

- Sin ' ' Sin O' Sin ♦') 


I.' 


( 48 ) 
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The transformation matrix is given by 


x 

y c 

z 




X 

y 

z 


or, in terms of the planar rotation matrices, the intermediate axes are elim¬ 
inated by 


X 


Cos y" 

Sin V" 

0 


1 

0 

0 

Coe 0" 

0 

-Sin e" 

y 

K 

-Sin y " 

Cos y" 

0 


0 

Cos (T 

Sin J?" 

0 

1 

0 

z 


0 

0 

1 


0 

-Sin f 

Cos 

Sin ©" 

0 

Cos e" 


The direction cosine elements of the transforiJJit 1 or: matrix are obtained by 
performing the Indicated mullipllcation. For the plteh-roll-ynw rotational 
sequence 


X 


(Cob y" Cos 0" 

+ Sin ♦ " Sin <f" Sin 0") 

(Sin *" Cos ^") 

(-Cos *" Sin ©" 

i Sin *” Sin Cos 0") 


y 

- 

(-Sin ♦" Cos 0" 

+ Cos ♦" Sin Sin 0”) 

(Cos y" Cos •f") 

(Sin 0" Sin y" 

♦ Cos y" Sin <f" Cos ©") 


% 


(Cos 0" S) , ©") 

(-Sin 0") 

(Cos Cos 0") 
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The directk ^.ines relating body and inertial coordinates are 
assigned the following symbols in the central program (see Equation 22). 


X 


4 

4 

h 


0 

y 


mi 

K2 



Y, 

a 

2 


n i 

n 2 

n 3 


K 

0' 


(50) 


By comparing identical positions in the matrix of Equation 50 with the 
matrices in Equations 47, 48, or 49, the direction cosines above are 
defined in terms of the appropriate sequence of Euler angles. 


(4) Platform Angles for a Flat-Planet Problem 

For a flat-planet problem, the orientation of the platform 
coordinate system will be assumed to coincide with the flat-planet 
coordinates. Therefore, the angles measured on the gimbals of this 
platform may be determined for the three gimbal arrangements considered. 
For the yaw-pitch-roll gimbal system, the following direction cosine 
relationships are obtained by comparing corresponding positions in the 


ma trice 

s used in Equations 

47 and 50. 

Five 

elements 

are sufficii 

ant 

to 

define 

these angles. 

4* 










-Sin 

a 








4 - 

Con 0 

Sin 

♦ 








Cos 0 

Cos 

♦ 




(51) 



m 3 = 

Sin (j 

Cos 

a 







n 3 " 

Cos y 

Coe. 

a 





The fir 

st equation defines 

the angle 
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angles 

and t may 

be 


defined 

explicitly by comb 
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:ond l 

fid 

thi rd 

equations an. 

> the 

fourth 
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For the flat-planet problem with the platform stabilized to 

coincide with the X -Y -Z coordinates, these angles represent the 
9 9 9 

angles measured on the gimbals and will be designated with a subscript p. 


*P 

*P 


-Sin"! 

Tan“l 4/4 
Tan -1 


(52) 


m 3/ n 3 


Similarly, the angles measured on a pitch-yaw-roll gimbal arrangement 
may be computed by comparing identical positions in the matrices used in 
Equations 48 and 50. 


Then 


merit are: 


h 

= 

Sin 

*P 

t. 

= 

Cos 

y * Cos G f 
p p 

4 

= 

-Cos y* Sin 0* 
P P 

m 2 

s: 

Coo 

tp Cor. V' 

n 2 


-Sin 4' Cos y* 
p 

Sin 

♦p 

c 

4 

Tan 

°P 



Tan 

<P 

rr 

-n^/nip 

problem 

, the gimbal , 

♦p 


Sin 

- 1 4 

°P 

w 


- 1 - 4/4 



Ten 

- 1 -np/r.'. ? 


(53) 


(54) 
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The appropriate direction cosines for the computation of the angles 
for a pitch-rol1-yaw system are: 

n2 = -Sin 0p 

mp = Cos yl! Cos 

(55) 

i 2 = sin yp cos 0p 

nx = Cos Sin ©£ 

n^ = Cos <jf" Cos ©p 

The platform angles are found from these direction cosines to be: 

0p ' -fxn' 1 n 2 

©p = Tan -1 nx/n^ (56) 

♦p = Tan' 1 l 2 /in 2 

For the flat-planet problem, the angles derived in Equations 51 

through 56 represent the attitudes of the vehicle with respect to the 
Xg'Yg~ Z g flat 'P^ anet coordinates and also with respect to a platform 

coordinate system whose respective X Y -Z axes are parallel to X Y Z . 

PPP 999 

The computation of these platform relations are summarized, along with 

the accelerometer indication in Figure 10. 

(5) Accelerometer Indicati .>. Let A be the vector sum of the 
platform accelerometer outputs and g be the mass attractive acceleration 
of the planet. The accelerometers are calibrated to read zero when they 
are unaccelerated and aligned such that the sensitive axis is 
perpendicular to g. The vector R will represent the displacement of 
the platform with respect to the inertial axes. It will be shown that 
A - R - g, wnere R is the inertial acceleration of the origin of the 
platform. Consider the vehicle accelerating vertically at lg. In 
inertial components, then 

R • -9 (5f) 
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Functional Flow Diagram - Platform Angles for 
Si x-De'ree-of-Freedom Flat-Planet Option 


Figure 10. 
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In the absence of a gravitational field, the accelerometer should read 
-Ig. Positive motion of the accelerometer mass along the Zg axis 
represents a negative acceleration in this case, and the vector g is 
equal to g TZg. Consideration of the gravitational field will cause an 
additional displacement of the accelerometer mass in the positive Zg 
direction, giving a total indication of -2g. The equation 

A * R - g (58) 

will be evaluated from the data 



so that 

A » - 2* 7 . (59) 

9 


This result is shown schematically in Figure 11. 



Figure 11. Accelerometer With Sensitive Axis Aligned With 
Local-Geocentric Vertical 
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The vector A is equal to the vector sum of the accelerations 

produced by the externally applied forces. The body components of the 

externally applied forces may be taken from the separate subprogram 

which aives the summation of forces and moments. F , F , and F are 

x y ^ 

the body components of the external forces plus the weight. The weight 
must then be substracted to determine the body components of A: 


A 


R-g 



F y -%g 




( 60 ) 

o x T x 4-o y Ty + a*l* 


3. VEHICLE CHARACTERISTICS 

The methods by which the aerodynamic, propulsive, and physical 
characteristics of a vehicle are introduced into the Takeoff and 
Landing Analysis computer program are presented in this section. The 
form and preparation of these input data are discussed together with 
methods by which stages and staging may be used to increase the 
effective data storage area allotted to a description of the vehicle's 
properties. 


a. Aerodynamic Coefficients 
(1) Form of Data Input 

The primary objective of the aerodynamic data input 
subprogram is to provide for a complete accounting of the various 
contributions to the aerodynamic forces and moments regardless of the 
flight conditions or the vehicle. !wo powerful techniques are 
available for use in digital computer programs: (a) an n-dimensional 
table look-up and interpolation and (b) an m-order polynomial function 
of n variables prepared by "curve fit" techniques. In the first 
method, the proper value for eai h term is obtained by an interpolation 
in "n" dimensions where the number of dimensions is taken to be the 
number of parameters tn bt varied independently plus the dependent 
variable. This method I, >> the advantage of accurately describing 
even the most nonlinear variations with a minimum of preparation effort. 
The amount of storage space whirl must be allocated to such a method, 
however, can achieve completely unrea enable proportions and may require 
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substantial computing time for interpolation as the number of dimensions 
is increased. The second method has essentially the opposite 
characteristics; that is, a large amount of data may be represented 
with a minimum amount of storage space and the computation time is 
held to reasonable limit? but the data variations which may be 
represented must be regular. A substantial amount of effort is usually 
required for the preparation of data by a curve-fit technique. Both of 
these methods are very convenient when the amount of data to be handled 
is moderate, but tend to become unmanageable when large amounts of data 
are required. This usually occurs when the program having several 
degrees of freedom, is committed to one or the other of these two 
techniques. Therefore, the Takeoff and Landing Analysis computer 
program will incorporate both of the techniques discussed as a compromise 
to take advantage of the more desirable features of both. To do this, 
a general set of data equations will be programmed which define each of 
the aerodynamic forces or moments. In general, the coefficients for 
these equations will be obtained from a curve-read interpolation. 

Several simplifications may be made to the equations, depending on the 
flight condition and vehicle to be considered. 

The effects of the following parameters will be considered: 

(a) Angle of attack and its time derivative (a , u) 

(b) Angle of sideslip and its time derivative (f, i-:) 

(c) Roll, pitch, and yaw control deflections (<S , A , I’y) 

(d) Rull, pitch, and yaw angular rates (p, q, r) 

(e) Center-of-gravi ty position (X^. ^ ) 
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The aerodynamic forces and moments considered with respect to each 
coordinate axis include the effects of angle of attack and sideslip, 
primary control deflection with respect to each axis, lag of downwash, 
and primary damping effects. In addition, the rolling moment due to 
yaw rate is included. Complete generality in the aerodynamic coupling 
effects has not been included in the present subprogram options since 
the descriptive terms required depend upon the particular problem. 

The storage space provided for the several existing options is considered 
to be adequate to accommodate other special problem formulations through 
substitution of terms. 

Quite often the particular application will not require Some of the 
terms listed in order to describe completely the flight path and vehicle 
under consideration. The subprogram will be arranged so that the 
computer will assign . constant value to any curve for which the data 
has not been supplied. For most curves, the constant value will be 
zero. This technique will reduce substantially the time required for 
the preparation of data. Values intermediate to those introduced in a 
tabular listing will be obtained by linear interpolation. The method 
of incorporating data for staged vehicles is discussed in the paragraph 
entitled Stages and Staging . 

(2) Flight Path and Vehicle Types 

In most of the cases discussed below, e "curve-fit" 
technique will be used to obtain all or a portion of the aerodynamic 
terms. For the purpose of this subprogram, it will be assumed that 
the curve fit has been selected to represent the variation of the 
coefficient about the trim conditions. This may have the effect of 
removing physical significance from some of the individual terms, and 
only the sum of the terms will represent the data. A typical example 
is Figure 12. 
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- Actual 

-— —- - Curva Fit 


Figure 12. Curve Fit Nonlinear Aerodynamic Characteristics 

In this case, the C^ and C^ values used in the equation for C are 

o a 

obviously different from the actual values of these parameters. 

A functional flow diagram for the solution ol the aerodynamic 
forces and moments is presented in Figure 13. It should be noted that 
the actual machine programming will not necessarily follow the sequence 
shown since certain computer operations have been omitted in this 
description of the problem formulation. 

A controlled aircraft represents the most general case that will be 
considered. In order to account for the many component forces, we made 
certain restricting assumptions. The assumptions will be made that the 
aircraft is confined to moderate variations in position angles and 
control deflections. Center-of-gravity shift along the x-axis is 
included. 'he coefficients can then be expressed as shown in Block 3 
of Figure 13. The functional computation sequence for this option 
proceeds from Block 1 to Block 5 in a straightforward manner. 

The change if. dynamic derivatives due to a change in the center-of- 
gravity location is programmed as a curve-read in order to avoid the 
complications of a transfer. The definition of u and ft as applied to 
the TOLA computer program (see Equations 45 and 46) is 

a ■ Ton ~' ) and P * Tan ~' 160 


48 



AFFDL-TR-71-155 
Part II 



Figure 13. Solution of Aerodynamic Forces and Moments Subprogram 
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Data supplied must correspond to this definition or an alternate 
computation of these angles must be formulated to agree with the method 
of data reduction. 

(3) Error Constants 

The use of error constants, designated by the symbol , 
to modify the aerodynamic data characteristics is shown in Figure 13. 

A detailed explanation of these error constants and their use is given 
later. 


b. Thrust Fuel Flow Data 

An airbreathing engine is. strongly affected by the environmental 
conditions under which it is operating. Engines which would be grouped 
in this classification are turbojets, ramjets, pulsejets, turboprops, 
and reciprocating machines. The parameters which will be considered 
of consequence in this program are: 

(a) Mach Number (M^) 

(b) Throttle setting (N - units defined by problem). 

To accommodate these variables, a three-dimensional tabular listing 
and interpolation will be used to obtain thrust. The thrust needs no 
further correction since the effects of all paiameters are included in 
the interpolated value. The functional computation sequence tor 
introducing these data is straightforward, as outlined in Figure 14. 

The same thrust table is used for each of i enyines, where i - 1, 
.... IN. Each throttle may be varied independently by the Throttle 
Autopilot. All thrust vectors are assumed to act:: 

IN 

T* • I Tj IN,, M n ) (62) 

1*1 


bO 





r 


AFFDL-TR-71-155 
Part II 



Figure 14. Thrust Forces and Moments Subprogram 


where T x is the x- component of thrust in the vehicle body-axes system. 
(A positive T produces a positive u.) This force will introduce 


moments: 


IN IN 

m T * X W Ti * 2, T i z Ni 

i»l <=1 

IN IN 

N T * Z N Ti * I " T i VNi 

i«l i*l 


(63) 


where My ami are the thrust moments about the vehicle y and z body 
axes, respectively, and y^ and z^. are the distances to the center 
of gravity of the i tn enyine from the reference center of gravity, 
where i = 1, 2, ... , IN. 


c. Physical Characteristics 

The methods to be employed for the introduction of vehicle 
physical characteristics i. m tne lakeoti and La.iui.tg Analysis computer 
program are outlined in this section. A table look-up and interpolation 
technique is used to determine those parameters which are variable. 

A provision is made for introducimj ercn. consi'ants ..Lc- several of th j 
parameters. 
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(1) Categories of Physical Characteristics 

Physical characteristics are introduced into the computer 
program in two groups: (a) characteristics used in the general 
solution of the equations of motion, and (b) characteristics used 
only in specific, or auxiliary, subprograms. The following items 
will be defined in the general vehicle characteristics subprogram: 

(c) Initial mass of the vehicle [V / 0 ), 

(b) Reference area (S), 


(c) Reference lengths (d,, d^). 


(d) Reference center-of-gravity location (X c ^ R e f)» 

(e) Rotating machinery pitch angle ( 9 r ), 

(f) Rotating machinery angular rate 


(g) Rotating machinery moments of inertia (I , I , I ), 

x r Y r z r 


(h) Vehicle center-of-gravity location (x ), 


I 


xz ’ 



(i) Vehicle moments of Inertia (I xx , I , I 2? , I , 
and 


(j) Reference jet-damping lengths (1 , 1^, 1^, l m> l n ). 

Items (a) through (g) wil 1 be constant throughout any stage. Items 
(h) through (j) will be variable during the stage due to the variation 
in mass caused by fuel consumption. Figure 15 presents a functional 
flow diagram defining the manner in which these characteristics are 
introduced into the compute 1 " program. 


5 ? 
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( 2 ) Reference Weight. The instantaneous mass is used in the 
computation of the body motion. The reference weight is obtained by: 

W T * W ( 32.174) (64) 

(3) Error Constants. The use of error constants, designated 
by the symbol c., to modify the general ''ehicle physical characteristics 
is shown in Figure 15. A detailed explanation of these error constants 
and their use is given later. 

d. Stages and Staging 

A problem encountered frequently in airplane performance work 
is that of staging to allow a set of aerodynamic data to be changed. 

Stage changes at constant weight, such as extending drag brakes or 
turning on afterburners, may also require revising the aerodynamic or 
physical characteristics of the vehicle. Another use of the staging 
technique is possible with the present computer program which does not 
involve physical changes to the configuration; this technique may be 
used to revise the aerodynamic descriptors as a function of aerodynamic 
attitude or Mach number. With this use of the stage concept, accurate 
descriptions of the forces and moments acting upon a vehicle may be 
maintained over wide attitude ranges if required. Other applications 
of this stage technique are possible. Normally it is not practical to 
stop the computer and manually insert a new set of data. A better 
approach is to have the computer do this automatically: new data will 
be loaded automatically on the basis of whether a specified variable 
has exceeded or become less than a preselected value. For generality, 
it is possible to test on four values in each direction. 

When the new data are read in, the conditions representing the 
last time step will be read in as initial conditions for the next stage. 
This avoids the discontinuity that would result from an infinite rate 
of change of center-of-gravity location. It also will cause the 
integration routine to be started over, which will reduce the computer- 
induced transients due to staging. 
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e. Error Analyses 

The Takeoff and Landing Analysis computer program will incorporate 
a provision for conveniently performing flight-path error and dispersion 
analyses by trajectory computation. This prob’em involves the 
determination of flight-path dispersion due to deviations of input 
quantities from their predicted nominal values. The usual approach to 
this type of problem requires that a series of trajectories be computed 
in which standard deviations, or errors, are systematically introduced 
for each parameter while the remaining parameters are held at their 
nominal values. These results are then combined to determine the 
"probable" dispersion. This approach will be implemented in the 
Takeoff and Landing Analysis computer program by providing a simple and 
efficient method of introducing the deviations. The capability of 
modifying a nominal value by either an error constant multiplier or an 
additive error constant is provided for many of the parameters as 
outlined below. The provision of these error constants will reduce 
substantially the number of tabular data listings that must be changed 
for an error analysis, thereby reducing the work of the analyst. The 
determine!I ion of the standard deviation of each of the parameters and 
the method of combining the trajectory variations are left to the analyst 
in view of the multiplicity of combinations possible. 

(1) Aerodynamic Data. The provision to modify the aerodynamic 
coefficients through the use of error constants, e-> is outlined in 
Figure 13. The constants are applied as follows: 

h F » < «, C N + « t )q*S 

q * < « 5 C A + < 4 ) q* S 

y * I« B Cy+«*)q*S 
I 8 + 

m 8 ( e, C m + « (0 ) q* Sd, 

n ■ < C n «„ 1 q* Sd f 


55 







AFFDL-TR-71-155 
Part II 


These error constants allow the total aerodynamic coefficient to be 
modified to account for configuration modification, experimental or 
analytical error, or misalignments. 

(2) Vehicle Physical Characterisitcs 

The provision to modify some of the vehicle physical 
characteristics through the use of error constants is outlined in 
Figure 15. The constants are applied as follows: 


Cb 

- 

f 


s 

f m +«, 9 

1 y y 

= 

f m + e t0 


= 

f m+c t , 

*:<y 

X 

f m + e tz 

Ixz 

X 

f tw)+« fs 

Iyz 

= 

f m) € t4 


In the application of error constants in the above equations caution 
must be exercised to insure that the units are consistent. Each of the 
error constants will he assigned a nominal value which will be used 
when no other value is specified. The constants which are multipliers 
will have a nominal value of unity, and those that are additive will 
have a nominal value of zero. 

(3) Additional Errors. Not all of the system input constants 
can be modified for error analysis studies as indicated above. In 
certain cases, it may be found unrealistic to modify the input data 
through the use of error constants because the actual deviation would 
not appear as simply a constant increment or percentag 0 change. An 
example of such a case would be the change in thrust-time history of a 
rocket due to temperature changes of the propellant since such a change 
affects both thrust level and burning time. For an accurate 
representation of such a case, it would be necessary to modify trie entire 
tabular listing accordingly. 
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4. VEHICLE ENVIRONMENT 

The models for simulating the environment in which a vehicle will 
operate are presented in this section. This environment includes the 
atmospheric wind. In the discussions which follow, the descriotions 
of vehicle environment pertain to the planet Earth. The environmental 
simulation may be extended to any planet by replacing appropriate 
constants in the describing equations. 

a. Atmospheres 

The concept of a modal atmosphere was introduced many years 
ago, and over the years several models have been developed. Reference 
16 outlines the historical background of the gradual evolution of the 
ARDC model. The original (1956) ARDC model has been revised to reflect 
the density variation with altitude that was obtained from an analysis 
of artificial satellite orbit data. This revision is the 1959 ARDC 
Model Atmosphere. 

The advantage of a model atmosphere is that it provides a common 
reference upon which performance calculations can be based. The 
model is not intended to be the "final word" on the properties of the 
atmosphere for a particular time and location. II must be realized 
that the properties of the atmosphere are quite variable and are affected 
by many parameters other than altitude. At the present time, the 
"state-of-the-art" is not advanced to the point where these parameters 
can be accounted for, and it may be several ypars before the effects 
of some parameters can be evaluated. 

(1) 1959 ARDC Model Atmosphere. 

The 1959 ARDC Mouel Atmosphere is SDecified in layers 
assuming either isothermal or linear temperature lapse-rate sections. 

This construction makes it very convenient to incorporate other 
atmospheres, either from specifications for design purposes or for 
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other planets. The relations which mathematically specify the 1959 
ARDC Model Atmosphere are as follows (Reference 17 ): 


The 1959 

ARDC Model Atmosphere is divided into 

11 layers as follows 

Laver 

Hh -Lover Altitude 


Upper Altitude 


(Geopotential) 


(Geopotential) 


Meters 


Meters 

1 

0 


11,000 

2 

11,000 


25,000 

3 

25,000 


47,000 

4 

47,000 


53,000 

5 

53,000 


79,000 

6 

79,000 


90,000 

7 

90,000 


105,000 

8 

105,000 


160,000 

9 

160,000 


170,000 

10 

170,000 


200,000 

ll 

200,000 


700,000 

For 1 ayers I, 

3 , 5 , 7 , 8, 9 , 10 , and 11 , 

a linear 

molecular-scale 

temperature lapse-rate is assumed and the 

following 

equations are used: 


H gp 


0.3048 h 

meters 

(65) 


1 + 0.3048/6356766 

T M 

S 

|T M : b['* K l IH «p- H b»] 

°R 

(66) 

T 

■ 

t„[a-b ,„b'( gc D )] 

°R 

167) 

P 

= 

P b [' + K,(Hg p -H b )]' K 2 

lb/ft 2 

(68) 

P 

s 

P b [l + K,( Hgp - 

slugs/ft 9 

(69) 

v» 

9 

49.020576(T M )''* 

ft /sec 

(70) 

V 

S 

0.0226988 X 10 * [ (T+ ,1,, 72 ,, ] 

ft*/sec 

(71) 
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For the isothermal layers 2 , 4 , and 6, the following changes are made 
in the above equations: 


P 

P 


Pb t 
6 


*Kg(Hgp ~ Hp) 

(72) 

-Kg(Hgp-Hb) 

(73) 


Values for temperature, pressure, density, and altitude at the base 
of each altitude layer are listed below, along with the appropriate 
values of , K^, and K^. 


Quantity 

1 

2 

3 4 

5 6- 

Ki 

-.225569 Jt 

0 

.138466-’* 0 

-.159202-4 0 

k 2 

-5.25612 

- 

11.3883 

-7.59218 

k 3 

- 

.157639*3 

.120869-3 - .206234-3 

*b 

518.638 

389.988 

389.988 508.788 508.788 298.186 

Pb 

2116.21695 

472.73 

51.979 2-5155 

1.2181 2.1080" 2 

Pb 

2.37692-3 

7.0620- 1 * 

7.7650-5 2.8304 

"6 1.39463-6 4.1189- 3 

»b 

0 

11000. 

25000. 47000. 

53000. 79000. 

Quantity 

7 

8 

9 

10 11 

Kl 

.241458-4 

.886289- 

■4 .754341-5 

.350715-5 .222129-5 

k 2 

8.54.120 

1.70824 

3.41643 

6.83296 9.7613'/ 

k 3 

- 

- 

- 

- 

1'b 

298.188 

4 o 6 .l 83 

2386.188 

2566.188 2036.133 

Pb 

?..l 809'3 

1.5562-’ 1 

7.5578-6 

5.8954-6 2.9759-6 

P b 

4.261-9 

2.232" 10 

].845-13 

1.333" 1? 6.113 _1 3 

Hb 

90000. 

105000. 

160000. 

170000. 200000. 
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Values of the appropriate constants to be applied 
equatior (Equation 67) are as follows: 

in the 

temperature 

Hgp(Km) 

A 

B 

C 

D 

0-90 

1 . 

0 . 

- 

- 

90-1S0 

0-759511 

0.17U16U 

220 

25 

180-1200 

0.935787 

0.273966 

ISO 

lho 


(2) Limitations. The validity of the 1959 ARDC model is 
limited to altitudes below 700 km, although the program is arranged to 
extrapolate the relationships to greater altitudes, if desired. 
Extrapolation to greater altitudes is accomplished by altering the cutoff 
altitude. At an altitude of 90 km (approximately 300,000 ft.) the 
subprogram normally ceases to calculate kinematic viscosity and speed of 
sound and assigns a value of zero to each of these parameters as an 
indication that the computation has stopped. This is done for the 
following reasons: (a) the molecular composition of the atmosphere is 
unknown; (b) the variation of the ratio of specific heats above 90 km, 
is not known; and (c) the numerical value of the speed of sound has 
little physical significance. The validity of Sutherland's empirical 
formula for viscosity is also reduced because of the extremely low 
pressures which exist. 

(3) Accuracy. Due to a lack of knowledge of the rounding-off 
procedures used to evaluate the constants in Reference 16, it was 
impossible to obtain exact agreement between the subprogram and the 
values tabulated there. A comparison of the results over an altitude 
range of 0 - 1,000,000 ft. revealed that the deviation of the computed 
from the reference values never exceeded one tenth of one percent and 
in most cases was less than one half of this value. 
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b. Winds Aloft 

The winds-a'ioft subprogram provides for three separate methods 
of introducing the wind vector - as a function of altitude, a function 
of range, and a function of time. This will facilitate the 
investigation of winds effects for the conventional performance studies. 
The wind vector \ ill be approximated by a series of straight-line 
segments for each of the methods mentioned above. Statistically derived 
profiles of the type presented in Reference 13 can be represented by 
this approach and it is presumed that the analyst will resort to sources 
of this type to obtain the wind irput data. The present subprogram 
will not be particularly concerned with the method used to determine the 
wind vector, as this is a separate problem outside the scope of the 
Takeoff and Landing Analysis computer program. 

Four options will be used to define the wind vector in the SDF 
computer program. The three components of the wind vector in the 
inertial coordinate system will be specified as tabular listings with 
linear interpolations (curve reads) in the following options. 

Wind Option (0) . In this option the wind vector is zero throughout 
the problem. This will allow the analyst the option of evaluating 
performance without the effects of wind. This option causes the wind- 
aloft subprogram to be bypassed in the computational sequence. 

Wind Option (1) . In this option the components oF the wind vector 
will be specified as a function of time for the estimated altitude. 

Wind speed will be specified in f eot per second and time will be 
specified in seconds. 

Wind Option ( 2 ) . The three components of the wind vector will be 
introduced as a function of altitude in this option. Wind speed will 
ba specified in feet per second ariu altitude will be specified in feet. 
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Wind Option (3) . In this option the components of the wind 
vector will be introduced as ? function of range for the estimated 
altitude. Wind speed will be specified in fec f per second and range 
will be specified in nautical miles. The range used in this computation 
will be the range from the starting poi"t, Rg. 

By staging of the wind option, it w'11 be possible to switch from 
one method of reading wind data to another during the computer run. 

Care must be exercised in this operation, however, as the switching 
will introduce sharp-edged gusts if there are sizeable differences in 
the wind vector from one option to another at the time of switching. 

This effect should be avoided except in cases where gust effects are 
being studied. 

Figure 16 presents a functional flow diagram of the wind-aloft 
subprogram. For Wind Option (3) the range from the starting point must 
be defined. Therefore 

Rg * 'v/* X g-Xg 0 )2 + (Yg-y g Q ) £ (74) 

5. SUF-2 CHANGES 

Nearly all of the autopilot logic requires additions to >.DF-2 anJ 
not changes in the basic 5DF-2 cuncepts. Because it is desired to 
simulate ground effect, multiple engines, and a drag chute, some 
direct changes in SDF-2 must be m, le. 

Aerodynamic" (vv'S) Cnanges 

The aerodynamics subroutine in SDF-2 (called SACS) will not 
simulate a smooth transition into ground effect and has several options 
that take up needed computer storage and are not needed in a landing 
analysis. Aerodynamic dependence on Mach number is also no longer 
needed (see Assumptions in Section II of Appendix III). Special 
requiiernents of the autopilot (i.e. aerodynamic force coefficients 
and Cp are needed in a pseudo-wind axes system, and the side Forte 
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coefficient, C , and the- aerodynamic moment coefficients are needed 
in the body axes system) can a iso L met at the outset by a rewrite 
of SACS. 


Ground effect data, if available at all, is usually given by 
giving the aerodynamic coefficients for no ground effect (i.e., 
aircraft is sufficiently away from the ground) and by giving the 
aerodynamic coefficients in ground effect (i.e., aircraft is essential 
on the ground). References 7 and 19 indicate that the aerodynamic 
coefficients affected begin to change when the aircraft reaches one wing 
span above the ground, and the transition into full ground effect 
appears t.o be exponential. For this reason, the following equation 
was adopted to simulate the ground effect transition: 


/ h G -« U R 2 >-» C 6 ^ 

A hG * A o + < A GF" A o )e V b ‘ h CG ' 


(75) 


where 


f(x R2> 




= runway perturbation at the displacement of gear number 
two from the runway origin (can represent an aircraft 
carrier deck - for takeoff) 

= aerodynamic coefficient at altitude h Q 

= aerodynamic coefficient for no ground effect 

= aerodynamic coefficient for full ground effect 

= altitude aircraft mass center is above the ground 

- sec "glide slope" 

- wing span 


Note that for > b, the value of A^ is essentially A q . As hg 
approaches h , the value of A^g is essentially Ag P . In this manner 
the ground effect transition is simulated. 
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The details of the rewrite of SACS will not be given here. Suffice 
it to say that the modified SACS subroutine will include the following: 

(1) Same equations and aerodynamics coefficients used in 
Option I (page 45) 

(2) Eliminate all other SACS Options and the aerothermoelastic 
calculation. 

(3) Aerodynamic coefficients will not be a function of Mach 
number anji two data points are needed for each coefficient: one for no 
ground effect and one for full ground effect. 

(4) Force coefficients C L and C^ must be input in the pseudo 
wind axes. 

(5) Side force coefficient, C , and the moment coefficients 
must be input in body axes. 

b. Multiple Engine Changes 

As originally designed, SDF-2 only allowed for the simulation 
of one airbreathing engine. The thrust table look-up routine, TFFS, 
included the effects of altitude on thrust and determined the fuel rate, 
both of which are no longer needed (see Assumption in Section II of 
Appendix III). Some provision must also be made to determine thrust 
for engine reversing, thrust for engine failure, and the net thrust 
forces and moments acting on the multiple engine aircraft. We begin 
by discussing the thrust table lookup. 

The thrust table is now a function of throttle setting and Mach 
number alone. Because of the throttle autopilot logic, the ranges for 
the throttle setting, N, must be as follows: 

N = -2 means f ulI reverse 

N = -1 means idle reverse 
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N = 0 means engine failure 
N = +1 means idle forward 

N = +2 means full forward 


Negative values of thrust are stored for the negative throttle settings. 
In this manner reverse thrust is simulated. The zero value of N is 
used for the data simulating engine failure (note that failure thrust 
may be zero or some negative value, depending on Mach number). The 
actual reversing is achieved in the throttle autopilot by requesting 
a negative value of desired throttle setting, N ... The change of sign 
on N is sensed in the control response (see Figure 29 of Appendix III) 
and this changes the sign of the actual N. used in the thrust table. 

Every time thrust is needed by the main program, the thrust table is 

entered IN times - once for each engine - and the actual thrust array 

T(IN) is obtained depending on the values in the actual throttle setting 

array N (IN). The engine thrust vector is assumed parallel to the 

longitudinal body axis (this is also assumed in the autopilot equations) 

and therefore causes no roll moments. The engine position arrays, 

3^ (IN) and Y (IN), along with the actual thrust array T (IN) 

determine the engine pitch moment array M-j-(IN) and the engine yaw 

moment array Ny(IN). These arrays are then used to obtain the net 

longitudinal thrust, T , (note T and T~ are zero by assumption) 

x y o 

the net engine pitch moment, My, and the net engine yaw moment, Ny. 

The net values T x> My, and Ny are then used in SDF-2 and the calculation 
proceeds as normal. 


c. Drag Chute Changes 


A drag chute simulation is included for analysis of those 
aircraft using chutes for deceleration on landing roll. The chute drag 
vector force, F^, is assumed opposite the relative airspeed vector 
’’ and is written as follows: 



-F 


DC T a 


(76) 
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where 

F DC " C DCH S SH 


and 



= magnitude of chute drag force 
= chute drag coefficient (assumed a constant) 

= chute reference area 
= dynamic pressure (see SDF-2) 


( 77 ) 


In terms of SDF-2 variables, the relative airspeed vector in body 
axes T x , T , ~T is 

V a « (u -U w ) T x + ( v-v w ) T y + (w-w w )T z (78) 

Equation 76 can therefore be written 

F OC * ~vp [ (u - u w , *X + iv -''w>»y + ( "'- w w> T z] (79) 

The vector, , from the aircraft mass center to the chute attachment 
point is written in body a v es as 

^CH * X CH *x + y CH + Z CM 

The vector moment, M CH , of F D( . about the aircraft mass center is 

M CH * ^CH * F OC 


If F^ is written 


F 0C * F cx *x ^cy *y + F cz *z 


(82) 


and M ch 


is written as. 


^CH " M Ox '* + M Oy *y * **Dz *7 


(83) 
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then the scalar components F cx , F^, F c3 , M DX> M Dy and M are as 
follows: 

F cx s — F( u - u w ) / Vg (84) 

F cy s -F dC (w-v w )/V Q '85) 

F cz “ -F oc (w-w w )/V 0 (86) 

^Dx * Y CH F cz ~ *CH F cy (® 7 ) 

**^Dy " ^CH F cx -X CH F cz I® 0 ) 

M Dz * X CH F cy “ y CH F cx l® 9 ) 

Figure 17 is a flow chart of the drag chute equations. The ICS 

indicator is normally input at a value other than one. As such, the 

chute body forces and moments (i.e., Equations 84-89) are all zero. 

The drag chute is deployed some time during the landing rollout (see 
Figure 11 of Apendix III), at which time the ICS indicator is made one, 
which indicates to the program that the drag chute is deployed. The 
chute body forces and moments are finally used in the SDF-? total 
summation of body forces and moments (i.e., F x> F , F , L, M, N) acting 
on the aircraft. 
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' IF \ 
ICS = 1 


F cx = 

0 

F cy = 

0 

F cz = 

0 

M Dx = 

0 

M D y = 

0 

M 0z = 

0 


f dc 

= C gen s 

SUP* 

F cx 

= -f dc ( u - 

U W )/V Q 

F cy 

-f dc <v- 

< 

* 

X 

< 

o 

F cz 

= -FdcIW- 

o 

> 

* 

F dx 

= y ch f cz 

-Zch f CY 

Mjjy 

= 2ch f CX 

" x CH f CZ 

Mdz 

= x ch f cy 

-y ch f cx 


SACS COMPUTATION 


CONTINUE CALCULATION 
AND ADD Fcx.Fcy.Fcz, 
Md*.M dy ,M dz TO BODY 
FORCE SUMMATION 


Figure 17. Drag Chute Forces and Moments 
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APPENDIX II 

LANDING GEAR EQUATIONS 
TO SDF-2 
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SECTION I 
INTRODUCTION 


This appendix shows how the present equations of motion in SD T -2 
must be modified to include the landing gear dynamics and ground 
reactions. In this formulation the main aircraft frame is assumed a 
rigid body; however, the gears are allowed to move, and the dynamics 
of this motion is included in the formulation. The formulation is 

generalized to a vehicle with any number of oleo struts (single or 
double chamber) with balloon tires. Though specific application is 
first intended for the C-5A aircraft, a modification of SDF-2 will 
allow takeoff and landing analysis of many aircraft. 


For a detailed understanding of the formulation, we reference the 
reader to two reports: Re f erence 1, which is the original documentation 
of SDF-2; and Reference 6, which is a derivation of the equations of 
motion for a series of nonrigid bodies. 
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SECTION II 
MODIFIED EQUATIONS 


1. SDF-2 

The equations of motion in Option 2 of the SDF are written in a body- 

fixed-axes system, T . T , 1 , which is located at the vehicle mass 

xo yo zo 

center. The mass center is located relative to an earth-fixed-axes 

system, 1^, 1^, 1^, O zg points in the direction of gravity, i.e., 

down) by a vector R which has coordinates X^, Y^, along the 

T xg , T yg , T zg axes, respectively. The body-fixed-axes T xQ , T yQ , T zq 

are oriented with respect to earth axes T , 7 ,1 through the 

xg yg zg 3 

j" 1 m J matrix of direction cosines. This is illustrated in Figure i8. 



Figure 18. SDF-2 Coordinate Systems 


As such, the displacement vector R can be written 

R • T xg +V yg + Z,i zg (90) 
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and the transformation matrix 
system is 


from the 1 


9 


axes 


system to the 1 axes 


r -i 


— “ 


“ — 

— i 
X 

o 




Ug 

lyo 

s 

m | ffi g nr* j 


1 y« 



n, n t n 3 


^zg 



- 


- 


(98) 


Before proceeding with deriving the modified equations of motion, 
it is instructive to examine the landing gear to be simulated. 

2. LANDING GEAR 

Reference 26 describes the C-5A landing gear. The basic gear 
configuration is shown in Figure 19. 



Upper Air Chamber 


Primary Metering Orifice 
Secondary Metering Orifice 

Secondary Piston 
Extension Stop 

Secondary Piston 

Secondary Piston 
Compression Stop 

Lower Air Chamber 


Figure 19. Strut Configuration 
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The strut configuration is unique in that it has a secondary piston 
and associated air chamber. Both the upper and lower air chambers are 
preloaded, with the lower chamber having the higher preload pressure. 

The secondary piston remains forced against its extension stop until the 
strut has deflected to the state where the upper chamber forces exceed 
the preload pressure force of the lower air chamber. Once this occurs, 
the secondary piston begins to move downward to seek an equilibrium 
position. With this design, the secondary piston allows continued 
strut movement without severely increasing the air compression load on 
the strut. Since most gear struts have only one air chamber, the final 
formulation will have an input indicator which will allow the secondary 
piston and lower air chamber to be removed from the problem. 

3. MODIFIED EQUATIONS OF MOTION 

a. Present SDF-2 Equations of Motion 

The present SDF-2 equations of motion (see Reference 1, pgs 5-7) 
are written with the mass center of the entire vehicle as the reference 
point and in a body-fixed-axes system. As such, the vector application 
of Newton's Law takes the form 

F ■ « m ( w + £3 x v ) (92) 


M - * tHT*) ■I w + S*lf C) (93) 

QT OT 

where 

F = total applied vector force 
m = total mass 

V = inertial velocity vector of mass center in body-fixed-axes 

system 

V = velocity vector rate as seen from the body-fixed-axes 

system 
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Hi ~ inertial rotation rate vector of body axes 
H = rigid body angular momentum about mass center 
I = moment of inertia matrix about mass center 
5 » = rate of change of SJ vector 


Equation 92 is the same as Equation 14 of Appendix I, and Equation 
93 is the same as Equation 17 of Appendix I. 

b. Modified Equations of Motion 
(1) Identification of Problem 

Although the problem definition of TOLA does not include 
body flexure as part of the problem, it does include the rigid body 
movement of the struts. Thus one has a problem where rigid bodies 
are in contact with and moving relative to one another. In the strict 
sense of the word, this is no longer a rigid body problem. 


The equations of motion for such a system have been derived in 
Reference 6. The general equations of motion (see Equations 48 and 
50 of Reference 6) are therefore: 


F T - I F k « m T R + £ {(R k ) Q + 25 0 x(S k ) 0 

k*o k»o 


+ x 


(ff h ) 0 + w 0 x [w 0 x{R k > 0 ] 


+ ,7 kc , U + 


(94) 


25) k* |r kc , k + u, k* 


? kc + 


u5 k x (o> k 


’■n.c 1 } 


M o * I "ho ‘ m T F c * ^ ♦ I m k [ ( Vo + 7 kc] * 
k*o k*o 

{(R k ) 0 + 2a 0 x(R K ) 0 -E £ 0 < (R k ) Q + ui 0 x [ffl 0 x(R k ) 0 ] 

(95) 

+ , \c , k + 2i> k xi7 kc ) k + i k* ? kc 


ui k x(a) k x 


he 


K 

+ z * 

h*o 


'k® 
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The derivation of Equations 94 and 95 is quite long and details may 
be found in Reference 6. The definition of each term in Equations 94 
arid 95 is, however, given here as a point of departure, and is 
illustrated in Figure 20. 


"h- 

R 


m 


k 


U) 

0 



r 


kc 


( 


CJ 


k 


M 

o 



r 

c 


H 


ko 


the total applied vector f orce to the system of K + 1 
bodies 

total applied vector force on the k*^ body 

a given body: k = 0 is the firs, body, and k K 
is the last body. 

total mass of K + 1 bodies 

th 

vector to reference point on (j body 
total mass of k^ body 

i. L_ 

inertial rotation rate vector of 0 tn body 
vector position of k tn body reference point from 

i L. 

(j body reference point 

vector position of x n body mass center from k L body 
reference point 

x L 

a symbol meaning "as seen by the k Ln coordinate 
system" 

inertial rotation rate vector of k n body axes 

total vector moment of all applied forces to the K 1 
bodies about the 0 ™ body reference point 

total ve<-co. moment of applied forces to body 
about oth body reference point 

vector position of mass center of the K + 1 bodies 
from otn body reference point 

vector moment of linear momentum of k th body about k th 
body mass center 


m 
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Equations 94 and 95 are very general and therefore very complicated, 
but have the added feature of being very complete and rigorous as regards 
their dynamics formulation At present. Equations 94 and 95 bear ver> 
little resemblance to the SDF-2 equations of motion (i.e., Equations 
92 and 93), and, at first examination, the necessary modification to 
Equations 92 and 93 may appear extreme! ' large. Specific application 
of Equations 94 and 95 will show this not to be true. We begin by 
giving physical meaning to the coordinate system in Figure 20. 

The \g’ ^yg* ^ zg axes system in Figure 20 is the same as the 
SDF-2 axes system in Figure 18 (i.e., the inertial axes system). 

The o th body is presently chosen to the vehicle structure less the 
moving gear struts. The other K bodies represent the moving gear 
struts. In reality, t.nere are two moving masses per strut, i.e., 
the main strut and secondary piston; the secondary piston mass, 
however, is so small in comparison to the main strut mass 
(approximately 1/50), its inertial effects on the overal 1 motion of the 
vehicle are neglected. 

Each strut mass has assigned to it a coordinate system 1 , 1 k, 1 , 

th _ _ xk y v. 

fixed to the o body (note that makesw. equal w , and convenient! 

a. u K 0 _ _ _ 

oriented with respect to the o u body coordinate system l xQ , 1 , 1 ^ q 

as shown in Figure 21. 


The strut coordinate system 1 . , 1 . , 1 , is located by the vector 
- ■ th xk yk zk 


Rjj (note since the o body is assumed rigid, this vector does not 
vary) and is rotated about the positive 1 ^ axis by the angle G 
The angle 9, is chosen so that the direction ot gear movement is along 


the unit vector 1 


zk - 


With these physical 
us now examine Equation 


meanings given to the coordinate systems, 
95. 


let 
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Figure 21. Strut Coordinate System 


(2) Application of Equation 95 

The terms of Equation 95 can be divided into two groups: 
terms which depend on position and terms which depend on rate of change 
of position. Equation 95 is first rearranged into this division. 

Mo - I [ m k[<Rk l o+ 7 kc] k{«o^R k ) 0 + 
k»o 

(S 0 * [a 0 x(&k> 0 ] + ‘‘'k* 7 kc + ®k*^k x 7 ltc*} + 

IkA + ^V'V] + £ [^k + V^k {96 ‘* 

k»o 

+ m ) ,[Ul k l 0 + 7 | ,c] »{(# k l 0 + JO 0 »l*»> 0 +lV» 

+ za,«ttjtin}] + «> T f c »« 
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4 


The first summation in Equation 96 includes the terms that depend on 
position, i.e., the rigid body contribution to rate of change of angular 

i.L. 

momentum about the o tn body reference point. The second summation is 
that which depends on the time derivatives of position, i.e., the 
nonrigid body contribution to rate of change of angular momentum about 
the o th body reference point. Note that has been replaced with its 
equivalent expression (Equation 37 Reference 6), i.e.. 


The entire first summation in Equation 96 can be replaced with 

I .to +to x (I . to where 1 is the familiar rigid body inertia 
tensor for the entire K + 1 bodies about the o tn body reference point. 
Note that l Q will vary in time ever so slightly because r^ is changing 
for each gear. The variation of is normally on the order of 10-12 
inches, as compared to vehicle dimensions on the order of hundreds of 
feet. The variation of r Q due to gear strut movement will therefore be 
extremely small and justifiably neglected. 

A. L 

We next choose to let. the reference point of the o body be the 
nominal mass center of the entire vehicle (i.e., the mass center for 
the gears in their average deflected position). As such, r c will be 
very small at any given time and can be considered approximately zero. 


Reference 6 shows that 1^ and 1^ . are identically zero since 
each mass is a rigid body within itself. The vectors (R^) and (R^) 

are also identically zero since the body is considered rigid. 


The foregoing specific application of Equation 96 has presently 
reduced Equation 96 to 


Me * ^o ■ J 

+ i B k K>o ^c] * { ( “kc , * + 2 C h l( *l«c , l«} 
h»o 


( 98 ) 
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Separating the gear struts from the main vehicle structure in the 
summation yields 

M 0 * f 0 • <5 0 + <S 0 x ( f o • tD 0 ) + m 0 [ ( R 0 ) 0 + 'oc ] * 


f(r oe ) 0 +2a) 0 x<F 0C ) 0 ] + I "* k [<R k > 0 + 'k C l x [ (7 kc , k (99) 

J k-1 


+ 2<D h 


x<7 &c>k] 


th 


The mid; le term in Equation 99 is identically zero since the o body is 

riqid, i.e., (r ) = ir ) =0. Equation 99 therefore becomes 

3 oc oc 

0 0 


M c * lo * ^ I 

K 


( 100 ) 


^ | m k[ ( ^»o + 7 kc]*[ (? kc ) k +2C k xt? kc , k] 


Realizing that the coriolis acceleration of the strut, 2u^, x (r kc ) , 

will be very small (mainly because the rotational velocity of the 
vehicle, ui^, at landing is very small) in comparison to the strut 

acceleration (r, ) , reduces Equation 100 to its final form 
Kc k 


Mo * Tq * o 0 +x a n ) 


0 0 1 


( 101 ) 


+ I w k [<*k>o+ 7 kc] x <*kc>k 
k»l 


Note that Equation 101 is the same as the SDF-2 equation (Equation 
93) except for the summation involving the struts. 



(3) Application of Equation 94 

We next examine Equation 94. Let us first realize that 
and (R^) are zero because the o tfl body is rigid. Next we 
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/ 


realize that (r ) and (r ) are zero for the same reason, 
o o 

— 4 - U 

(R ) is zero by definition (i.e., the displacement of the o n body 
0 o 

reference point with respect to itself is zero). Next we realize that 
£3^ = <D q since all coordinate systems are attached to the o^ body. 

This reduces Equation 94 


F T -m T R + I m k {®0*[t*k> 0 + 7 kc] 
k«o 

+ <S 0 x(w 0 *[tR| t ) 0 + 7 |, c ]>} (102) 

k r I 

+ I m k [(r kc ) k + 2ffi 0 x(7 kc ) k J 

Equation 102 may be rearranged 

“■ a ^ K f I 

F t - m T R + wxI m k [(R k ) 0 +r kc J 
k»o 

K r 1 

+ Oo»««o*f "^[(IVo+'xcJ) (,03) 

+ £ l w h[ ( “kc , k + 2S o* (t kc , k] 
k«l L J 


Since Equation 47 of Reference 6 is the definition 


m T T c 


K 

I 

k a o 


m, 


[«•„>.♦»*] 


(104) 


and since r can be considered approximately zero, then Equation 103 
reduces to 

u K r i 

F t • m T R + £ B » k [(“|, c ) k +2O 0 x(f kc ) k J (105) 
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Neglecting the small coriolis term, as was done in Equation 100 yields 
the final form of Equation 94. 

K 

F t * m T R+j£ m k (r kc ) k U 06 ) 


Note that Equation 106 is the same as the SDF-2 equation (see 
Equation 92 and realize that R is except for the summation 

involving the struts. 

Equations 101 and 106 are therefore the modified equations of motion 

for a vehicle with moving struts. Note that if t.ie struts are not 

accelerating, i.e. (r, ) =0, the modified equations of motion 

kc k 

reduce to the original SDF equations of motion for a rigid body. 

(4) Discussion of Modified Equations of Motion 

The vector F-j. in Equation 106 represents the total applied 
vector force on the vehicle. Because the gears have been added to the 
problem, F^ is composed partly of the ground reaction forces, and 
partly of the forces making up F,, which are gravity, aerodynamics, 

and thrust. The summation ^ m k^kc^k can 156 ^ 00 * cec ^ u P° n as the 

k»l 

inertial force of the accelerated gear struts. 

The vector M in Equation 101 represents the total moment of all the 

O x L. 

vehicle applied forces about the o Ln body reference position (which was 

chosen as the nominal vehicle mass center). M includes the moments 

o 

of the ground reaction, aerodynamic, and thrust forces. The summat , 'n 

X m k [^k^o + *kc ]“Wk can be looked upon as the moment 

k«l 1 

about the o^ body reference position of the gear strut inertia forces. 
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Considerable effort was expended in reducing the general Equations 94 
and 95 to a reference point, which was defined as the "nominal mass 
center." This extra ejfort was not without cause, for it simplifies 
the equations of motion and allows the present moment of inertia definition 
in SDF-2 (which was defined about the vehicle mass center) to be 
applicable. 

Vector Equations 101 and 106 yield six scalar equations for the second 
derivatives of the six rigid-body degrees of freedom for the vehicle. 

Before these equations can be numerically integrated, the ground 
reactions and the strut accelerations must be defined from known quantities. 

c. Auxiliary Equations 
(1) Coordinate Systems 

Figure 21 shows the relationship between the strut coordinate 
system T^, I . , 1,^ and the body fixed coordinate system T xq , T yQ , T zo> 

The coordinate transformation is 

'xk = cos &k I* o -sin0klzo 

Tyk - T y0 U07) 

l Z k = sin^k *xo + cos ^k •ro 


Equation 107 can be conveniently written in matrix form as 


where 


- 


- - 


“ - 

*xk 


a kn 0 a kis 


~*o 

> 

= 

0 1 0 


1 yo 

•*k 


a ksi 0 a k33 


'*0 





- 


akn = cos@k 
Q kis = - sin#k 
Qksi ■ 4lfl $k 
Okas = co * ^k 


( 108 ) 


( 109 ) 
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The landing runway and its coordinate system is shown in Figure 22. 



Figure 22. Runway Coordinate System 

The runway coordinate system T ,1 , T is located down the 1 

J J xr yr zr xg 

axis a distance R . The runway also has an elevation angle, E D , 

9' .. k 

with respect to the 1^ axis (E R allows the calculation to examine landing 

on elevated runways if desired). The runway length is R^. The 

coordinate transformation between the 1 axes system and the 1 axes 

r J g 

system is 


•xr = cosEpl^g-sinER l zg 
•yr = lyg 

T zr = sinER7 xg +cosErT z5 


( 110 ) 


Equation 110 in matrix form is 


r i 




“ 

'xr 


r Gi i 0 Rgi 3 


~*9 

V 


0 1 0 


iyg 

1 *r 


R G» 0 r Gss 


•ig 


(III) 
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where 

r Gi« = cosEr 
r Gi 3 = -sinER 
R G3 , = sinER 

Rq 33= cosE R 


( 112 ) 


The position of the runway origin. 


'gr 


R , can be written 
gr 

R gr Ug 


(113) 


The position of each strut coordinate system (R. ) can be written 

* o 

(Rk)o “ R kx * xo + R ky I yo + R kz * io (H4) 

Figure 23 shows the strut coordinate system and associated strut 
displacements. 


I 


I 


xk 



Axis Position - 

Fully Extended 
Position of Axis 



I S l 


x—i- 


Figure 23. Strut Displacements 
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The definition of each strut displacement is: 

Fj <c = displacement of k th moving strut mass center 


or 


r kc = 1 r Fk ^kc ^k^ 1 3k 


( 115 ) 


where 

= fully extended position of axle 
S k = main strut displacement 
$ kc = distance between strut axle and strut mass center 
? k2c = displacement of secondary piston 


or 


r k2c = ( r Fk " s k2S + s k2 " s k> 1 zk 


( 106 ) 


where 

S k2S 

extension stop 


stop 



displacement between axle and secondary piston 
displacement of secondary piston from its extension 


(2) Strut Fquations of Motion 

The modified equations of motion, Equations 101 and 106, 


require the variable (r, ) . Note that (r, ) 

KC » KC , 

^ Fk ' ^ h 

of the k n strut mass center as viewed from the k u 
(i.e., the vehicle). Again we apply Newtor's Law 
frame). 


is the acceleration 

coordinate system 
(note in an inertial 
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1 .L — 

The total displacement of the k n strut mass center, R^, can be 
expressed as 

R kc = R + (R k > 0 + U kc , k (II 


Newton's Law states that 


where 


£ R k = m k R t 


F k = one of many vector forces applied to the k 1 " strut 

»» X L. 

R. = inertial vector acceleration of k tn strut mass center 
kc 

The first derivative of Equation 117 is: 

JL X. — _ -S- — — 

R kc = R + ou 0 x (Rk>o + ( r kc>k + w o * 1 r kc>k 

where (R,) = 0 because the o th body is rigid. 

K o 


The second derivative of Equation 117 is 

Rkc - R + ui x (Rklo ♦ w o x [ <"o x * R k>oJ + < r kc >k 

+ u> 0 X (r kc )k + 2u/ 0 X (r kc ) k + u 0 X | uj 0 X (r kc >k ] 


Since <u and cu will be very small at landing and (R. ) >> (r, ) 

oo * o kc k 

the last three terms of Equation 120 are neglected which yields 


R kc - R + w x(R k ) 0 + Wo x[o Jo X(R K ) 0 ] + U hC >l 


Note that the term R + ou 0 x (r k ) 0 + w 0 x | u 0 x (R k ) 


is the acceleration the strut coordinate system origin and the term 
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(r ) 

v kc' is the acceleration of the strut mass center as viewed from the 
k 

strut coordinate system (i.e., the term required in the modified 
equations of motion). 

Substituting Equation 121 into Equation 118 yields 

£F k = m k {R +w 0 X ('R k )o + ^o x [^o x (’Rk)o] + (< r kc)k} 022) 


Since the strut is physically constrained to move in only one 
direction, 1 k> we obtain that scalar equation of motion by the dot 
product of Equation 122 and the unit vector 1^. 


5/k * 1 zk 


= m k | R + ou 0 X (R k ] 


+ m k < r kc)k ''zk 


( 123 ) 


The term F k . 1^ represents the forces applied to the strut 
along the unit vector T zk and will be replaced by Note also that 

(r, ) can be obtained from Equation 115 and is 
kc k 

<“kc>k =-s'kljk tl24) 

Equation 123 therefore becomes 
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Solving for the strut displacement acceleration , yields 

Sk - | R + <UqK ( )q + UVq K [ OJ 0 K ( Rk)o ]}’ 1 ZK 
- £ Fkz 


( 126 ) 


ihe only term not defined from known quantities in Equation 126 is 



is composed of several kinds of forces: air compression 


forces, orifice drag forces, friction forces, tire forces (i.e., the 
ground reaction), secondary piston force, end gravity forces. Let us 
examine each in turn. 


(a) Air Compression Force 


Figure 19 shows the C-5A strut configuration. There are 
two air chambers. The pressure. P^, in the upper air chamber is 
transmitted through the oil to the upper surface of the secondary piston 
extension stop and the upper surface of the secondary piston. The 
pressure, P k2 , in the lower air chamber is transmitted to the lower 
surface of the secondary piston and the bottom of the main strut. 


Let A^ be the area of the main piston and A^ 
secondary piston. The air compression force. 


be the area of the 
P Ak> acting directly on 


the main strut is therefore P Ak = P^ (A k - A k2 ) + P k9 A k2 (127) 


The pressures in the upper and lower air chambers depend upon the 
main strut displacement, S k> the secondary piston Displacement S^, 
and initial conditions. When the strut is fully extended and the 
secondary piston is against its extension stop, the initial conditions 
are: 

P Qk = preload pressure of upper chamber 
V pk = preload volume of upper chamber 
P o k 2 = preload pressure of lower chamber 
V , 0 = preload volume o f lower chamber 
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• k 

The relationship between P and V in general is PV = constant 

where k depends on the kind of compression process. Since there is 

little time for heat transfer, one might think that k = — (i.e., 

adiabatic process) would be a good representation. Many references 

indicate, though, that vaporization of the oil in the orifices makes 

the process approach isothermal (i.e., k = 1). Experiments show 

that k = 1.06 is a very good representation. For all practical 

purposes, then, isothermal compression will suffice. The equations 

for and for any combination of and are therefore 


_ p ok v ok 

h ‘ V ok -A k S k + A k| S k ; 


(128) 


p 3 P 0k2 v okt 
k * Vok,“A k 2 S kl 


(129) 


Note that the air compression force always acts along the positive 1 k 
direction. 

(b) Orifice Drag Forces 

Figure 24 shows a typical orifice configuration. The 
force, F , required to move the piston at a velocity V generally takes 
the form F = CV^ 
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where the value of C depends on the fluid arid orifice design (i.e., 
shape, roughness, etc). Depending on the orifice design, C can depend 
on the direction which the fluid is moving through the orifice. Note 

that V is the velocity of the fluid upstream of the orifice. 


The secondary piston (see Figure 19) is the surface which pushes 
the oil through the orifices, and it is this surface which feels the 
.najor resistance to the oil flow. The only resistance to oil flow that 
the main strut feels directly is that associated with oil flow through 
the secondary orifice. This resistance can be formulated C, „ S. „ I S 
where C 

that for positive S 
positive l zk direction. 


k2 


This resistance can be formulated 
is the secondary orifice coefficient of the strut. Note 
the orifice drag on the main strut is in the 


(c) Friction Forces 

There are two obvious friction forces associateo with the 
strut: friction between secondary piston and walls of main strut, and 
friction between main strut and wing gear root. 


The friction force exerted on the main strut by the movement of the 
secondary piston can be expressed as S k2 , where C k21 is the linear 

drag coefficient associated with the kth strut secondary piston movement. 
Note that is is positive, this force is in the positive direction 


of 1 


zk’ 


The friction force exerted by the outer sleeve of the strut due to 

,th 


movement of the main strut, 


F^ k , for the k 


strut can be written 


F fk * Ptk N fk 


( 130 ) 


where 


^sk 


N 


fk 


coefficient of friction at wing gear support 
normal force at wing gear support 
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The normal force, N fk , can be caused by many forces such as fit 
forces, inertia forces, or binding loads caused by ground reactions. 
Many references neglect the friction force because the contribution 
of fit normal forces and inertia normal forces to are small in 
comparison to other forces acting on the strut. The references, 
however, indicate that for struts which are not perpendicular to the 
runway, the friction force resulting from ground reaction binding loads 
at the gear support can be considerable. For this reason, the normal 
force, N^, resulting from ground reactions will be simulated. N^ k 
will be formulated once the ground reactions are determined. Since 
F^ k will oppose the strut movement S^,, Equation 130 becomes 


F fk 



S k 

“ ^ N *k \J^\ 


( 131 ) 


Note that when is positive, 

of V 



is in the positive direction 


(d) G.ound Reaction Forces 

i. L 

The force of the ground pushing up on the k strut along 
the strut will be given the symbol F tk . Its formulation will depend on 
the formulation of the ground reaction. Note that F tk is positive in 
the negative direction of 1^. 


(e) Gravity Forces 


written 


The gravity force. 



, acting on the k 


th 


strut can be 


F kg • m k g 


( 132 ) 
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where g is the gravity vector. SDF-2 already calculates the gravity 
vector in the o*'* 1 body axes system as 

5 ■ 9x T xo + «y T yo + 9z T zo « 


However, it is the component of g along the strut direction of movement, 
l zk that is needed. The gravity force, F^ , along the direction 
1 , can therefore be written (see Equation 108) 


where 


F kgz " m k «zk 


9zk * °kzl + ®z °kzz 


(f) Stop Contact Forces 

Whenever the secondary piston is against either of its 
stops, it can transmit a force to the main strut. This contact force 
is given the symbol F ^ and is assumed to act positive on the main strut 
in the negative T zk direction. 

Having looked at the major forces acting on the strut, we now 
return to the main strut equation of motion, Equation 126, and 
substitute the major forces in the summation Z F |< z ' Since F fcz was 
defined positive along the direction l zk and was defined negative 
along the direction 1 ^ the result is 

Sk ■ | S + u 0 x (Rj, ) 0 +• C 0 x [ui 0 * < R|( ) 0 ] } * T x k 


( P Ak“ F ck2 + c k2 s k2 S k2 + C k2L S k2 (,36) 


Fn + 


m k«zk} 
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Equation 136 is rewritten 

S k « {r + tf 0 x(R k ) c + 53 0 x[w 0 x(R k ) 0 ]} • lzk + °k^^zk <*3 7 > 


\ 



(138) 


i, L 

From Equation 138, the sum of the forces resisting the k strut 
movement may be obtained as 

(139) 

s k 

S Fk * m k °k " F Tk * “ p Ak + F ck2 “ C k2 3 k2| S k2|- c k2LSfc2“ F fk |5^| 


In developing the general equations of motion (Equations 101 and 
106), the inertia effects of the secondary piston were neglected 
(valid assumption) because its mass was very small in comparison to 
that of the main strut. In developing the equation for S^, however, 
we find that the position and velocity and a contact force of the 
secondary piston are needed to determine the major forces acting on the 
main strut. Therefore, we must deve^p an equation for (which 
can be numerically integrated to obtain and S^) in a manner similar 
to what was done for S^. The derivation of follows. 


The inertial vector 

Is the same as Equation 

(r, ) is replaced with 
KC k 


acceleration, R^, 
121 except R kc is repl 

( W k y ield1n 9 


of the secondary piston 
aced with R^g and 


k2c 


8 + w 0 x(R k » 0 


+ w 0 x 


£® 0 x(R 


k*0 


] + (7 k2c , k 


(140) 
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Application of Newton's Law yields 

X F k2 • m k2^k2c 


( 141 ) 


where 

_ th 

F ^2 = applied vector force on secondary piston in k n strut 

X. L. 

= secondary piston mass in k tn strut 


Since the secondary piston is also constrained to move along the 
7^ direction, the equation of motion in that direction becomes 


I* 


k2 ' 'zk 


m k2 *k2c *' 


zk 


( 142 ) 


Equation 142 therefore becomes 




k2z 


m 


k2 


{ [ R ♦ w 0 x ( R k ) 0 + w 0 x [ffl 0 x( R k ) 0 ] ] • T zk 


+ (r k2c ) k 


'zk 

043 ) 


where F k2z is the component of F^ in the positive l zk direction. The 
vector (f* k 2 C ) is obtained from Equation 116 to be 


,F k2c J k 


( s k2 ” $ k > ! zk 


044 ) 


Substituting Equation 144 into Equation 143 and solving for S^ yields 




♦ s 


k 
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Note that the acceleration of the secondary piston displacement from its 
extension stop, S^, is coupled with the acceleration of the main 
strut, S k , which is determined by Equation 137. The only unknowns in 
Equation 145 are the applied forces, which are composed of air 

compression force, orifice drag forces, friction force, gravity force, 
and contact forces at the stops, 

(g) Air Compression Force 

The net air compression force (see Figure 19) acting on 
the secondary piston is simply (P k - P^) A^. Note that it acts in 
the positive T zk direction. 

(h) Orifice Drag Forces 

As was mentioned previously, the secondary piston is the 
surface which pushes the oil through the orifices. If the velocity of 
the secondary piston is zero (i.e-, = 0). the force required of 

the secondary piston to push the oil through the primary orifice is 
simply C k S k |s k | where S k is the velocity of approach of the oil. When 
the secondary piston is moving, however, the equivalent velocity of 
approach is - X S k2 (use continuity of flow). The general 
equation then for the force required of the secondary piston to push 
the oil through the primary orifice is therefore 

*kz| . 

Note that is S k and S k2 are positive, the force on the secondary piston 
is in the positive direction of 1 

The secondary orifice also requires a force of the secondary piston 
to push the oil through this orifice Since the velocity of approach 
is S k2 this force is C k? For positive S k2> this force acts 

in the negative T zk direction. 
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(i) Friction Force 

The friction force between the secondary piston and main 
strut walls is simply S k2 - This force acts in the negative l zk 
direction for positive § k2 - 


(j) Gravity Force 


The gravity force is simply m^ g^ is given by Equation 135. 


(k) Stop Contact Forces 

The stop contact force F^ exists only when the secondary 
piston is against either of its stop. Its value is that which is 
necessary to put the secondary piston in equilibrium as seen by the 
strut (i.e., § k2 = 0). Since its equal and opposite reaction on the 
main strut was assumed positive in the negative l zk direction, the force 
of the stop on the secondary piston must be positive in the positive 
1 k direction. 


The summation of forces on the secondary piston is therefore 

In.2* ■ <P|,- p i l 2' A l t 2 + C l,'VTcf *U'|*»-TS7*||*| 

“ c k2 *k2 | I ~ c k2L ^k2 + m k2 9zk + F ck JL 


(146) 


Substituting Equation 146 into Equation 145 yields 


$k2 * [ 


*^k“ ^k2^ A k2 C k 


<4 kZ t 
(S k ^ s k2 



" c k2*k2|*k2| c k2L ^k2 + F ck2J /m kZ + 9zk 
- [R ♦ « 0 * «Rk>e + fl o*[Oo*^k l o]] * 7 zk 


(147) 
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Whenever the secondary piston is at a stop, 
which is necessary to make S k 2 equal to zero. 


the value of F^ 
Therefore F^ 


is that 
is 


F ck2 = {[ R + tu 0 X (R k )o +a b x[tu 0 x (Rk^o] ]' l zk 
"9zk ~' s *k} m k2 - [ (P k -p k2> A k2 
+ c k ls k~ S k2 ) | S k S k2 | “C k 2 S k2 |s k2 | 

— c k2L^ k2] 


or 


F ck2 “ m k2 { ~ S k2 STOP > 


where \ 2 STOP = the value Equation 147 with F removed. 


One could substitute Equation 137 for S k into Equation 147 and obtain a 
much simpler equation. There are some physical constraints in the 
problem, however, that make the simpler equation for more difficult 
to use. Therefore, Equation 147 will be left in Its present form. 

Note that when all secondary piston terms are deleted, F ^ reduces to 

F ck 2 = " C k \ bk | ' 

(3) Physical Constraints 

Equation 137 and 147 will be numerically Integrated to 
obtain the velocity and displacements of the main strut and secondary 
piston. The range over which and are allowed to vary, however, 
are constrained. The constraints are as follows: 


°S s k 1 S kb 
o < S k2 <S k2 7 


(148) 
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where S kb is the maximum allowable displacement of the k th main strut 
(i.e., strut bottoming) and S^y is the maximum displacement of 
secondary piston (i.e., resting against the secondary piston compression 
stop). As the numerical integration of $ k and S k2 approaches these 
constraints, the values of S k , S k , S^, and must be appropriately 
constrained to ensure that the integrated values of S k and S k2 meet the 
constraints of Equation 148. The constraints on S k2 , S k2 m S k , and S k 
are therefore as follows: 

(a) If | | < Egk »then S^ >0 ond > 0 


(b) If (S kb — E sk ) <S k < (S kb +E sk > then S k < 0 

•* e 

end S k < 0 to exist 

(c) j S k 2 | < Esk2 > S k 2 £0 and *S k 2 2 0 
to exist. 


where 


(d) 


lf (S k2T ~ E sk2> ^ s k2 i < s k2T + E sk2 > S k2 <0 

M 

and S |,2 <0 to exist 


E k = error allowed in integrating S. around the positions zero 
and S kb . k 


E 


sk2 


= error allowed in Integrating S 
and S k2T . 


k2 


around the positions zero 


If the values of S k do not fall into constraints "a" or "b", 
then the numerical integration of S k proceeds normally. If the values 
of S k2 do not fall into constraints "c" or "d", then the numerical 
Integration of ^k2 proceeds normally. In this way the constraints of 
Equations 148 are mot and the "real system" is simulated. 
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(4) Ground Reactions 


The ground reactions interface the landing simulation in 
two main areas: strut applied forces (i.e., N.^ and F +u ); ground 


tk ; 


reaction contributions to M Q in Equation 101; and ground reaction 


contributions to Fy in Equation 106. 
reactions in general. 


We begin by examining the ground 


(a) Ground Reaction Discussion 

The ground reactions, obviously, do not exist until 
any or all wheels touch the runway. Examine Figure 25. 



As the aircraft settles onto the runway, the tires are depressed 

and a certain portion of the tire surface (I.e., the tire footprint) 

Interfaces the runway. The foot print Is allowed to skid on the runway 

with a vector velocity V (note V„is in the runway plane). 

9PTk 9PTk 

The vector Vgp depends on the wheel axle velocity and the rotation 
rate of the tire. 
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The ground reaction normal to the runway, Fy ^ depends on the 
tire deflection, 8^, and can be written 

F TRrk * F TR*k **R (,49i 

where 

^'Rzk = “ n k* k(S|<) 050) 
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The variable is the coefficient of friction between the tire and 
runway and, in general, is not a constant. The predominant dependence 


is en 


Vlk 


and the footprint pressure, P 


of 

fx j. is found to decrease with increasing j 
increasing P pj< . The nominal footprint pressure, P 
as 


In general, 


Pk = 


l F TRzk i 
h k Apk 


pk‘ 

and decrease with 


pk’ 


can be written 


(155) 


where A pk = the footprint area for one tire. 

A j, is mainly a function of the tire geometry and the tire deflection, 
5^. Page 62 of Reference 20 gives the empirical formula for A pk as 

A pk s 2 25(8 k -.03W k )/2r ok W k 056) 


where 

w fc = tire width 
r Qk = tire outer radius 

Note that A pk < 0 for < .03 W^, which has no apparent physical 
meaning. Knowledgeable people at WPAFB (ASNFL) indicated that P pk was 
essentially the airpressure of the tire, which remains approximately 

i 

constant during landing. The whole point of this is that it may not be 
necessary to include the variable effects of P pk on during landing. 
Therefore, will be assumed as a function of | ^GPTk J 


The total ground reaction force, 


can therefore be written 


F TRk - F GPk + F TR/k 


(•57) 
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If F Gpk is written 

F GPk = F TRXlJxR + F TRykT y R U58) 

then substituting Equations 149 and 153 into Equation 157 yields 
Fypk in the runway coordinate system as 

F TRk - F TRXk 1 XR + F TRyk^yR + F TRzk~zR (159) 

Furthermore, if the tire footprint velocity V GpTk is written 

v GPTk = v TXk l XR + v Tyk>yR (IwO) 

then Equation 154 becomes 

F GPk^k F TRzk l^f k | 'XR +^k F TRzk|7^T k| 'yR <^> 

Comparing Equation 158 and 161 yields 

F TRXk : /*k F TRzk i^^ k | (162) 

F TRyk =Mk F TRzk f^; k( (163) 

Note that | V GpT)< , J is 

1 VGPTk l= v GPTk = / v TXk + V 7yk 064) 


Equations 150, 162, and 163, therefore, define the total ground reaction 
Fypk ( see Equation 159) as seen by the runway coordinate system. 
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(b) Formulation of F.^ 

F.^ is the ground reaction component along the negative 
T zk direction in the gear coordinate system. Since the strut is not 
necessarily perpendicular to the runway, the formulation of F.^ in 
terms of Fj RXk , F TR y k , F TRZ|< , is not intuitively obvious. The 
matrices defining the transformation have already been defined in 
Equations 91, 108, and 111. In particular, the transformation from 
the runway coordinate axes 1 R to the strut coordinate axes T R is 

pxkl I a kn 0 ^kisl I ^2 ^3 | I r Gii 0 Rgs,] f UrI 


m i rrio mi 


I 0 


a k3i 0 


a kS3 I n 

J L 


I "2 "3 


R Gl3 0 R G3J *zR 


<yR 065) 


The product of the two matrices to the right make Equation 165 


'xk t*kn 0 a ki3 r L.i r L i 2 R Li3 *xR 


yk = 0 I 0 


R L 2 I r L 22 r L 23 1 yR 


where 


! zk |~ a k3l 0 Qk33 R L3I R E32 R L33j |jzR 


R L 11 : R Gl I + A 3 R Gi3 
R Li 2 ' ^3 

R l_!3 ’ ^1 R G3I + ^3 R GS3 
R|_ 2I : nt|R(j|, + m 3 R Gl3 

R L.22 = m 2 

RLlS^^lRGSl + m 3 R GS3 
R L 31 = n l R Gn + n 3 R G I 3 
R L 32 1 n 2 

R L33 = n l R G3I + n 3 R G33 


! 
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The final expansion of Equation 166 is 

f"'xkl I R luk R I izk R li3k 


yk ] = R l2ik ^122k R 123k 


•zk R i3ik R 132 k R 133k 



= a ki i R l_i i + a ki3 R l_3i 
= Qki i R L 12 + Qkl3 r L 32 
= a kn R Li3 + otki3 R L33 


R 121 k = 


= R (_22 


R l3lk = a k3l R Lll + Q k33 R L3l 
R 132k = a k3l R Ll2 + Qk33 R L32 
R 133k = Qk3l R Ll3 + a k33 R L33 


By definition of F Tk> it can be fcrmulated as follows 

F Tk = F TRk* t ~ 1 zk > 

From matrix 168, 1 ^ can be written 

•zk = R Isik 'XR + R l32k'yR + R iS3k 'zR 


Substituting Equations 159 and 171 into Equation 170 yields the formula 
for F Jk . 


F Tk z ~ < F TRXk R lsik +F TRyk R H2k + F TRzk R lssk> 


( 172 ) 
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Equation 172 is used to determine F Tk , which is needed in Equation 138 
to finally determine S^. 


(c) Formulation of N.^ 


N^ k is the ground reaction component normal to the direction 

If F™, (see 


normal to T . ) at the axle. 


of strut movement (i.e., nur ina i i QL LMC aA,c « 11 1 
Equation 159) is written in the strut coordinate system as 


F TRk s F dxk Ixk + F dyk 'yk + F dzk 1 zk 


then N^ k is 


N fk = / 


F dxk + F ' 2 


dyk 


(173) 


074 ) 


The forces F dyk and F dyk are determined from F JR)(k , F JRYk , F TRZk 
and Matrix 168 to be 

F dxk = F TRXk R tnk + F TRyk R ll2k + F TRzk R li3k * ,75 > 


F dyk = F TRXk R l 2 ik + F TRyk R 1 22 K + F TRzk R l23k < ,76 > 


Equation 174 is then used to define N fk in Equation 130, which defines 
the friction force F^ k used in Equation 138, which is ultimately used 
to help define S k in Equation 137. 

(d) Ground Reaction Contribution to F^ 

The vector F^. in Equation 106 is the total applied vector 
force on the vehicle. The thrust, aerodynamic, and gravity force 
contributions to Fy have already been ionnulated in SDF-2. To the 
SOF-2 force forruylation must be addea the contributions of the ground 
reactions of each gear. Since SDF-2 ultimately needs the forces 
expressed in the body fixed axes system (l xo » ^ y0 * 1 20 )» the total 
ground reaction must finally be written in this system. 
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The total ground reaction vector, F^. R , may be written in the 
runway coordinate system as 

F TR = ^ F TRk= F TRX<XR + F TRy~yR +F TR2^zR (* 77 > 

where 

f TRX = \ F TRX k 

F TRy=^ F TRyk 0 7 8) 

F TRz - ^ F TRzk 
y 

and £ implies a summation over all of the gears on the vehicle. 
Equation 177 is then transformed to the body axes through the matrix 

-| p 

R Ul R Ll2 r Li3 'xR 

R L2l R L22 R L23 *y R (179) 

R L 3t R L32 R L33 'zR 

which has already been defined by Equation 167. If Fy R is written in 
the body axes system as 

P TR = F TRA ! xo + F TRB •yo + F TRC'zo t ,80 > 

then F tra , F trb , and F TRC become 

f TRA‘ r Lii f TRX + R Lt2 F TRy + R Us F TRz 

f TRB = r L2i f TRX +R L 22 F TRy + R L23 F TRz (,8I > 

f TRC= r Lsi f TRX ♦ R l32 F TRy+ R Ls3 F TRz 

Fj R a» f jrb* an<J ^trc are components of the total ground reaction, 

F tr , for all the gears written in the body axes system. 
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(e) Ground Reaction Contribution to M q 

The vector M o in Equation 101 is the total moment about the 
nominal mass center of all the forces applied to the vehicle- The 
thrust and aerodynamic contributions to M q have already been formulated in 
in SDF-2. To the SDF-2 moment formulation must be added the 
contributions of the ground reactions of each gear. The position 
vector, R . , to the point of application (i.e., the footprint) of 

pK _ 

the k tn gear ground reaction force, Fy^, can be written (see Figure 
23 and Equation 114) as 

R Pk = ( R k>o + r k+< r ok'8k>'zR 0821 


where 


th 

ok = outside radius of the tire on the k strut 


The moment vector, RJ< , of the k th gear ground reaction, Fy Rk » is, 
therefore 

MTRk= r PK x F TRk 083) 

In o*-der to evaluate the cross product in Equation 183, R pk and Fy R(< 
must be written in the same coordinate system. Since F TRk is already 
expressed in the runway coordinate system, R^ k will be transformed to 
this system. 


The vector ? k (see Figure 23) can be written 

^k = ( r Fk-S k )T zk 

Transforming r. to the body axes system through matrix 108 yields 


( 184 ) 


r k = i r Fk - S k MQ k j|lxo + Q kss *zol 


( 185 ) 
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Adding Equation 185 to Equation 114 yields 


< R k>o + r k = [ R kx + dk3i f r Fk“ s k)] >xo + [ R ky] 'yo + ( ,88 ) 

[ R kz +Q k33^ r Fk~Sk)]'zo 


Equation 186 is then written in the inertial axes system 1^, 1 , l zg , 

as 

* R k)o + r k = R AXk *xg + R Ayk *yg + R Azk *zg (187) 


where (see matrix 91) 


R AXk 


£ 1 ftij n | 


[ R kx + Qk3i (r Pk“ S k )] 

R AYk 

= 

i. z m z n 2 


M 

R AZk 


^3 m 3 n 3 


[ R kz + Qk33i r Fk _s k>] 


(188) 


Finally, Equation 187 is transformed by matrix 111 to the runway 
coordinate system, T R , and added to (r Qk - S^) l zR to yield R p( , 
expressed in the runway coordinate system as 

R pk = Dxk T kR + D yk 7 yR + D Z k 'ZR * ,89 > 


where 


Dxk 1 r Gii R Axk + r Gis R Azk 

D Yk - R Ayk l ,9 °) 

°Zk 1 r Gsi r Ax'. + R G33 R Azk + r ok" S k 


Now that 
aS ^TRk’ 


(Equation 189) is written In the same coordinate system 
Equation 183 can be written in determinant form as 


^TRk 


>xR >y R 'zR 
^xk ®yk 9 zk 
F TRxk F TRyk F TRzk 


( 191 ) 


no 
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Expansion of Equation 191 yields 

M TRk = MTRxk ( xR + M TRyk ! yR + M TRzk>zR 

where 

M TRxk * °yk FTRzk-°2k F TRyk 

M TRyk * D zk F TRxk“ D xk F TRzk 

M TRzk * D xk F TRyk“°yk F TRxk 

The total moment, of all the gear ground reactions can 

be written 

**TR ■ Z M TRk * M TRx UR + M TRy *yR + M TRz 7 
k 

where 


m TRx 

* I M TRxk 


k 

M TRy 

* Z M TRyk 

M TRz 

* Z M TRzk 


If Equation 194 is written in the body axes system as 



'Tx 


*xo + M Ty *yo + 


M 


Tz 'to 


111 


* 


(192) 


(193) 


therefore 

ZR <‘ 94 > 


(195) 
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then by Equation 167 Kj.y and become 


*Tx~ 


R UI R L!2 R LI3 


M TRx 

M Ty 

m 

R L2I R L22 r L23 


M TRy 

m Ti 


R L31 R L32 R L33 


M TRt 

_ 


097 ) 


r ,, and M.^ are the total ground reaction moments (written in 
body axes) that make a contribution to M q in Equation 101. 


(f) Tire Deflection, 8 k 

In developing Equation 159 for the ground reactions, two 
new variables, 8, and V_ nT . , were introduced. The formulation of the 

.r K br I K 

strut tire deflection 8,., follows: 


The tire deflection arises because of a physical constraint between 
the position of the wheel axle and the runway, i.e., if the height of 
the wheel axle above the runway, Z Qk , is less than the unloaded 
outside radius, r Qk> of the tire, the difference must be the tire 
deflection. 


8 k * r ok + x ok 


098) 


The value of Z Qk is added in Equation 198 since Z ok will be negative 
when the axle is above the runway. Note that Equation 198 assumes that 
the tire shape is a sphere--which is good, provided the wheel axle is 
not far from being parallel to the runway plane (i.e., wings parallel to 
runway). 


Equation 198 is modified by letting the runway have an arbitrary 


profile « 
th 




where « is a function of the position, X 


k strut down the runway, 
compute the effects of 


The « 




Rk’ 


of the 


allows the calculation to 


8 k * r ok ♦ *ok ♦ <u Rk ) 


(199) 


112 






AFFDL-TR-71-155 
Part II 


step inputs (or other runway profiles) on gear dynamics. Note the 
_ th 

following: 8. < 0 means the k tire surface is above the runway, 

K iL 

8. = 0 means the k tn tire surface is just touching the runway; 

8 k > 0 means the k tire is on the runway and is deflected. The. 

formulation of Z Qk and X Rk follows. 

— fK 

Let the vector *Rk represent the displacement of the k n gear axle 
* as seen by the runway coordinate system origin. The vector R Rk can 

therefore be written (see Figure 18, 21, 22, and 23) 

R Rk » ft - RgR r <Rk) 0 + 7 k < 200 > 

where from Equation 90 

R • Xg T„g + Yg Tyg + Zg J Z g (201) 

and from Equation 113 

R gR a R gR *xg 

and from Equation 114 

<R k*o ’ R kx 'xo 4 * R ky^yo + R kt^*o 

and from Equation 184 

F k * (r Fk”V T zk < 204 > 

Now R Rk may also be written in the runway coordinate system as 

^Rk * *Rk Ur + y Rk*yR + z ok *IR (205) 

The formulation of X Rk> Y Rk and Z Qk from known quantities Involves several 
coordinate transformations which follow. 

The vector (R k ) Q + r k has already been formulated in the inertial 
coordinate system 7^, 7 , (see Equation 187) as follows: 

(R k l o +7 k * R A*k **9 4 ’ R Ayk^y9 4R A*k'sg < 20 ®> 


J 


( 202 ) 

(203) 
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Substituting Equations 206, 202 and 201 into Equation 200 yields 

R Rk s (X 9 ' R qR + R Axk ) 

+ (Yg + R Aylt Hyg (207) 

+ ( Z q + RAzk*^9 


> 


/ 


Equation 207 can be transformed to the 1 R axes system (i.e.. Equation 
205) through matrix 111 as follows 


x Rk 


R GM 0 r GI3 


,x a" B 5R + R A«k l 

y Rk 

C 

0 1 0 



*Ok 


&G3I 0 r G33 


IZ a +R A,l.* 


The value of X Rk is used to determine « (X Rk ). The values of Z Qk and 
r , (a constant) then completely determine the tire deflection, 8 .,, 
for the k th strut. 


(g) Tire Ground Plane Velocity, 

The total velocity Ry k of the bottom surface of the k th 
strut tires as seen by the runway can be written 

R Tk * ^GPT h + v T«k T zR (20S) 

Substituting Equation 160 for Vgpy k yields 


ft 


Tk 


v T*k **R ♦ v Tjrk T yR +v T*k ( *R 


(2IO) 
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The formulation of V TXk , V.^ and Vj Z(i , 


from known variables follows: 


_ > l. 

Let the vector R.^ be the position of the k strut bottom tire 
surface as seen by the runway coordinate system origin. R^ k can 
therefore be expressed as 


R n . S-R, R + (H |l ! 0 +r ll + (, 1)k -8 k )T lR (211) 

where ^ r ok' 8 k^ ^ZR is the vector from the wheel axle to the bottom 
surface of the tire. The first derivative of Equation 211 yields 


*Tk ' *-*«R + <«kN> + ^+£Kk“ 8 k>‘*R] <2121 

where from Equation 90 

R * Xg T„g + Yg Tyg + Zg I 2 g ( 213 ) 

and from Equation 113 

RgR “ 0 (214) 

and from Equation 114 

'*»'• . 3„.(S k ) 0 . Co.lR,,,!,,,* R k y ) y0 + R kl ) J0 ) 12)5) 

and from Equation 184 

' '*k T *). + "o* (r Fh-»k ,T I». <2 ,6 > 

and 

Tt [ <r oh~ ^k^Tj^] ■ -&j, Tj R +( w 0 + «Tk U,r ok "^'iR (2,7) 


where<3-^ Is the rotational vector velocity of the tires on the k th 
strut as viewed by the vehicle body axes. Equations 213-217 each must 
be written in the runway coordinate system. We begin with Equation 217. 
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First S k is assumed negligible and » w 0 .Therefore 
Equation 217 becomes 

li [ <r ok“ 8 k , ^zR]“ 5J Tk x(r ok”^k , ^zR ' 2,8) 


Since the wheel is constrained to spin about its axle, 2^ can be 
written 


S Tk * w Tk ! yk 

Since the vector l^ k can be written (see matrix 168) 

*yk 8 R I2lk *xR +R l22k *yR + R I23k*zR 

then Equation 218 can be written in determinant form as 


(219) 


( 220 ) 


£[l'o;,- V'.r] • « 


Tk 


xR 


yR 


ZR 


R I2Ik R I22k R I23k 

0 0 (r ok -8„) 


( 221 ) 


Expansion of Equation 221 yields 

dt [^ofc-^k^zR] " tt, Tk (r ok“ 8 k ) R I22k*xR 
“ w Tk (r ok“ 8 k * R I2lk‘yR 


( 222 ) 


Equation 215 and 216 will be expanded together. First the vector T k 
in Equation 216 is written (see matrix 108) 


1 *k * °kS) ‘*o + 0 k33 ‘*0 


(223) 
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Therefore Equation 216 becomes 



r k * " $ k (a k3l + °k33 ^*0 ) 



•xo ^yo ^zo 

+ (. Fk - Sjt ) P P r 

°k3l 0 a k33 

(224) 

where 

f <o * P T x0 + Q ^yo + f ^zo 

(225) 


On *, panding, Equities 224 becomes 


K * 

i' 

i“*k °k31 + “ a k33 <r Fk “*k 

>] '«,> 


+ 

[ * r Fk “ *k ^ f 0 !;3I ~ P °k33 *] 


(226) 

+ 

[“ *k 0 k33 _q0 k3l ,r Fk”*k ) ] 

| ^zo 



Equation 215 on expanding becomes 


<«*lo * 


*xo *yo ! zo 
P Q r 

R kx R ky R kz 


* <0 R hz- rR ky ,T xo 

♦ ,,R k»-P R kz iT yo 
+ ( P f ~ q r, fcx ^jo 
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The vector summation (R,) + r. can therefore be written 

k o k 

* R kt + ^k * R Dxk ^xo R Dyk ^yo + R Dzk Uo 

where 


R Ox k 

s “*k a k3l +£,a k33 (r rk“ 8 k 1 


+ (qR kz -r R ky ) 

R Dyk 

s 1 r Fk ~ s k H r °k3l “ po k33 l 


+ (rR kR - pR kz ) 

R Dzk 

s “*k°k33 _qo k3l U Fk” s k l 


+ (? R ky - qR ks ) 


Equation 228 is now transformed to the inertial axes 1 


9 


Dy Matrix 



*\ m l "I 


R Dxk 

(R k ) 0 F7k » 

#2 m g n 2 


R Dyk 


^ 3 m 3 n 5 

MW 


R Dzk 


and added with 



to yield 


R - R gR "MFVo+r* 


R DXGk *xg 4 R DYGk 'yg 
+ R DZGk ‘z© 


9 


( 228 ) 


(229) 


91 as 


(230) 


( 231 ) 
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where the inertial axis system components of the velocity of the axle 
as seen by the runway origin are 

R DXGk ' x g + *\ R Oxk + m l R Dyk +fl l R Dzk 

R DYGk = *g + R Oxk + m 2 R Dyk + n 2 R Dzk (232) 

R OZGk " 2 g + £ Z R Dxk + m 3 R Dyk + n 3 R Dzk 

Equation 231 is now transformed to the runway coordinate system through 
Matrix 111 and is added with Equation 222 to yield Ry k (see Equation 212) 
expressed in the runway coordinate system. The result is 

«Tk ■ S-S GR + (R k l 0 +r k + £[l, 0k -8 k >T lR ]. 

] R GJi R DXGk * R GI3 R DZGk + w Tk tr ok“ 8 k ,R I22k bxR 
L (233) 

+ [ R OYGk~ “'Tk {,r ek~ ^k } R 12lk]^yR 
+ [ R G3I R DXGk R G33 R DZGk] 'zR 

Comparing Equations 233 nnd 210 yields 

v Txk * R GII R DXGk + r GI 3 *DZGk + ^Yk {r ok ~ 8 k 1 R I22 k 
v Tyk ' R DYGk _ w Tk ( r ok " R I2lk <234) 

v Tzk “ r G 3I R OXGk + R G33 R 0Z 3k 
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Equations 234 for V Txk and Vyy k are then used to calculate the ground 
plane forces Fy RXk and F TRYk ^ see Ec l uat1ons ^62 and 163). Note that 

Equation 234 for is not used but is useful when one realizes it is 

the vertical velocity (i.e., sink rate) of the k^ gear axle as seen 

by the runway. 


(h) Wheel Equations of Motion 

In deriving an expression for the ground plane velocity, 
VfipTk* of the k th tire footprint, the rotational velocity, u;^, was 
introduced. since u»y k changes during landing, another equation of 
motion is needed. Applying Newton's Law about the wheel axle yields 

2>k • n k l k*Tk <235) 

where 

- applied moment about k th axle 

I. _ fmoment of inertia of a tire, wheel, and anything else 
1 constrained to rotate~with that wheel about the axle. 


‘.he applied moments are predominantly the moments of the ground 


Let be the 


reaction force* and braking moments (if applied), 
moment of the ground reaction forces about the axle and let M pk be a 
braking moment which will be determined by the Brake Autopilot as shown 
in Appendix 111. Equation 235, therefore, becomes 


M 


Ak 


M 


Bk 


^Tk 

l"Tki 


n k x k w Tk 


(236) 


Note that the oraking moment always opposes the motion, Solving 

Fqu.3t.ion 235 fordi^ yields 


1 1 M Ak “ M 0k w 


Tk 


Tk 


) r I fc n 


k "k 


(237) 
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Equation 237 is then solved numerically to obtaincu.^. The variable 
must be defined however, in terms of known quantities before tne 
integration can proceed. 


The vector from the axle to the footprint of the tire is (r ^ - 8^) 
T zr . The vector moment of the ground reaction Fy Rk (see Equation 
159) about the axle is 

M Ak = < r ok ~8 k ) 'ZR * F TR k (238) 


Substituting Equation 159 for F TRk 
form yields 


and writing Equation 238 in determinant 


M Ak - 


*XR 

0 

F TRxk 


! YR *zr. 

0 i r ok~^k) 

F TRyk F TRzk 


(239) 


Expanding Equation 239 yields 

MAk = - F TRyk< r ok _ Sk)T K R + FTRxk( r ok “8k> 'yR 


(240) 


Since (see Equation 236) is the component of along the axle 
direction T ^ Matrix 168 performs the desired transformation of 
Equation 240 and shows to be 


M Ak r “ t TRyk ( f ok ~8k> R lzik + F TRxk< r ok _ 8k ) R| 22 k 


(241) 


This completes the derivation of all quantities needed for all the 
equations of motion (see Equation 101, 106, 137, 147 and 237). 
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d. Expanded Equations of Motion 

Five kinds of equations of motion were needed to describe the 
landing vehicle with gears. Each equation is repeated here for 
convenience. 

M 0- i 0 *^o + £ *>o x < l o ,w o> + £ m k[< R k>o + r kc] X(rkc>k OOll 

k=l L J 


F t = m T R + 2, m k U kc 
k*l 


(106) 


S k = { R + cu 0 X (R k ) 0 + oj 0 X [ a» 0 X (R k ) 0 ] } • I zk 

+ a k " 


(137) 


Ski =[<Pk-Pk 2 >A k2 + C k (S k - •&! S k2 ){ 5 k - S k2 | 

“ ^k?Sk 2 [S k2 ] - C k2 i_ S K2 + ^ckz]/ m k2 + dzk (147) 

- [ R !- X ( R k ) 0 + X f w 0 X (R k ) 0 J] -T zk + S k 

^Tk = < M Ak " Wgk ]/ 1 K n K (237) 


Equations 101, 106, 137, and 147 contain vector operations (though 
completely defined) that are yet to be expended before numerical 
integration can proceed. We begin with Equation 101. 

Since the counterpart to Equation 101 in SDF-2 is written In body 
fixed axes, the vector operations in Equation 101 must finally be 
expressed in body fixed axes. The vector (R k ) 0 (see Equation 114) 

is first transformed to the T k axes system through Matrix 108 and then 
added ^ r Fk ' S kc^zk ( see Equation 115) to yield 

(R k ) 0 +(rp k “Skc^zk sR RlCGXUk + R RiCGX ’ yk + R RlCGzIrk ( 2 **) 
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where 


m RICGX = ^kii R kx + a ki3 R kz 
r RICGV = R ky 

r RICGZ‘ <2kai R kx + a k33 R kz + r Fk -s kc 


(243) 


Note that Equations 243 are constants. The vector (R^) Q + r kc can 
therefore be written 

( R k)o + r kc = R RICGX^xk +R RICGY^yk'H R «ICGZ - s k>'zk (244) 
Since (r. ), can be written (see Equation TIB), 


( r kc 1 k = “ ^k ! zk 


(245) 


Then the vector cross product in Equation 101 can be expressed as 


'XK 'yx 'ZK 

[< R k>o + r kc] x$r kc>k = R RICGX r RICGY (R RICGZ _s k) 


-SxRRiCGY'xk + S k RRiCGX'yk 


(246) 


Transforming Equation 246 back to body axes (see Matrix 108) yields 


( R k )q + r kc] x » r kc>k - ~Ok,, S h RricgY 'xo + 


(247) 


"Sk r RICGX Iyo -axis S^RricGY' zo 
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Substituting Equation 247 into Equation 101 yields 
M 0 - lo ■ ujq + u»o ^ (to * *^0 ^ 

K _ _ (248) 

£m k S k (“Q k |, R RICGY 1 xo + ^RICGX^yo — ct ki3 ^RICGY *zo 1 
k= I 

Equation 248 (i.e., also Equation 101) is now expressed in body axes. 

Equation 106 is written in body axes by substituting Equation 245 

f °r (r kc ) k and realizing that l zk = T y0 + o k33 l z0 (see Matrix 

108). The result is 

.. K .. 

Fy = m y TF — £ m k S k ( a k3( I xo + a k 33 *zo 1 (249) 

k= I 


Equations 137 and 147 contain the vector expression 
{r +i 0 x (R k ) 0 +J 0 x[w 0 x<R k ) 0 ]}-~zk 
Each expression is taken in turn. 

0) I. i !k 

The vector R is obtained directly from SDF-2 as 

R - Fx Ixo +■ Fy Iyo “F Fz Ijq 
my my my 

Since by Matrix 108 1 2k - a fc3 , T x0 * a k33 

u - 

then the scalar R . l zk is 

R l lk : ^<>kJI + *f T a kS» 

(2) [u X(R k ) 0 ] -T lk 


(250) 


(251) 


.4 
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The vector^ is obtained directly from SDF-2 as 

<*> 0 - P ' no 0 I yo r "^ 2 o (252) 

The vector (R^) is Equation 114. The vector cross product & o x (R^ 
is therefore 

*xo 'yo 'zo 

w o p q r = (qRkz ~ ' R ky> 'xo + 

R kx R ky R kz (255) 

R kx “P R kz ^~yo +( P R ky~q R kx)^zo 


Using the expression previously used for 1 . yields the scalar 
expression 

[ w 0 x (R k ) 0 ]-Tzk to be 

^ a» 0 X (R^ ) 0 j -1 = Q k3i id R k —rR ky) (254) 

+ a k 35 !pR ky -qR kx ) 

(3) cu 0 X[a»o* ( R k *o] *' zk 
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The vector w Q x u Q x (R k ) 0 J is 


co 0 x [a) 0 X(R k ^J= 


xo 


yo 


zo 


p q r 
(q R ta~ rR kyK r R kx ~ pR kzMP R ky -qR|tx)| 


[q(pRky-qRkx ) ~ r(rR kx~P R kz)j Uo 
+ [ r (qRk z - r Rky)-P(PRky-q R kx)] 'yo 


(255) 


+ [p( rR kx - P R kz>-q(q R kz" rR ky>] 



Once again using the expression for yields the scalar u Q 
<\>ob T zk t0 be 

Uq X [* 0 X (R k ) 0 j T zk = a k9 i [q(P R ky -< l R kx) 
-r ( rR kx'P R kz )] + a k95 [p( rR kx“P R kz) lq R kz~ rR ky)] 



x 


(256) 


Summing the results of paragraphs (1), (2) and (3) and letting the result 
be S Rk which is the l zk component of the inertial acceleration of the 
strut axes system, yields 

SRk i { R+w o*( R k)o +w o*[ w o* : ' R k)o]}' , zk :: 


a ksi + a kS 3 ■*" a ksi fq R kz - ^ R ky) 

+ akis <P R ky“P R kx> + °k3i [p<P R ky-P R kx> 

- r(rRk x -pRk z )] + Qkjs [ p(r R kx -pRkx>- q (Q R kz _r R ky)] 


(257) 
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Equations 137 and 147 therefore become 


s k = s Rk + Qk -9zk 


( 258 ) 


Sk2 r [< p k -p k2> A k2 + c k<S k - S k2 ) S k - S k2 


^k2 §k2 


“ c k2L s k2 + F ckt] /m k2 


(259) 


+ 9*k ~ s Rk '''Sk 


Equations 237, 248, 2 , 258, and 259 are the five equations of motion 
written in the proper coordinate systems and are in a form convenient 
for numerical integration. 
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SECTION III 
DISCUSSION 


The SDF-2 equations of motion are written with the total applied 
force Fy and total applied moment M q written in body axes as 

* F x Uo + Fy ^yo + ^z^zo (260) 


Mq ■ LIjo + Mlyj + NIjj 


(260 


By rewriting the modified equations of motion (i.e., Equations 248 and 

249) with the added terms as applied forces and moments (i.e., on the 

left sides of the equations), the added terms can be looked upon as 

changes in the SDF-2 applied forces and moments brought on by adding the 

landing gears and ground reactions to the simulation. If the new Fy 

and M are written as 
o 


Ft 

M 0 


+ <F y +F yfn ) 

*yo 

+ <F I + F„,) 

^zo 


<L + L m )T „ 0 + (M + M m H y0 
+ IN + N ffi ll l0 


(262) 
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where F xm’ F ym’ F zm’ L m, M m’ and N m are the changes in F x’ F y’ F z’ 

L, M, N, respectively, which are presently in the SDF-2 formulation, 

then the changes are 

F xm * F TRA + Z m k *k °k3l 
F ym * f TRB 

F zm ' F TRC + Z m k *k 0 k33 

* (263) 

Lm * M Tx + I •% * k o k|| R k 
k 

Mm * M Ty ~Z m k*k R RICGx 
k 

Nm • M Tz + I "»n *|, °klS R ky 


This concludes the formulation of the equations defining the SDF-2 
modification to include landing gear dynamics and ground reactions. 
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SECTION IV 

LANDING GEAR FLOW CHARTS 

The following flow charts show the previous equations in their 
required order of calculation. The corresponding equation numbers 
are given on the right s.de of the equation to serve as a handy reference 
in finding a desired derivation. 


130 



AFFDL-TR-71-155 
Part II 


CALCULATIONS DONE ONCE 

°kn = 

cos 6 k 

a ki3 = 

-sir 9 k 

°k3i r 

stnd k 

a k33 = 

cos 6 k 

r Gii ' 

cosEr 

R Gl3 = 

-sinE R 

R G3I : 

smER 

R G33 “ 

cosEr 

R RkcGX ' 

a ki i R kx + °ki3 R kz 



[rl] MATRIX ELEMENTS 

r Lii = ^II R G|| + ^3 R G13 
R Ll2 : ^2 

R Ll3 = k I R G3I + ^ 3 R Gs3 
R L2! = m l R Gll +m 3 R Gl3 
R L22 =m 2 


R L23 :m l R G3l 4 m 3 R G33 
R L 31 ■ n I R G 1 1 + n 3 R Gl3 
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_ £ _ 

[Rl] MATRIX ELEMENTS 

R lnk = J kn r Lii + a ki3 R L3i 
R l 12*\ = °kl I R Ll2 +a kl3 R L32 
R l 13k = Q kll R Ll3 +Q kl3 R L33 
R l2lk = R L2I 
R 122 k = R L22 
R 123k - R L 23 

R l3lk = a k3l R LII +Q k33 R L3l 
R l 32K = a k3l R Lt2 +0 k33 R L32 
R 133k ' a k3l R Lt3 + a k33 R L33 


♦ _ 

R kx + a k3i ( r Fk S|<) 


R AXk 




m, 


R Avi 

R AZk 



R ky 

R kz + a k33' r Fk“ ^k ) 
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R DXK r ' s k a k3i + ( f Fk - s k) £ J°k33 + 9 R kz ~ r R 
R OVk = ' r Fk ~S K ) l ra k3t -po k33 ) +rR kx~P R kz 
R DZk r “ s k a k33 _ < r Fk -S k ) qo k3l +pR ky -qR k » 







































Ska = [< p k - p ka! A kz + C k (S k - ^ S k2 i| S k - ^ S k2 1 4 
■^(12 S k 2 I S k j | - C k j^S k2 j /rn k2 4- g 2k - Sp k + S k 
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_ 1 _ 

ADD F «M. F yM- r «ll* L M. M M- N M T0 S0F - 2 
CALCULATIONS OF F x , F y , F x , L, M, N AS FOLLOWS 

F x * T x-o + > ^F, + F xM 

F y « Ty >y + AF y + F y || 

Fj * Tjj *■ np + i#9j+ AFj + F xM 

L « L t + AL t + / + L m 
M • Mj+ Am y + (n + M|| 

N > N T +AN T + n + N M 


( CONTINUE SDF- 2 A 
V CALCULATIONS J 
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APPENDIX III 

CONTROL MANAGEMENT EQUATIONS 
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SECTION I 

AUTOPILOT PROBLEM DEFINITION 
1. PROBLEM DISCUSSION 

The aerodynamics subroutine in SDF-2 (called SACS - see pgs 41-45 
of Appendi < I) can simulate the aerodynamic effects of a single control 
surface in each of the axes pitch, yaw and roll. The input aerodynamic 
data can also be "staged" (i.e. changed during the running of the 
program) to simulate the deployment of flaps, slats, spoilers, etc. 

The thrust Subroutine in aDF-2 (called TFFS see pgs 46-47 of Appendix I) 
can simulate the thrust and fuel flow of a single airbreathing engine as 
a function of altitude, Mach number, angle of attack, and throttle 
setting. Therefore, SDF-2 as originally formulated has sufficient 
control variables, with limited modification, to control the aircraft 
in six degrees of freedom. The question which remains unanswered, 
however, is that concerning the magnitude of the control variables at 
any time. This is the function of the autopilot which can be stated iri 
questior form as follows: Given the state of the aircraft, what 
values should the control variables be? This Appendix concerns itself 
with the answer to *his question. Figure 26 shows this basic interface 
between SDF-2 and the autopilot. 

In finally arriving at the control variable values, the autopilot 
must necessarily answer the following questions: 

a. What is the desired state of the aircraft? 

b. Given the present state of the aircraft, does an error in 
state exist and if so what maneuver will be done to correct the state 
error? 

c. What control variable values (i.e., elevator, rudder, and 
aileron deflections and throttle settings) will result in the desired 
maneuvers? 


144 


AFFDL-TR-71-155 
Part II 



Figure 2S. Autopilot - SDF-2 Interface 


2. DESIRED CAPALILITIES 

The following is a list of the major capabilities required of the 
autopilot: 

a. Control the aircraft in six degrees of freedom during the 
glide slope, flare, landing roll, and takeoff roll within ‘he capability 
•of a given aircraft. 

b. Perform aircraft control with changing winds. 

c. Be capable of analyzing multiple engine aircraft with 
engine failures and engine reverse. 

d. Simulate <ontrol as aircraft transitions into ground effect. 

e. Examine first order effects of control surface and engine 

lags. 
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f. Examine braking capability and selected braking failures. 

g. Include a drag chute simulation. 

h. Be able to start at any point in the landing or takeoff sequence 
and terminate calculation at selected points. 


3. ASSUMPTIONS 

The basic assumption is that initially the aircraft is on or near 
the desired spatial position and is trimmed in all three axes and power 
such that if no further parturbations occur, (such as wind changes, 
engine failures, changing aerodynamics due to ground effect, etc.) the 
aircraft will remain on the nominal glide slope position and inertial 
velocity. It is further assumed that the approach speeds are lew 
enough (such as 1.2 - 1.3 ^ tall) S0 t ^ iat the deroc *y nam '' c coeff cients 
are predominantly a function of aircraft control surface configuration 
and altitude (i.e., ground effect) and not Mach number. It is also 
assumed that the aircraft is low enough so that runway level thrust 
properties are valid. This eliminates engine performance dependence 
on altitude changes. Variation of aircraft weight during landing and 
takeoff is also considered negligible. 
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SECTION II 
MANEUVER LOGIC 

The maneuver logic concerns itself with the answer to questions 1 
and 2 in Section I. As shall be seen, the maneuver logic determines 
the desired values of angle of attack, angle of sideslip, roll angle, 
and thrust (along with other commands) which define a maneuver to 
correct a state error. The organization of this logic is divided into 
four areas: glide slope, flare, landing roll, and takeoff roll. 

1. GLIDE SLOPE 

The basic requirements of the aircraft while in the glide slope 
phase are two: to be vectorially near the glide slope within an allowed 
error; and to maintain the inercial speed down the glide slope a constant 
We begin by defining the position error. 

a. Postion Error 

Examine Figure 27, part of which is extracted from Figure 22 



Figure 27. Glide Slope Geometry 
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The nominal glide slope position is defined in the vertical plane 

(i.e., plane containing the unit vector T and T ) and has the glide 

zg xg 

slope angle, measured up from the horizontal plane. The glide 

slope origin begins at a distance, h CQ , above (not indicated in 
Figure 27) the runway origin. The variable, h^, is the fixed 
distance between the aircraft mass center and a line parallel to the 
longitudinal body axis which is tangent to the bottom tire surface of 
the main landing gear. As such, if the aircraft were on the glide 
slope at the runway origin, the tires would just touch the runway. 
Figure 28 shows this geometry for the vertical plane. The allowed 

qlide slope position error in the vertical plane, hi (see Figure 28), 

is defined by the fixed angular perturbation, , and the position 

V 


vector R n . The allowed glide slope position error in the horizontal 


plane, h 


Da 


(see Figure 29) is defined by a similar fixed angular 
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Figure 29. Horizontal Plane Glide Slope 


perturbation, Sq- , and the position vector R^. Note that as the 
aircraft gets closer to the runway origin, the allowed position errors 



and 


and hi 

pa 

and h 
' pa 


reduce gradually to a very small value. The equations for 
will be approximated (i.e., realizing that is small 



as follows: 


eo (Rr| 


(264) 

po * |«r! 

So- 

(265) 


where 
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and 8* and Srr are in radian units. The variables X and Z are 
e 9 9 

available from the existing SDF-2 formulation and the distance R_ n is 


available from the landing gear modification (see pg 81). 
glide slope altitude. 


gR 

The nominal 


hg S , (see Figure 28) is dependent on the position 


coordinate X as follows 
9 

h GS 


Xg- R gR 


tan « GS + h CQ 


(267) 


The vertical glide slope position error, h , (see Figure 28) is 
therefore 

h e * h-h G$ (268) 


Note that h is positive when the aircraft is above the nominal glide 
slope and is negative when the aircraft ib b^low the nominal glide slope. 


The nominal glide slope position in the horizontal plane (see 
Figure 29) is identically zero (i.e., the aircraft ideally should 
remain in the vertical glide slope plane). The horizontal glide slope 
position error is therefore the position coordinate, Y , which is also 
defined in the existing SDF-2 formulation. Note that Y^ is positive 
to the r'ght of the vertical plane and negative to the left. 


In controlling the aircraft to remain within the allowed spatial 
error (see Figure 27 for a pictorial view), it is helpful to know the 
rate at which the glide slope position error is changing (i.e., rate 
feedback). The rate of the vertical glide slope position error, h g , 
can be expressed as follows: 

h e * - Zg + Xg Ton « G $ (269) 

Note that h is coupled with the horizontal velocity X and h is zero 
c y * 

only when the aircraft sink rate, Z , is equal to the apparent rate 
at which the glide slope is falling. The rate of the horizontal glide 
slope position error is simply Y^. The variables X g , Y^, and , are 
all available from the existing SDF-? formulation. 
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Tha total vertical glide slope position error, h y, and the total 
horizontal glide slope position error, hpy, are written as the sum of 
a position error and a rate error as follows: 

h eT s h e + RF h (270) 


hp T = Y g + RF y Yg 


(271) 


The RF^ and RF are input constants and determine the amount of rate 
feedback in the vertical and horizontal planes respectively. This 
completes the formulation of the glide slope position error. 

b. Velocity Error 

If the aircraft is controlled such that h and hpy fall 

within the allowed errors h and h , respectively, this only 

ea pa J 

guarantees that the aircraft is vectorially near the nominal glide 
slope position and does not control the inertial velocity with which the 
aircraft comes down the glide slope. The inertial velocity error, V , 
must therefore be sensed and is formulated as follows: 

V e « V 9 -V d (272) 

Vg is the inertial velocity magnitude (as formulated in the existing 
SDF-2) and V d is the desired inertial velocity down the glide slope 
which is input as a constant. 

c. Glide Slope Dynamic 

The two basic requirements for the glide slope (see Pg 147) 
define a steady-state descent maneuver at constant ground speed. The 
nominal values of angle o. attack, angle of sideslip, roll angle, and 
thrust to perform this maneuver ..an be obtained by a steady-state 
summation of forces acting on the aircraft. We begin by examining 
forces in the aircraft pitch plane. Examine Figure 30 and the 
following definitions. 
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Figure 30. Nominal Forces in Glide Slope 


desired inertial vector velocity down the glide slope 
wind vector 

desired airspeed vector 

pitch angle of aircraft from horizon 

desired angle of attack 

desired flight path angle of V . 

0 

desired thrust vector 
lift vector 
drag vector 
weight vector 
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p = atmospheric density 

S = aircraft reference area 


Cpp = required lift coefficient 


Summing forces in the vertical direction (assume aircraft wings 
approximately level) yields: 

L cos () + T H sin + Dsin(-yo)-W«0 


V 


(273) 


This may be solved for the required trim value of lift coefficient 
for the desired static equilibrium gliding condition as: 


•LR 


[ M gREF+ D R sin *yR* -T d sin <0R>] 


Or S COS( y R ) 


Summing forces in the horizontal direction yields: 

cos 0 r + Wsin ( ~Xr l — Dcos ( —y^) • 0 

solving for T, yields: 

Wsin y^ + Dcos y^ 


cos 


e c 


(274) 


(275) 


(276) 


Equation 274 can also be written as 


W 


% 


LR 


«REF 


?P V od S 


0 R S 


(277) 


wiiere and y R are both set to zero as a first approximation. The 
value of C^p in Equation 277 is used to find an estimate of trim 
and 8 , as will be shown subsequently. T, is then calculated 

using Equation 276, followed by a more accurate estimate of C LR through 
Equation 274, and of Tj through equation 276. 
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C^ R is dependent on | ^ | , which is also dependent on the wind 

vector V through the following relationship (see Figure 30): 


V fld = -^/vTcoT^g-Xqw^ + Yg^+iVjjSineGs-Zgwl* (278) 


The variables X , Y and Z are components of the wind vector, 

S w 3 W gw 

V on the T ,1 ,1 , axes, respectively, and are available from the 

w xg’ yg zg’ 

existing SDF-2 formulation. 


The coefficient of lift, C^, is primarily a function of angle 
of attack, a, and secondarily a function of elevator deflection, Sq, 
through the following relationship: 

C L ■ C Lc + C La a + C La *|a|ji+C L g g Sg + C LSg *| 8g| Sg (279) 

where C Lq , C L(J , C Lq2 , and C L g q2 are the predominant aerodynamic- 
coefficients contributing to lift. Therefore, for a given elevator 
deflection. Equation 279 ultimately defines the desired angle of attack 
a d , to meet the requirements of Equation 274 for the glide slope. The 
nominal elevator deflection for use in Equation 279 is obtained using 
Fquation 302 which will be discussed under Pitch Autopilot. 


Equation 276 defines the required thrust magnitude, T^, to maintain 
the glide slope as a function of W, y R , U, :id# R . The drag force, D, 

is evaluated at and 

0 1 c ojp\ v 0d|* s ,280) 

where 

Cq * c 0o + c Da a rf + c D a * a d* (28)) 

and C Dq , C Qa , and C Qa 2 are the predominant aerodynamic coefficients 
contributing to drag. 


ad 


through the following equations: 
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The equations defining y p and can be obtained by referring to 

Figure 30. 



ton 


V d sin < GS - Z gw 
v d cos e GS“ *9# 


( 282 ) 


9 


R ‘ 



+ Q d 


( 283 ) 


The a d resulting from Equation* 277 , 302, and 279 and the T d 
resulting from evaluating Equation 276, as indicated, are the nominal 
airplane requirements in the pitch plane to satisfy the requirements of 
the glide slope in the vertical plane. Summing forces in the horizontal 
plane yields the basic requirements for the aircraft to remain in the 
glide slope vertical plane and not drift horizontally. The nominal 
requirements are: wings level (i.e., no aircraft roll), and sideslip 
angle of attack zero (i.e., crab aircraft into wind). 

The preceding analysis allows the autopilot to determine the nominal 
trim requirements for a steady ‘descent at constant ground speed for an 
arbitrary wind vector. Unfortunately, changing winds, engine failures, 
changing aerodynamic coefficients due to ground effect, etc., require 
the aircraft to transition between nominal trim requirements which can 
put all six degrees of freedom into andesired oscillations. Of primary 
concern in the glide slope are the long-period oscillations of the mass 
center position about the nominal glide slope position. Though these 
long-period oscillations present no real control problem to the pilot 
in real life, they must be sensed in the autopilot simulation and 
appropriately controlled. A means of sensing these oscillations has 
already been provided through the variables h ^ and h pT in Equations 270 
and 271. 
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The logic for the modification of comes from an understanding 
of the nature of the longitudinal long-period mode, but can be briefly 
stated as follows: when the aircraft is rising above the allowed vertical 
glide slope error, lower the a^; when the aircraft is falling below 
the allowed vertical glide slope error, increase the a^. The magnitude 
of the modification to is determined by the input constant PG S and 

will depend on the particular aircraft simulated. 
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Briefly stated, the horizontal logic is as follows: When the 
aircraft drift to the right of the glide slope vertical plane exceeds 
the allowed value h pa> roll the aircraft - <£ c (left); when the aircraft 
drift to the left of the glide slope vertical plane exceeds the allowed 
value -hp a> roll the aircraft <£ c (right). Since the <£ c command 
will cause a component of the large lift vector to project on the 
horizontal plane, can be small (3°, etc.) and still control the 

horizontal oscillations. 


This concludes the equations defining the maneuver logic for the 
glide slope phase of the problem. Figure 31 presents a summary of the 
equations and logic. 
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2. FLARE 

The basic requirement of the flare maneuver is to transition the 
aircraft from an arbitrary vector position and vector velocity in the 
glide slope, to a desired touchdown vector position and touchdown 
vector velocity, simultaneously, on the runway. For a realistic 
landing, the touchdown vector position is limited K y runway length and 
landing distance required to stop the aircraft. The touchdown vector 
velocity is also limited by acceptable landing speeds and aircraft sink 
rates. Since there always exists the possible effects of engine 
failures, wind changes, and limited aircraft flare capability for a given 
situation, some basic logic is necessary to sacrifice the desirea 
touchdown constraints logically and still make an acceptable landing. 

The problem is typically a guidance problem subject to constraints. 

The development of basic guidance laws to perform this function follows: 

a. Guidance Laws 

The flare guidance laws will be based on a constant acceleration 
maneuver. Flare initiation will be staged on a particular altitude, 
hp, above the ground (see Figure 32). At the time of flare initiation, 
it is assumed that the glide slope cross range control has the aircraft 
sufficiently close to the vertical plane so that the wings may be leveled 
(note aircraft will still remain crabbed into the wind). As such, the 
flare maneuver is essentially a planar problem and occurs near the 
vertical plane. 

Since the runway can have an elevation angle, E^, the guidance 
laws will be developed in the runway coordinate system X^, h^ where 

Xp = aircraft scalar distance down runway 
■ aircraft altitude above runway 
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Figure 32. Flare Coordinate System 


The desired touchdown conditions are as follows: 

x RF = distance down runway for touchdown 
h R p -- mass center altitude normal to runway for touchdown 
x R p = landing speed parallel to runway for touchdown 
h pf - = altitude rate normal to runway for touchdown 


Consider the following general derivation for any scalar variable, 

S, with constant acceleration, S: 

S f :St s +Sf (2B4) 

where 


t = elapsed time 
S. = initial S 
S f = final S 
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Integrating Equation 284 again yields 


S f = S 


.i + s t • 

2 s 


( 285 ) 


where 

= initial S 
Sj. = final S 


Solving for t in Equation 284 yields: 



(286) 


Substituting Equation 286 into Equation 285 and expanding yields: 


Sf-S|= -f-(-^) + s,(^) = 

- Sf-sf _ (S t -Sj)(Sf+Si ) 

2 S ' 2 S 


(287) 


Solving Equation 287 for S yields 


S 


'Sf-Sj) 


(SfSj) 
2(S f -S i ) 


(288) 


Substituting Equation 288 into Equation 286 yields 


t 


s 


ztgiiSii 

§ f + §i 


(289) 


Equation 288 expresses the constant acceleration, S, required to 
transition from the initial state S.. , S^ , to the final state S^, S^., 

In terms of the initial present conditions and desired final conditions. 
Equation 289 expresses the time required, based on constant accelerations, 
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to make the transition in state. Comparable equations for the scalar 
variables and h R are therefore: 


a xr = 


< x rf~ x r' 


(Xrf+Xr) 

2(Xr F -Xr) 


(290) 


2(Xr F -Xr) 
* !^rf +x r) 


(291) 


. _ : . ( h RF+ h R> 

AhR MhRF ~ hR feFw 


(292) 


2(hRF~bR) 
( hRF + ?>R ) 


(293) 


where 



required acceleration of x^ 

time to perform the x^ state change 

required acceleration of h^ 

time to perform the h^ state change 


If Equations 290 and 292 can be satisfied at all times during the flare, 
they guarantee that the x^p, x^p touchdown constraints will occur 
simultaneously and that the h^p, h^r. touchdown constraints will occur 
simultaneously. However, this does not guarantee that the two 
Individual sets of scalar constraints will occur simultaneously. This 
latter constraint can be met by requiring t to equal t^. This is 
achieved by the following logic. 
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Calculate t and t^ from present conditions (i.e., x R , x R , h R , h R ) 
and the desired touchdown condition (i.e., x RR , x RR , h RR , h RR ). 

Compare the values of t and t h> If t h is greater than t x> Equations 
291 and 293 suggest four possible sacrifices of the desired touchdown 
constraints to make t^ and t x equal: 


( 1 ) 

( 2 ) 

(3) 

(4) 


Increase x RR (i.e., land further down the runway), 
Reduce the landing speed, x RR . 


Increase h oir (note h D > h Rp ). 


'RF 


Increase h RR in the negative direction. 


Sacrifice 3 is eliminated since this would require the aircraft to flare 
out above the runway. Sacrifice 2 is eliminated since the desired 
landing speed will already be cl r se to the stall speed. Sacrifice 1 is 
possible provided there is suf' '’ent runway length left after touchdown 
to stop the aircraft. Sacrifice means an increase in sink rate at 
touchdown; this is undesirable, but can be tried once the possibilities 
of Sacrifice 1 have been exhausted. This logic is formulated as follows 
Let 

X RF 3 X TD + D m (294) 


where 

x Tn * distance down the runway before which the aircraft must 
u not touchdown. (This is normally zero but can have a 

positive value which in effect places the closest allowable 
touchdown point ahead of the effective runway beginning.) 

D m = distance down runway from Xyp to desired touchdown point. 

Apply Sacrifice 1 by equating t and t^ in Equations 291 

and 293, substitute Equation 294 for x RR in Equation 

291, and solve for the required L) . 

^ m 

D m MX RF 4X R )f^E_W-- + X R -X TD (296) 
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Equation 295 expresses the D that will allow t to increase to the 

m x 

value of t^. This value of D^, however, may be too large i.e., 
must be constrained as follows: 

Lo+D m <R L (296) 

where 

Rp = runway length 

Lp = required landing distance to stop aircraft 

If D m is too large to meet the constraint of Equation 296, the maximum 
value of D m is Rp - Lp. If this occurs, t^ will still be greater 
than t , then Sacrifice 4 is applied as a last resort. Sacrifice 4 
is applied by equating t^ and t h in Equations 291 and 293, by substituting 

X RF = x td + r l - l D 

as the limiting value for X R p in Equation 291, and solving for the 
increased h R p. 

• • 

^RF : Si L (h R F - h R)-^R (297) 


If a landing situation is such that Sacrifice 4 is necessary, the 
sacrifice is one of deciding to land the aircraft harder versus landing 
further down the runway and risk running off the runway end. 


The preceding sacrifices are for the case in which t h is greater 

than t . For the case in which t is greater than t. , Eouations 291 
x x n 

and 293 also suggest four possible sacrifices of the desired touchdown 

constraints to make t K = t : 

h x 


0 ) 

( 2 ) 

(3) 

( 4 ) 


Decrease x R p, that is, touchdown shorter 


Increase the landing speed, x 


RF 


Decrease h 


RF 


Make h R p less negative 
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Sacrifice 3 here is eliminated for the same reason as previously. 

Sacrifices 1 and 4 are applied in the same order as before except with 

different constraints. The constraint on D is that it cannot be 

m 

negative i.e., one does not want to touchdown prior to the effective 


runway beginning. 


In the event the required D m is negative, then 


Sacrifice 4 is applied and h DF is made less negative; however, h 


'RF 

cannot be positive. The equations for and h RR 


RF 

are the same as 


Equations 295 and 297, respectively. In the event h RR 


is reduced to 

zero, Sacrifice 2 (though undesirable) can be applied as a last resort. 
The equation for x Dr is obtained by equating t v to t h , substituting 

as follows: 


X RF " X TD and h RF 


Rp .. — — '■m— '-"■a - x — » h 

0, and solving for the required x 


RF 


X RF = h R 


(X rf-xr) 

(hRp-hR) 


-Xr 


(298) 


At this point, the desired touchdown conditions have been logically 
manipulated so as to occur at the same time within specific constraints 
of the runway and aircraft. These updated values cf x R p, x R p, h RF , and 
h R p (not changed) are then used with the present conditions X R , 

X R , h R , li R in Equations 290 and 292 to calculate the required 
accelerations A D and A. n to make the transition. 


b. Flare Dynamics 


, The preceding paragraph merely defines the required accelerations 

A xR and A^ R to perform the flare maneuver. The values of angle of 
attack, angle of sideslip, roll angle, and thrust to achieve the 
particular A^ r and A^ R are yet to be found. As indicated in the 
preceding paragraph, the flare maneuver will be done with wings level 
(desired roll angle = 0) and the aircraft crabbed (desired angle of 
sideslip = 0) into th r> o’.t • «»---d Assuming the aircraft yaw is 

S’ .a 11, A^ r and A^ R are t fr. • r 1 11 y in the aircraft pitch 
plane aid are predominantly ■ termined by angle of attack and thrust. 


The relationship between the scalar accelerations A^ r and A^ R and 
the independent variables a and T is as follows. Examine Figure 33, 
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Figure 33. Flare Forces 


velocity vector 
wind vector 
airspeed vector 
angle of attack 

elevation angle of V, relative to runway 

u 

elevation angle of V relative to earth 
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runway elevation angle 
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Since A xR and A^ R are defined relative to the runway, they are: 

A x r = j^Tcos(a+y 0 ) + Lsin(-y Q ) 

(299) 

- Ocos(-ya) -\A -psinERj /m 
At-,p = |^Tsin(a + 7 q) t Lcos (-^q ) 

(300) 

+ Dsin!-/□) — WfcosEpj /m 


For a given airspeed vector, L and D are predominantly quadratic 
polynomials in a (see Equations 279 and 281). Equations 299 and 300 
are therefore nonlinear functions of a and T whose solution is not 
immediately obvious. The observation that Equation 300 is predominantly 
a function of L and and is not appreciably affected by T and 0 
indicates that A^ )R is primarily dependent on a . This also indicates 
that the final value of A xR in Equation 299 is achieved through T. 

Both of these observations are helpful in constructing a numerical 
solution. Suffice it to say, that Equations 299 and 300 can be solved 
for the desired angle of attack, a^, and desired thrust, T rf , to 
achieve the required accelerations A xR and A hR . It should be noted 
that both a d and T^ have upper and lower bounds so that a situation 
can arise in which the requested accelerations are outside the capability 
of the aircraft. In such cases, the a d and T^ will be those that 
give the least vector error in the requested acceleration. A flow 
chart of the equations and logic which determine a d and T d for the 
flare maneuver is presented in Figure 34. A pit»~h rate calculation is 
made based on the desired angle of attack rate and the nominal rotation 
rate of the velocity vector. This pitch rate is called q d , the 
desired pitch rate in the flare. 


c. Hold - Decrab Maneuver 

As the aircraft approaches touchdown, it must be "decrabbed" 
(i.e., aligned with the runway center line) to allow the tires to roll 
and not skid. This is accomplished in the yaw autopilc'. (see yaw 
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autopilot development) by commanding the Euler yaw angle, ^-p, (see 
Appendix I) to be zero. The Euler roll angle, $ p , is also 
commanded zero for a "wings level" impact. More aircraft also have a 
"tail down" constraint which limits the maximum allowable angle of attack 
near the ground. 

The "hold-decrab" maneuver is staged on the mass center altitude 
above the runway. When this occurs, the desired angle of attack, 
a^, is held constant at the last value requested in the flare unless 
it must be limited by the tail-down constraint. The desired thrust, 

T d , is also held constant at the last value requested in the flare 
unless the "ki11-power" option (see throttle autopilot development) is 
exerc'sed. Most of the logic for the "hold-decrab" maneuver is directly 
in the autopilots. That logic which is necessary prior to autopilot 
entry is summarized in Figure 35. 


3. LANDING ROLLOUT 


The landing rollout begins at the instant any one of the tires 
touches the runway. The following maneuver will automatically occur 
in the landing rollout: control the aircraft yaw, <^ p , to be zero 
(i.e., keep the aircraft aligned with the runway); reduce roll control 
surface deflections to reutral position. The following events can be 
sequenced (on time after impact, t ) if so desired: 


( 1 ) 

( 2 ) 

(3) 

( 4 ) 

(5) 

final value 


Actuate ground spoilers (t p > t ) 

Kill power and reverse engines (t £ * ) 

Actuate drag chute (t p £ t ch ) 

Actuate tire braking (t r £ t^) 

Change elevator deflection from value at impact to a desired 
(in Pitch Autopilot, if t > t ^). 
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Most of the logic for the landing rollout maneuver is contained 
directly in the autopilots. That logic necessary prior to autopilot 
entry is shown in Figure 36. 

4. TAKEOFF ROLL 

The takeoff roll maneuver is simply one of rotating the aircraft in 
pitch to the desired takeoff angle of attack once the takeoff airspeed 
has been reached. Throttle setting is held constant through the roll 
at the takeoff value. The aircraft, yaw angle, is commanded zero 

in the yaw autopilot to keep the aircraft aligned with the runway, and 
the roll control surface deflections are kept in the neutral position. 
Elevator deflection is kept at a fixed value until the command to rotate 
the aircraft to the takeoff angle of attack is given. The takeoff roll 
maneuver is terminated on a particular mass center altitude, h s> above 
the runway. Most of the logic for the takeoff roll maneuver is contained 
directly in the autopilots. That logic which is necessary prior to 
autopilot entry is shown in Figure 37. 

5. PROBLEM PHASE LOGIC 

Given appropriate input, the phase logic determines which maneuver 
logic (i.e. glide slope, flare, hold-decrab, landing roll, takeoff 
roll) to use. The phase logic also determines where the problem is to 
terminate. Figure 38 shows the phase logic. 


The indicator ITO determines whether the problem is a takeoff or 
landing problem. If the altitude of the aircraft is above hp the 
aircraft is in the glide slope phase. When the altitude hp is reached, 
an option is provided to stop the program (see NF indicator) or go to 
the flare phase. The hold-decrab phase begins when the mass center 
altitude above the runway, h p , is less than or equal to h^ p + 8^ or 
if the distance down the runway X p has exceeded the desired value x^ p . 
The input value of h pp is normally h^, so that 8^ represents the 
vertical distance between the runway and tire bottom surface at which 
the hold-decrab phase begins. The KP indicator Is used to start the 
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Figure 38. Problem Phase Logic 





AFFDL-TR-71-155 
Part II 

program in the landing roll phase. If KP = 1, the time of impact, 

Tj, must be provided by initial data input. If KP * 1 the problem 
remains in the hold-decrab phase until any one of the tires receives a 
deflection, 8.. As soon as a tire deflection is received, the time 

of impact is stored in Tj and is never changed (even if the aircraft 
bounces on the runway). Immediately on impact, the NLRI indicator 
allows the program to be stopped or to continue with the landing roll 
phase. The landing roll phase can be terminated on a velocity, position, 
or time constraint, whichever occurs first. 


6. FAILURE LOGICS 

Prior to going to the autopilots, each maneuver logic checks for 
possible engine failures. In the event the aircraft is in the landing 
roll phase, a check is also made on possible brake failures. 

a. Engine Failure Logic 

Figure 39 is a flow diagram of the engine failure logic. The 
indicator array IC determines the failure (IC = 0) of up to four 
independent engines. While in the glide slope phase (i.e., IAP = 1), 
engine failures can be sequenced on two different altitudes h^ and h^. 
While in the flare or hold-decrab phases (i.e., IAP = 2,3), engine 
failures can be sequenced on two different runway altitudes h^ and h R £. 
While in the landing roll phase (i.e., IAP = 4) engine failures can 
be sequenced on two different times after impact t^ and t r p. While 
in the takeoff phase (i.e., IAP = 5, 6) engine failures can be 
sequenced on two different positions down the runway x R p^ and x R p^. 

The IC array is used in the throttle autopilot. 


b. Brake Failure Logic 

Figure 40 shows the brake failure logic. 


The Indicator array 


Ig determines the conditions of each gear brake. The Ig array can 
be changed from its initial array on two different times after impact 


t bkl and t bk 2. 


The Ig array Is used in the brake autopilot. 


This concludes the maneuver logic. 


180 





Figure 39. Engine Failure Logic 
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SECTION III 


t 


b 


i 




AUTOPILOTS 

The maneuver logic of Section II defined control of the aircraft by 
specifying desired values of angle of attack, angle of sideslip or yaw 
angle, roll angle, thrust, and status of the engines and brakes. It 
then becomes the function of the pitch, yaw, roll, throttle, and brake 
autopilots to achieve, within the capability of the aircraft, the 
desired values requested. Each autopilot development is discussed in 
turn. 


1. PITCH AUTOPILOT 

The basic function of the pitch autopilot is to control the aircraft 
angle of attack to the desired value requested by the maneuver logic. 

The aerodynamic surface for this control is the elevator on the horizontal 
stabilizer. Equation 38 expresses the major static moments acting about 
the body mass center pitch axis (excluding landing gear pitching 
moment, M ). 


[ C m 0 + C ma a+C ma 2|a| c-C N fp * 

c mg,8q +c m$ < , , |&q|&q + C mq 0 Sd l +M T +M M 


( 301 ) 


where 

^mo = pitching moment coefficient about the aerodynamic 
reference station at zero angle of attack and zero 
elevator deflection 


a = angle of attack 


C , C 

mo ma 


pitching moment coefficients for angle of attack 
dependence about aerodynamic reference station 


M = ground reaction moments in pitch (net moment) 
m ' 

C,. = none,a) force coefficient 

N 

AX - distance from aerodynamic reference station to 
mass center 
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= elevator pitching moment coefficients 
= elevator deflection 

= pitch reference length (usually mean aerodynamic 
chord) 


S = reference area 


q* = dynamic pressure 

Mj = engine thrust moments in pitch (net moment) 
r 

rriq = dynamic damping derivative in pitch 

= pitch rate based on d d and flare accelerations 
V, = airspeed 

G 


The nominal pitch trim requirements (i.e., the control surface 
deflection required to attain a particular desired angle of attack under 
equilibrium conditions) are obtained by requiring Equation 301 to zero. 
As mentioned in Section I, the aircraft is initially trimmed for the 
desired glide slope conditions. If this initial pitch trim is done by 
some surface other than the elevator, as with a moveable trim tab or 
flap, this is done through appropriate data input of C mQ . All pitch 
trim changes from the initial are assumed to be done by the elevator. 

The nominal elevator deflection, 8^, f° r a particular desired angle 
of attack, a^, is therefore obtained by solving Equation 301 set to 
zero. 


[c-rro + ^ma Q d + '“ma2| a d | a d f-N ^ 

+ Crn 8q® QN * Crn Sq l |^ q N + c mq q d ("fv^)] q Sd, +Mj +M m = 0 


( 302 ) 


Note that the equation is quadratic in 8^ and is also dependent 
on engine conditions. 
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Because there may be insufficient natural damping and because the 
a d command may constantly change (as in the flare), the 8 qN 
command to the elevator must be augmented to obtain good aircraft 
response and the desired angle of attack control. The pitch autopilot 
flow diagram is given in Figure 41. 


Three distinct routes are provided in the pitch autopilot, depending 

on the valu> of the indicator, IAP. When IAP = 5, this indicates a 

takeoff roll problem in which the takeoff airspeed has not been reached. 

For this case the desired elevator deflection, 8^, remains fixed at 

the input, 8^. When IAP = 4, this indicates a landing roll problem 

in which one may choose to change the impact value of 8 qd (i.e., 8 q j) 

to some final value 8^, at a rate 8p, and begin the change at time, 

t st , after impact. This impact maneuver can be used to change the 

horizontal stabilizer loads during the landing roll. The final route, 

IAP = 1, 2, 3, 6, indicates a maneuver in which a particular angle of 

attack, a d , is required. Here the pitch autopilot computes the 

angle of attack error <* e , the rate of change of error a Re , and 

includes rate feedback Rp fl , into the total angle of attack error 

signal, a g j. The allowed error in a g j is a fixed input, Aa a . 

If a T is within the allowed Aa , 8„j is set to the trim value 

61 3 qa 

8g N . If a gT exceeds the allowed error A<* a , 8 qd is set to the 
trim value 8 qN plus some overcontrol, which is determined by the 
product of a eT , and an input constant PS^. The desired elevator 
deflection, 8 qd , is finally limited by the aircraft constraints 8 qL 


(lower limit) and 8^ u (upper limit) 


The rate of feedback constant, 


R 


Fa 1 


and the overcontrol constant PS H> allow the pitch autopilot to 


be adjusted to a particular aircraft configuration. 


ALPDL prevents discontinuities in a d from entering the pitch 
autopilot thru a d - 8 q „ adds a ,: bang-bang" control capability to 
the pitch autopilot. 
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2. YAW AUTOPILOT 

The function of the yaw autopilot is to control either the sideslip 
angle of attack or the Euler yaw angle. In both cases the desired 
value of the angle is zero. The aerodynamic surface for this control 
is the rudder on the vertical stabilizer. Equation 303 expresses the 
major static moments acting about the body yaw axis. 


where 


[ c n/3/3 + C n £2 |/^j (3 c y dg +c ng r ^ r + 
Cn Sr*|^ r |® r ] q * Sdz + Nt 


(303) 



= sideslip angle of attack 

= yaw moment coefficients for dependence about 
aerodynamic reference station 

= body side force coefficient 

= yaw reference length (usually wing span) 

= rudder yaw moment coefficients 

= rudder deflection 

= engine thrust yaw moment 


The body side force coefficient is also predominantly expressed by 

Cy - Cy/3)3 +C y/Q 2|/3| /3 (304) 

where C ~ and C _2 are the side force coefficients for Q. 

y(B y/3 H 

Substituting- Equation 304 into Equation 303 yields 


[<c n /3-c y ^-^-)^+(c rV3 «-c yj9 *^-)|^|a + 

^ n 8r^ r + ^ n Sr |S f !& r ] q ^*T 


(305) 
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When Equation 305 is set to zero, it can be used tc solve for the 
nominal rudder deflection, 8 to trim the aircraft at a particular 
desired sideslip angle of attack, (3^. Note that Equation 305 is 
quadratic in 8, and is dependent on engine conditions. Examine 
Figure 42, which is the flow chart for the yaw autopilot. 


Two distinct routes exist in the yaw autopilot, depending on the 
value of the indicator, IAP. If IAP =1,2 (i.e., glide slope or 

flare) the desired value of jS is zero. The solution to Equation 305 

for the nominal rudder deflection, 8 will also be zero, except in 

the case where the engine thrust yaw moment is nonzero. Rate feedback, 
Rpg, and overcontrol, PS R , are provided for the same reasons discussed 
in the pitch autopilot. If IAP = 3, 4, 5, 6, the aircraft may have to 
operate at a non-zero steady sideslip (only for a cross wind situation) 
to make the Euler yaw ana^, typ, zero. The nominal rudder deflection, 
S rN , for this case is obtained from Equation 305 by substituting 
Equation 306 for (3 Equation 306 expresses approximately 

• 

B(j = sin 1 - (306) 


the steady sideslip /9 at which the aircraft must operate in order to 
align with the runway). Rate feedback, Rp^, and overcontrol, PS^,, 
are also provided in this route of the yaw autopilot. The desired 
rudder deflection, 8^, is finally limited by the aircraft constraints 
8 rL (lower limit) and 8 pu (upper limit). As with the pitch auto¬ 
pilot, the yaw autopilot can also be adapted to a particular aircraft 

Ipg and Rp^, 


configuration through the rate feedback constants, R™ and R p .,., ar.d 


the overcontrol constants, PS R and P3^ 


3. ROIL AUTOPILOT 

The function of the roll autopilot is to control the Euler roll 
angle to the desired value requested by the maneuver logic. The 
aerodynamic suiface for this control is the aileron. The roll autopilot 
differs from the pitch and yaw autopilots in that it assumes there is no 
major static moment developed about the roll axes for the maneuvers 
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requested. As such, the nominal aileron deflection, S p ^, to trim 
the aircraft about the roll axis is always zero. Other than this trim 
requirement, the roll autopilot is similar in concept to the pitch and 
yaw autopilot. The roll autopilot flow diagram is given in Figure 43. 


Two distinct routes exists in the roll autopilot, depending on the 
values of the indicator, IAP. If IAP = 4, 5, 6, the aircraft is on 
or has just left the runway and the roll attitude is determined primarily 
by the landing gear runway interface. As such, the desired value of 
aileron deflection, 8 pd> is zero. If IAP - 1, 2, 3, the aircraft 
is in flight and a particular desired roll angle, is commanded by 

the maneuver logic. Rate feedback, R p ^ , and over control, PS^, 
are provided, as was done in the pitch and yaw autopilots. The final 
desired aileron deflection, 8 pd , is limited by the aircraft constraints 
8 pL (lower limit) and 8 p y (upper limit). As with the pitch and 
yaw autopilots, the roll autopilot can be adapted to a particular 
aircraft configuration by appropriate selection of the constants R p ^ 
and P$^. 


4. THROTTLE AUTOPILOT 

The basic function of the throttle autopilot is to achieve the 
desired thrust, T^, requested by the maneuver logic. It performs 
this function by commanding desired values of the throttle settings. 

The throttle settings must be arrived at under any combination of the 
following conditions: 1, 2, 3, or 4 engine aircraft, all engine failure 
combinations practicable, engine reverse under selected throttle 
constraints, capability to let some engines carry more load than others. 


The throttle autopilot is divided into four logics = one for each 
of the engine-aircraft combinations (see Figure 44). Each of these 
logics is built up from two other basic engine logics called, "comnon 
engine logic" and "common two-engine logic." 
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a. Single-Engine Logic (SEL) 

The single-engine logic determines the desired throttle setting, 

(1), for a single-engine aircraft. A "fixed throttle" mode and a 
failure mode is provided. Figure 45 is a flow diagram of the SEL, In 
this logic, a check is first made to see if the engine is in the "fixed 
throttle" mode. This is done by examining the value of the first 
element in the fixed throttle indicator array, TF(I) - the TF(I) array 
is read on initial data input. If TF(1) = 1, this indicates a fixed 
throttle mode for which (1) is given the value N^p(l) (the N^p array 

is read in on input). A landing reverse check is then made to determine 
if this is desired. This is followed by a failure check (i.e. IC(1) = 0, 
as shown in Figure 39). If the engine has failed, N d (1) is set to 
zero; otherwise N d (1) remains unchanged. If the engine is not in a 
fixed tnrottle mode (i.e., TF(1) = 1), a failure check is again made. 

If thete is no failure, the data set up (use first element of input 
arrays) for the CEL is done. Entry into the CEL is made and the 
desired throttle setting, N^O) is determined. Note that the SEL is 
predominantly the CEL except that the fixed and failure modes are 
determine outside the CEL. 

b. Common Engine Logic (CEL) 

The CEL is applicable to any "variable throttle" engine on the 
aircraft. Given a set of variables, the CEL assumes the engine has 
not failed and finds a desired throttle setting. The set of variables 
required by the CEL is: 

IAP = problem phase indicator; available from maneuve** logic 

N c = actual throttle setting; available from calling program. 

T c * desired thrust; available from calling program. 

ILR = reverse engine signal; available from maneuver logic 

IRC - indicator for engine reverse capability; available from 
calling program 

N BC = bottle setting above which reverse should not be 

calling program actuated; available from calling program. 
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N. Rf = throttle setting for landing reverse; available from 
calling program. 

N 

TOC = throttle setting for takeoff; available from calling 
program. 

KT = kill engine indicator; available from calling program. 
Figure 46 is a flow diagram of the CEL. 

The IlR indicator is initially input to the program as a value other 
than 1. The landing reverse engine signal (i.e., ILR =1) is niven 
in the landing rollout phase of the maneuver logic (see Figure 36;. 

If ILR =1, a check is made on IRC (a value of 1 indicates that the 
engine has a reverse capability). If IRC * 1, the desired throttle 

setting, N d( ,, is set at 1.0 (this value of throttle setting is 
assigned to forward idle - see Section VI) which is the lowest possible 
forward throttle setting. If IRC = 1, a check is made to see if the 
actual throttle setting, N c , is below the value required, N 0C , to 
engage the reverse. If the check is not true, the desired throttle 
setting, N dC , is made 1 so that will reduce below N B( .. When 
N c < Ng C , the reverse is engaged by requesting N dr to equal the 
throttle setting for landing reverse, N LR( .. 

If the landing reverse signal has not been given (i.e., ILR #1), 
the.CEL has three routes depending on the value of the indicator, IAP. 
If IAP = 5, 6, this indicates a takeoff problem and N dC is set to the 
takeoff throttle setting, Njq C> If IAP = 3, 4 (i.e., hold-decrab 

or landing roll maneuver) an option is provided to kill power. This 
is done during data Input by making the kill engine indicator KE(I) = 1 
If KT y 1, t.he problem goes to the IAP - 1, 2 route where thrust will 
be maintained during the hold-decrab maneuver and will be maintained 
during the landing rollout maneuver until the landing reverse signal is 
given. If IAP =1,2 (glide slope or flare), a route is desired to 
find the throttle setting for a particular thrust, T . Here an 
interpolation is done in the thrust tables from N = -2.0 (full reverse) 


V 
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Figure 46. Common-Engine Logic 
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to N = -i.O (idle reverse) and from N = +1.0 (idle forward) to 
N = +2.0 (max forward) to find the throttle setting, N, that gives 
T c or comes closest to T . If the final interpolated value of N is 
negative (reverse thrust requested), the same throttle constraint is 
applied as was done for the landing reverse. It should be noted that 
the throttle constraint on reverse thrust applies only when engaging 
reverse and does not limit disengagement (N^ can always have a 
positive value). 

c. Function ENGREV 

As originally designed, any engine in a "fixed throttle" mode 
could not reverse when the reverse signal (ILR = 1) was given. A 

function called ENGREV was added, which allows a "fixed throttle" 
engine to change to reverse landing provided ILR has a value of one, 
the engine has a reverse capability (if not, forward idle is commanded), 
and the throttle constraint for reverse is met. This addition was 
necessary since the "fixed throttle" engines do not go through the 
"common engine logic" (CEL). The function performs essentially the same 
flow as that portion of CEL which occurs if ILR =1, as shown in Figure 47. 



Figure 47, Function ENGREV 
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d. Two-Engine Logic (2EL) 

The two-engine logic determines the desired throttle settings 
N d (l) and N d (2) for a two-engine aircraft. One fixed throttle mode 
and all possible failure inodes are provided. Figure 48 is a flow 
diagram of the 2EL. 


First a check is made on the one fixed throttle mode allowed 
(i.e., if one engine is fixed the other is assumed fixed also). If 
the fixed throttle mode is indicated, a failure check follows and 
N d (l) and N d (2) are de.ermined in a manner similar to that done in the 
SEL. If the fixed throttle mode is not indicated, the data set up for 
entry into the CTEL is done. Here the T dx of the two engines is the 
desired thrust, T^, which comes from the maneuver logic. The data 
set up for engine A is the appropriate data from tne arrays for engine 1. 

The data set up for engine B is the appropriate data from the arrays 
for engine 2. The thrust fractions and Kg are determined by the 

constant k ^ 2 ) 121 ( see s y ,T * 0 ^ s f° r meaning of numerical subscripts). 

Note that the sum of K^ and Kg is numerically 1. The output of the 
CTEL is the desired throttle settings N dA and N dg . Note that except 
for the fixed throttle mode, the 2EL is predominantly the CTEL. 


e. Common Two-Engi n e Logic (CTEL) 

The CTEL determines the desired throttle settings N dA and N dg 
for two "variable throttle" engines, A and B, which can have any possible 
failure mode (i.e., A and B failed, only A failed, only B failed, 
no failures). The set of variables required for the CTEL is much the 
same as that required for the CEL except that engine A is distinguished 
from B by a suffix in the variable name. The variables required for 
the CTEL are: 


'dx 
ICA, ICB 


K A* k b 
IRA, IRB 


total thrust required of the two engines 
failure indicators for engines A and B 
fractions of T d for each engines A and B 
reverse capability indicators for engines A and B 
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= actual throttle setting for engines A and B 


^BA’ N BB = throttle settings above which reverse should not be 
actuated for engines A and B, 


N LRA’ n lrb 


throttle settinas for landing reverse of engines A 
and B 


= takeoff throttle settings for engines A and B 
KEA, KEB =■ kill engine indicator for engines A and B 


Figure 49 is a flow diagram of the CTEL. 


First a check is made to see if both engines have failed (i.e., 

T CA = 0 and ICB =0). If both have failed, both and are set 

to zero (this is used for engine failure in the thrust table look=up 
which is stored under the N = 0. data - see Section VI). If both engines 
have not failed, a failu"e check of engine A is done. If engine A has 

failed, is set to zero and the thrust table lookup is entered at 

N = 0 to determine the desired thrust, T^, of engine A (in this case 
is actually the failure thrust). The required thrust of engine B, 
T dB’ t0Lmd subtracting T dA from T dx (engine B carries the full 
thrust load). The input variables for the CEL (IRC, N c . . . etc.) 
are set up using the data for engine B af n-, which the CEL is entered. 
The output throttle setting of the CEL, N d( ,, becomes the desired 
throttle setting, N^g, of engine B. If engine A has rot failed, a 
failure check on engine B is done. If engine B has failed, engine A 
is assumed to carry the full thrust load. The determination of N dA 
when engine B has failed is similar to the previous case where N d g was 
determined when engine A failed. If engine B has not failed, this 
indicates that both engines are working. In this case, the thrust 


f raction 
The CEL, 


and Kg determine what portion of T dx each engine carries, 
with appropriate engine input data, is used to determine the 


desired throttle settings and N d g of each engine. 
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F. Three-Engine Logic (3EL) 

The three-engine logic determines the desired throttle settings 
N d (l), N d ( 2 ) and (3) for a three-engine aircraft. Three-throttle 

fix modes and all possible failure modes are provided. Figure 50 is a 
flow diagram of the fixed-throttle mode section of the 3EL, 


k 


k 

\ v 


The first fixed-throttle mode allowed, that is checked, is the mode 
where all engines are fixed TF(1) = 1 and TF(2) = 1 and TF(3) = 1. Here 

the desired throttle settings are arrived at in a fashion similar to 
that done in the SEL and the 2EL. If this fixed-throttle mode is not 
indicated (note all this takes is For one of the three TF indicators 
not to be equal to 1), the next allowed fixed-throttle mode is checked. 
This mode is one where the outboard engines (engines 1 and 3) are 
assumed fixed, and the center engine (engine 2) is assumed variable. 
Both of the two fixed engines are allowed to fail. The thrust table 
look-up is then called to evaluate the desired thrusts of engines 1 and 
3, be they failed or not. Failure of the variable engine (engine 2) 

is also allowed. If engine 2 is not failed, it is assumed to take the 

remaining thrust load required to meet the T d requirement of the maneuver 
logic. Data set up for the CEL is made using array data for engine 2. 
Entry is made to the CEL and N d (2) is determined. The final fixed 
throttle mode allowed is where the center engine (engine 2) is fixed 
arid the two outboard engines (1 and 3) are assumed variable. Again, 
engine 2 is allowed to fail and the thrust table look-up is called to 
get the desired thrust of engine 2, (2) „ be the engine failed or not. 


Data set up is then made to enter the CTEL, where T dx - . d - , d 


= T d - V 2 >- 

Engine 1 data is assigned to engine A and eng>ne 3 data is assigned to 
engine B. The thrust fractions k^ and kg are determined by the constant 
k(3) 131 * ^ ot:e a 9 a ' n that plus Kg numerica1 'y 1 • The CTEL is 

entered and the desired throttle settings N d (l) and N d (3) are determined. 
If none of the three fixed throttle modes is indicated, all engines are 
assumed variable throttle. For this section of the 3EL, the flow 
diagram is given in Figure 51. 
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This section of the 3EL allows for all possible failure combinations 
of the three engines. The logic is as follows: if engine 1 has failed, 
call the CTEL for engines 2 and 3; if engine 2 has failed, call the CTEL 
for engines 1 and 3; if engine 3 has failed, call the CTEL for engines 
1 and 2. CTEL, in turn, checks for failure in the remaining two 

engines. Note that in any engine failure combination, the required 

thrust load is carried by the remaining working engines. If no engines 
have failed, the CEL is called for each of the three engines separately. 

g. Four-Engine Logic (4EL) 

The four-engine logic determines the desired throttle settings 
N d (l), N d (2), N £j (3) and N rf (4) for a four-engine aircraft. Three 
throttle fix modes and all possible symmetric failure modes are provided 
(failure of engines 1 and 2 is similar to failure of engines 3 and 4, 
etc., and therefore only one mode is simulated). Figure 52 is a flow 
diagram of the fixed-throttle mode section of the 4EL. 

The first fixed throttle mode is for all four engines fixed. The 
logic here is much the same as was done in the 2EL and the 3EL. The 
second fixed-throttle mode is for the two inboard engines (i.e, engines 
?. and 3) fixed and the outboard engines (1 and 4) variable. The 
desired throt + le settings, N^(l) and (4), for the two variable 
engines are obtained by calling the CTEL for engines 1 ano 4. The third 
and last fixed-throttle mode is for the two outboard engines fixed and 
the two inboard engines variable. The logic is similar to the previous 
case. If none of the three fixed-throttle modes is indicated, all 
four engines are assumed variable. Figure 53 is a flow diagram of this 
section of the 4EL. 

The symmetric failure modes allowed are as follows: four-engine- 
failure, 1-2-3-4; three-engl.ie-failure, 1-2-3 and 1-2-4; two-engine 
failure 1-2, 1-3, 1-4, and 2-3; and single-engine fai’ure 1 and 2. 
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Figure 63. Four-Engine Logic 
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Figure 53. Four Engine Logic (Font) 
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The four-engine failure mode is checked first. This is followed 
by a check on the three-engine failure mode 1-2-3. Ir. this mode engine 
4 is assumed to carry the thrust load. If the three-engine failure 
mode 1-2-3 is not indicated, a check is made on the last three engine 
failure mode 1-2-4, where engine 3 is assumed to carry the thrust load. 

If failure mode 1-2-4 is not indicated, the two-engine failure mode 1-2 is 
checked, where engines 3 and 4 carry the thrust load. This is followed 
by successive checks of the remaining two-engine failure modes 1-3, 1-4, 
and 2-3. In each case, the remaining engines carry the required 
thrust load. If none of the two-engine failure modes are indicated, 
the two single-engine failure modes (engine 1 and/or engine 2) are 
checked. If a single-engine failure mode is indicated, the remaining 
three engines carry the thrust load. If no single-engine failures are 
indicated, all engines are working, and the desired throttle setting of 
each engine is obtained by calling the CEL for each of the four engines. 
This completes the 4EL and the throttle autopilot discussion. 


5. BRAKING AUTOPILOT 

The function of the braking autopilot is t.o control the braking 
moments applied to each landing gear on the aircraft. The braking 
autopilot provides the following four options for each landing gear: 

(1) No braking (i.e., braking moment zero) 

(2) Constant braking moment (comparable to constant braking pressure) 

(3) LocKed wheel (i.e., wheel angular velocity zero) 

(4) Controlled braking 

The wheel equation of motion is documented in Appendix II. This 
equation. Equation 23C is repeated here for convenience. 


Wj [ 

m AI " M B1 * "i -| “>TI 


( 307 ) 


n 
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where 

M^. = ground reaction moment along axle of gear 

i. L. 

Mg. = braking moment applied along axle of i gear 

u 1 (only + values allowed) 

= wheel rotational velocity on i^ gear 

= number of tires on i tn gear axle 

I. = moment of inertia of a tire wheel, and anything else 

1 constrained to rotate with that tire about the i th gear 

axle. 

• X L. 

Wj^ = wheel angular acceleration on i cn gear axle 

The function of the braking autopilot, therefore, becomes one of 
specifying the value of Mg^ in Equation 307. Figure 54 is a flow 
diagram of the braking autopilot logic. 

The braking autopilot logic is contained in a "do loop, 1 ' which is 
repeated I times. The indicator, I, is the number of gears on the 
aircraft. When the "do loop" is finished, the M g .. array contains 
the braking moment for each gear. The first check inside the do loop 
is on the value of the brake actuate indicator, IBS. IBS is normally 
input at some other value that 1 so that M gi remains at the initial 
values read into the program (the initial values of M g ^ are usually 
zero). The braking signal (IBS =1) is given in the landing rollout 
phase of the maneuver locic (see Figure 36). When the brake signal is 
given, the brake condition indicator, (see Figure 40) is examined 

to determine which of the four brake options is to be applied. If the 
value of I gi is 0, brake option 1 is exercised and is set to 0. 

If Ig^ has the value 2, brake option 2 is exercised and M g ^ is set 
to the constant value M 0 ~., which is read in on input. If the value 

DL 1 

of Ig^ is -1, brake option 3 is exercised and Wy. Is made 0 to 
simulate brake locking If the value of I g . is 1, controlled braking 
option 4 is exercised. 
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Controlled braking option 4 controls the wheel angular speed, , 
to that value which results in a particular desired value of "percent 
skid," Pp. Since the coefficient of friction between the tire and 
runway is predominantly a function of "percent skid," controlling 
"percent skid" actually controls the braking forces applied to the 
aircraft. "Percent skid" is defined as follows: 


P 


st 



where 

P $i = percent skid of i wheel 

V-^ = tire footprint velocity of wheel 

V . = axle velocity of i th wheel 

3 i 


(308) 


The tire footprint velocity can be written as: 


V Ti » V oi + w Ti (r oi - 8,1 (309) 

where 

r . - undeflected radius of tire 

01 

S i = tire deflection 


Equation 308 can therefore be written as 


■ti 


V ai + w Ti U oi" S i ) 


r oi 


(310) 


Note that, when w is 0 (full skid), the value of is 1; and 
when the tire is fully rolling (Vj, = 0), the value of P $i is zero. 
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Let w jpj b® the required value of ou that 
percent skid P^. Since V . is always positive. 


results in the desired 
Equation 310 yields 


<'-p d ) 

“TRi ’ - V ai „ oi . 8l , 


(31 I) 


where 


'at 




611 


R DXGk + 


613 


r dzgk’ 


R DYGk 


(312) 


Assuming the tire deflection reaches a nominal value, the value of 
reduces because the axle speed V . reduces during the landing 


roll 


The rate at which 


TRi 


should reduce to maintain P~ is 


D 


therefore 


ai 


TRi 


<< -Ppl 

( r 0 j - Si ) 


(313) 


The variable u is the body longitudinal acceleration and is approximately 
Fx/m (see Appendix I) where 

F^ = sum of forces acting in x direction 

in = aircraft mass 


(Note that F includes not only ground forces, but also forces of 
drag, reverse thrust, drag chutes, etc. Equation 313, therefore, 
becomes 


u) 


TRi 


F x (l~Pp> 

W I'oi-V 


(314) 


(See Figure 54 and the option 4 branch). 
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First the desired value of wheel angular speed, is calculated 

using F^uation 311. The required rate of change of cuyp^» - ’ s 

also calculated using Equation 314. This is followed by a calculation 
of the error, in wheel angular speed. If the error, w.^-, 

is less than a certain fraction (Aw is the allowed fractional error 
in w-j-p.j) of w jR.j> Then the braking moment, , should be such 
that the wheel angular acceleration u> is maintained (this helps 
to avoid further error development). The required value of is 
obtained by rearranging Equation 307 and substituting Wyp-j Tor u^... 


M 


Bi 


I^Tijr 

w Ti L 


m ai" 


w TRi Ij n 


.] 


(315) 


If the allowed error is exceeded and the error is positive, this means 
the value of needs to be more negative (i.e., mo^e rolling and 

less braking). This requires a negative The absolute magnitude 

of this control acceleration is picked a constant value w c (chosen at 
data input). The required braking moment in this case is found by 
substituting -w c into Equation 307 for and solving for the 

required braking moment. Note in this case 

M Si'^ [ M Ai + "cli w i] (3,6) 

will be loss than (J . If the allowed error is exceeded and 
the error is negative, this means the value of needs to be less 

negative (i.e., more braking). This requires a positive Wj. and, 

as previously, will have a control magnitude of u» c . Substituting 
+ oj r into Equation 307 for yields the required braking moment. 

m BI * 15'TI 

Note in this ca^e will exceed || . Finally, the requested 

value of Mg. is checked against ’he braking moment limits M^.j (lower 
limit - not. less than zero) and M,^. (upper limit). This completes 
the braking autopilot discussion. 
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SECTION IV 

CONTROL VARIABLE RESPONSE 

The main aircraft control variables {i.e., elevator deflection, 
rudder deflection, aileron deflection, and engine throttle) do not in 
reality respond instantaneously to the desired values requested by the 
autopilots. To get a first order approximation of the effects of 
control surface and engine lags on landing performance, a "constant 
rate" control variable response is built into the autopilot simulation. 
Figure 55 is a flow chart of the control variable response logic. This 
logic allows the control variables to move at a fixed rate (rate depends 
on initial data input) toward the desired values requested by the 
autopilots. No overshoot is allowed; that is, as the actual value of 
the control variable approaches the desired value, the response logic 
■locks the actual value to the desired value until a situation arises 
where the rate of change of the desired value exceeds the rate at which 
the control variable can respond. In this way, the first order effects 
of control variable lags are simulated. Note that a system with 
essentially instantaneous response can be built by input of large 
numbers for the control variable rates. 
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SECTION V 
SDF-2 CHANGES 

I. MULTIPLE ENGINE CHANGES 

This section was written in Appendix I. It is repeated here for 
emphasis. 

As originally designed, SDF-2 only allowed for the simulation of 
one airbreathing engine. The thrust table look-up routine, TFFS, also 
included tne effects of altitude on thrust and determined the fuel rate, 
both of which are no longer needed (see Assumption in Section I). 

Some provision must also be made to determine thrust for engine 
reversing, thrust for engine failure, and the net thrust forces and 
moments acting on the multiple engine aircraft. We begin by discussing 
the thrust table lookup. 

The thrust table is now a function of throttle setting and Mach 
number alone. Because of the throttle autopilot logic, the ranges for 
the throttle setting, N, must be as follows: 

N = -2. means full reverse 
N = -1. means idle reverse 
N = 0. means engine failure 
N ~ +1. means idle forward 
N - +2. means full forward 

Negative values of thrust are stored for the negative throttle settings. 

In this manner reverse thrust is simulated. The zero value of N is used 
for the data simulating engine failure (note that failure thrust may be 
zero or some negative value, depending on Mach number). The actual 
reversing is achieved in the throttle autopilot by requesting a negative 
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value of desired throttle setting, N^... The change of sign on N^ 
is sensed in the control response (see Figure 55) and this changes 
the sign of the actual used in the thrust table. 

Every time thrust is needed by the main program, the thrust table is 
entered IN times - once for each engine - and the actual thrust array 
T (IN) is obtained depending on the values in the actual throttle setting 
array N (IN). The engine thrust vector is assumed parallel to the 
longitudinal body axis (this is also assumed in the autopilot equations) 
and therefore causes no roll moments. The'engine position arrays, 
z N (IN) and y^ (IN), along with the actual thrust array T (IN) 
determine the engine pitch moment array Mj-(IN) and the engine yaw 
moment array Ny(IN). These arrays are then used to obtain the net 
longitudinal thrust, T , (note T and T are zero by assumption) the 

X y Z 

net engine pitch moment. My, and the net engine yaw moment, Ny. The 
net value T x , My, and Ny are then used in SDF-2 and the calculation 
proceeds as normal. 


2. AUXILIARY COMPUTATION 


Several state variables used by the autopilot must be defined in 
terms of the state variables existing in SDF-2. The variables that 
need to be defined are as follows: x R , y R , z R , x R , y R , z R , h R , ti R , 
t//p, and <$> p Their definition in terms of SDF-2 variables is as 
fellows: 



Gil 

0 

631 

Gil 

0 

G3I 


0 R Gj 3 
I 0 

0 R G33 

0 R GI3 
I 0 

R G33 


“ R ?r1 


*3 



(318) 


( 319 ) 
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h R * " Z R 


(320) 


h 


R “ 



(321) 



(322) 


(323) 


The variables x p , x p , h p , and are used in the flare maneuver logic 

(see Section II). The RG matrix (see Equation 111 of Appendix II) is 

a trains format!on from the earth coordinate system to the elevated runway 

• * 

coordinate system. The Euler angular races p and are used for 
the yaw and roll autopilots, respectively. They are developed for a 
pitch-yaw-roll sequence in terms of the direction cosines and direction 
cosine rates which are available in SDF-2. 
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SECTION VI 


PROBLEM ORGANIZATION 


Figure 1 shows the basic autopilot interface with SDF-2. The 
detailed autopilot interface and the logic organization is shown in 
Figure 56. A comparison of Figures 26 and 56 shows that the autopilot 
begins with the auxiliary computation and ends with the control variable 
response. Several changes were also required internal to SDF-2. 


This concludes the formulation of the equations and logic for the 


autopilot modification to SDF-2. 
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